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Many state-of-the-art QCD calculations for multileg processes use helicity amplitudes as their 
fundamental ingredients. We construct a simple and easy-to-use helicity operator basis in soft- 
collinear effective theory (SCET), for which the hard Wilson coefficients from matching QCD onto 
SCET are directly given in terms of color-ordered helicity amplitudes. Using this basis allows 
one to seamlessly combine fixed-order helicity amplitudes at any order they are known with a 
resummation of higher-order logarithmic corrections. In particular, the virtual loop amplitudes can 
be employed in factorization theorems to make predictions for exclusive jet cross sections without 
the use of numerical subtraction schemes to handle real-virtual infrared cancellations. We also 
discuss matching onto SCET in renormalization schemes with helicities in 4- and d-dimensions. 

To demonstrate that our helicity operator basis is easy to use, we provide an explicit construction 
of the operator basis, as well as results for the hard matching coefficients, for pp —>• H + 0,1, 2 
jets, pp —>• W/Z/ 7 + 0,1, 2 jets, and pp —> 2, 3 jets. These operator bases are completely crossing 
symmetric, so the results can easily be applied to processes with e + e“ and e~p collisions. 
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I. INTRODUCTION 

The production of hadronic jets is one of the most ba¬ 
sic processes at particle colliders. Processes including a 
vector boson (W, Z, 7 ) or Higgs boson together with jets 
provide probes of the Standard Model (SM), and are also 
dominant backgrounds for many new-physics searches. 
Optimizing the precision and discovery potential of these 
channels requires accurate predictions of the SM back¬ 
grounds. Furthermore, the growth of the jet substructure 
field has sparked a renewed interest in the study of jets 
themselves, both for an improved understanding of QCD, 
and for applications to identify boosted heavy objects in 
and beyond the SM. 

Precise predictions for jet production require pertur¬ 
bative calculations including both fixed-order corrections 
as well as logarithmic resummation. QCD corrections to 
processes with jets are typically enhanced due to phase 
space restrictions. Such restrictions often introduce sen¬ 
sitivity to low momentum scales, p, of order a few tens 
of GeV, in addition to the hard scale, Q , which is of 
order the partonic center-of-mass energy. In this case, 
the perturbative series contains large double logarithms 
a" ln m (p/Q) with m < 2 n. To obtain the best possi¬ 
ble perturbative predictions, these logarithms should be 
resummed to all orders in a s . 

There has been tremendous progress in the calculation 
of fixed-order perturbative amplitudes in QCD using the 
spinor helicity formalism , color ordering techniques 
[I|-Q and unitarity based methods H, Hcj. NLO predic¬ 
tions are now available for a large number of high mul¬ 
tiplicity final states, includi ng pp —» V+ up to 5 jets 
mm, V V 7 ~ ^ U P 5 jets [UEil], and pp —> H+ up 
to 3 jets mm, an< ^ there are many efforts (39l - l56| to 
fully automatize the computation of one-loop corrections 
to generic helicity amplitudes. 

For high-multiplicity jet events, the resummation of 
large logarithms is typically achieved with parton shower 
Monte Carlo programs. Here, the hard process enters 
through tree-level (and also one-loop) matrix elements 
and the QCD corrections due to final-state and initial- 
state radiation are described by the parton shower. The 
parton shower resums logarithms at the leading logarith¬ 
mic (LL) accuracy, with some subleading improvements, 
but it is difficult to reliably assess and systematically im¬ 
prove its logarithmic accuracy. 

The approach we will take in this paper is to match 
onto soft-collinear effective theory (SCET) {57| - f60j |. the 
effective theory describing the soft and collinear limits of 
QCD. In SCET, the QCD corrections at the hard scale 
are captured by process-dependent Wilson coefficients. 
The low-energy QCD dynamics does not depend on the 
details of the hard scattering (other than the underlying 
Born kinematics), similar to the parton shower picture. 
Resummation in SCET is achieved analytically through 
renormalization group evolution (RGE) in the effective 
theory, allowing one to systematically improve the loga¬ 
rithmic accuracy and assess the associated perturbative 


uncertainties. For example, for dijet event shape vari¬ 
ables in e + e _ collisions, SCET has enabled resummation 
to N 3 LL accuracy and global fits for a s {mz) [ 6 ll - [ 66 j . The 
analytic higher-order resummation can also be used to 
improve the Monte Carlo parton-shower description m- 
[(39). Furthermore, SCET allows for the direct calculation 
of exclusive jet cross sections, eliminating the need for 
numerical subtraction schemes for real emissions up to 
power corrections. 

An important prerequisite for employing SCET is to 
obtain the hard matching coefficients, which are ex¬ 
tracted from the fixed-order QCD amplitudes. The 
matching for V + 2 parton and H + 2 parton processes is 
well known from the QCD quark and gluon form factors, 
and is known to three loops [H, Ell 01 • The matching 
for V + 3 partons |72| - [73| . and H + 3 partons [76i - [79j |. 
has been performed at both NLO and NNLO. Partonic 
processes with four external quarks have been studied in 
SCET in Refs. [80l487j| . and the matching for all massless 
2 —> 2 processes has been obtained at NLO in Ref. (88j | 
and recently at NNLO in Ref. [h!)| . 

For high-multiplicity processes, the usual approach to 
constructing an operator basis with explicit Lorentz in¬ 
dices and gamma matrices is laborious. In this paper, we 
introduce a convenient formalism, based on helicity op¬ 
erators, which allows for a seamless matching for higher 
multiplicity processes onto SCET. A first look at the for¬ 
malism discussed here was already given in Ref. !)() . In¬ 
deed, results for helicity amplitudes are already employed 
in the SCET matching calculations mentioned above, 
though without the construction of corresponding SCET 
operators. 

In the spinor helicity formalism, the individual helicity 
amplitudes (i.e. the amplitudes for given fixed external 
helicities) are calculated, as opposed to calculating the 
amplitude for arbitrary external spins in one step and 
then summing over all spins at the end. One advantage 
is that the individual helicity amplitudes typically yield 
more compact expressions. And since they correspond to 
distinct external states, they can be squared and summed 
at the end. Helicity amplitudes remove the large redun¬ 
dancies in the usual description of (external) gauge fields, 
allowing for much simplified calculations particularly for 
amplitudes with many external gluons. 

As we will see, this helicity-based approach is also 
advantageous in SCET. In SCET, as we will review in 
Sec. IIIB1 collinear fields carry label directions corre¬ 
sponding to the directions of jets in the process, which 
provide natural lightlike vectors with which to define 
fields of definite helicity. As we will demonstrate, the 
construction of an appropriate operator basis becomes 
simple when using operators built out of fields with defi¬ 
nite helicity. Furthermore, using such a helicity operator 
basis greatly facilitates the matching of QCD onto SCET, 
because one can directly utilize the known QCD helicity 
amplitudes for the matching. Together, this substantially 
simplifies the study of high-multiplicity jet processes with 
SCET. 
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A. Overview 


Consider a process with N final-state jets and L lep¬ 
tons, photons, or other nonstrongly interacting particles, 
with the underlying hard Born process 

Ka(q a ) Kb{qb) -A Ki(gi) • • • K N+L (q N+L ), (1) 


where K a ,b denote the colliding partons, and Ki denote 
the outgoing quarks, gluons, leptons, and other particles 
with momenta qi. The incoming partons are along the 
beam directions, q£ b = x a ,bPa b i where x a , b are the mo¬ 
mentum fractions and P^ b the (anti)proton momenta. 
For definiteness, we consider two colliding partons, but 
our discussion of the matching will be completely cross¬ 
ing symmetric, so it applies equally well to ep and ee 
collisions. 

In SCET, the active-parton exclusive jet cross section 
corresponding to Eq. o can be proven to factorize for 
a variety of jet resolution variables!]] The factorized ex¬ 
pression for the exclusive jet cross section can be written 
schematically in the form 


d(J = / dXadXbd ® N + L ( qa+qb ’ qi '--^ M ({ qi ^ ( 1 2 ) 

x y]tr \H K ({qi})S K ] ® B Ka B Kb J K 


Here, d$Ar + i(g a + q b \qi,...) denotes the Lorentz- 
invariant phase space for the Born process in Eq. 0, 
and M ({</,;}) denotes the measurement made on the hard 
momenta of the jets (which in the factorization are ap¬ 
proximated by the Born momenta qi). The dependence 
on the underlying hard interaction is encoded in the hard 
function #«({&}), where {q,} = {qi, ..., q N +L}, the sum 
over k = {n a ,Kb, ■ ■ ■ kn+l} is over all relevant partonic 
processes, and the trace is over color. Any dependence 
probing softer momenta, such as measuring jet masses 
or low ptS, as well as the choice of jet algorithm, will 
affect the precise form of the factorization, but not the 
hard function H K . This dependence enters through the 
definition of the soft function S K (describing soft radia¬ 
tion), jet functions J KJ (describing energetic final-state 
radiation in the jets) and the beam functions Bi (de¬ 
scribing energetic initial-state radiation along the beam 
direction). More precisely, the beam function is given by 
Bi = (JTL, P-u' ® fi' with fi the parton distributions of 
the incoming protons, and In’ a perturbatively calcula¬ 
ble matching coefficient depending on the measurement 
definition j94[ . The ellipses at the end of Eq. © denote 
power-suppressed corrections. All functions in the factor¬ 
ized cross section depend only on the physics at a single 


1 Here active parton refers to initial-state quarks or gluons. Proofs 

of factorization with initial-state hadrons must also account for 
effects due to Glaubers ftTij . which may or may not cancel, and 
whose relevance depends on the observable in question 192 , 93;1 . 


scale. This allows one to evaluate all functions at their 
own natural scale, and then evolve them to a common 
scale using their RGE. This procedure resums the large 
logarithms of scale ratios appearing in the cross section 
to all orders in perturbation theory. 

The explicit form of the factorization theorem in 
Eq. (0, including field-theoretic definitions for the jet, 
beam, and soft functions is known for a number of exclu¬ 
sive jet cross sections and measurements of interest. For 
example, factorization theorems exist for t he IV -jet cross 
section defined using iV-jettiness Izzl Mm- These 
have also been utilized to include higher-order resumma¬ 
tion in Monte Carlo programs |67H69j . and are the ba¬ 


sis of the ./V-i ettiness s ubtraction method for fixed-order 
calculations |l02l . 03. In addition, there h as b een a 
focus on color-singlet production at small q? jl04l - tl0^ |. 
as well as the f actorizati on of processes defined with jet 
algorithms 0, H 0 M 25 , jet shape v ariables 
or fragmentation properties |l37l - ll45 | for identified jets. 


The same hard functions also appear in threshold resum¬ 
mation factorization formulas, which are often used to 
obtain an approximate higher order result for inclusive 
cross sections. ^ 

The focus of our paper is the hard function H K ({qi}) in 
Eq. 0 , which contains the process-dependent underlying 
hard interaction of Eq. 0, but is independent of the 
particular measurement. In SCET, the dependence on 
the hard interaction is encoded in the Wilson coefficients, 
C, of a basis of operators built out of SCET fields. The 
Wilson coefficients can be calculated through a matching 
calculation from QCD onto the effective theory. The hard 
function appearing in the factorization theorem is then 
given by 

£«({*}) = E^ 1 -(-A„)(te})^ 1 .. ( .. A „ ) (te})- (3) 

{AH 


Here, the {Ai} denote helicity labels and the sum runs 
over all relevant helicity configurations. The C are vec¬ 
tors in color space, and the hard function is therefore a 
matrix in color space. 

For processes of higher multiplicities, the construction 
of a complete basis of SCET operators, and the subse¬ 
quent matching calculation, becomes laborious due to 
the proliferation of Lorentz and color structures, simi¬ 
lar to the case of high-multiplicity fixed-order calcula¬ 
tions using standard Feynman diagrams. The use of 
SCET helicity fields introduced in this paper, combined 
with analogous color management techniques as used in 
the calculation of amplitudes, makes the construction of 
an operator basis extremely simple, even in the case of 
high-multiplicity processes. Furthermore, with this basis 
choice, the SCET Wilson coefficients are precisely given 
by the Ill-finite parts of the color-ordered QCD helic¬ 
ity amplitudes, rendering the matching procedure almost 
trivial. Combining the results for the hard function with 
known results for the soft, jet, and beam functions, then 
allows for the resummation of jet observables in higher 
multiplicity processes, which are ubiquitous at the LHC. 
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The remainder of this paper is organized as follows. In 
Sec. Ill Al we review the notation for the spinor-helicity 
formalism. Additional useful lrelicity and color identi¬ 
ties can be found in Appendix [A] We provide a brief 
summary of SCET in Sec. IIIBI In Sec. Mil we introduce 
SCET helicity fields and operators, and describe the con¬ 
struction of the helicity and color basis, as well as its 
symmetry properties. In Sec. IIVI we discuss the match¬ 
ing from QCD onto the SCET helicity operators, includ¬ 
ing a discussion of the dependence on the regularization 
and renormalization scheme. We then demonstrate the 
matching explicitly for H+ 0,1, 2 jets in Sec. El V+0,1, 2 
jets in Sec. IVII and pp —> 2,3 jets in Sec. IVIII Explicit 
results for the required helicity amplitudes are collected 
in the appendices. In Sec. IVIIII we discuss the general 
renormalization group evolution of the hard coefficients, 
which involves mixing between different color structures, 
to all orders. We give explicit results for the anomalous 
dimensions for up to 4 colored particles plus an arbitrary 
number of uncolored particles. We conclude in Sec. m 


II. NOTATION 
A. Helicity Formalism 

We will use the standard notation for the spinor alge¬ 
bra (for a review see for example Refs. |l46l . Il47j b Con¬ 
sider the four-component spinor u(p) of a massless Dirac 
particle with momentum p, satisfying the massless Dirac 
equation, 

i>u{jp) = 0 , p 2 * = 0. (4) 

The charge conjugate (antiparticle) spinor v(p) also sat¬ 
isfies Eq. ©, and we can choose a representation such 
that v(p) = u(p). The spinors and conjugate spinors for 
the two helicity states are denoted by 

\p±)= l -±J^u{p), 

(P ±I = sgn(p°) u{p) 1 ^ 75 . (5) 

For massless particles chirality and helicity agree while 
for antiparticles they are opposite, so |p+) = u+(p) = 
V- (p) corresponds to positive (negative) helicity for par¬ 
ticles (antiparticles). The spinors |p±) are defined by 
Eqs. © and © for both physical (p° > 0) and unphysical 
(p° < 0) momenta. Their explicit expression, including 
our overall phase convention, is given in Appendix IA II 
The spinor products are denoted by 

(pq) = {p-\q +), \p<l\ = <P+k“) • ( 6 ) 

They satisfy 

{pq) = ~{qp), [pq] = ~[qp }, (pq)[qp\ = 2p-q- ( 7 ) 

Additional relations are collected in Appendix I A II The 
minus sign for p° < 0 in Eq. © is included so the spinor 


relations are invariant under inverting the signs of mo¬ 
menta, p —» —p, when crossing particles between the ini¬ 
tial and final state, e.g. ((— p)q)[q{— p)] = 2(—p) • q. 

If there are several momenta pi, it is common to ab¬ 
breviate 


IPi±) = l*±) > {PiPj) = { ij ), [ PiPj} = [ij] • (8) 


The polarization vectors of an outgoing gluon with mo¬ 
mentum p are given in the helicity formalism by 


£+{p,k) 


(p+|7 M |fc+) 

y/2{kp) 


£-{p,k) = 


(p-\Y\k~) 

V2[kp\ 


(9) 

where k is an arbitrary reference vector with k 2 = 0, 
which fixes the gauge of the external gluons. Using the 
relations in Appendix IA II it is easy to check that 


P ' £±(j>, k) = k ■ e±(p, k) = 0 , 

£±(p, k) ■ £±(p, k) = 0, 

£±{p,k) -e T (p,k) = -1, 

£±{p,k) = e T (jp,k ), ( 10 ) 

as is required for physical polarization vectors. With 
p^ = E(l, 0,0,1), the choice k M = £7(1,0,0,— 1) yields 
the conventional 

4(p,fc) = -^( 0,l,Ti,0). (11) 


B. SCET 

Soft-collinear effective theory is an effective field the¬ 
ory of QCD that describes the interactions of collinear 
and soft particles [5^ - l60| in the presence of a hard inter¬ 
actional Collinear particles are characterized by having 
large energy and small invariant mass. To separate the 
large and small momentum components, it is convenient 
to use light-cone coordinates. We define two light-cone 
vectors 


n^ = {l,fi), n^ = (l,-n), (12) 

with n a unit three-vector, which satisfy n 2 = n 2 = 0 and 
n ■ n = 2. Any four-momentum p can be decomposed as 

fl A 4 fi A 4 * * 

P 11 = n-p — +n-p — +p^ ± . (13) 

An “n-collinear” particle has momentum p close to the 
n direction, so that p scales as (n-p,n-p, _l) ~ n-p 


2 Throughout this paper, we will fo r simplicity use the notation of 

SCETi. The theory SCETn [l48|| is required for a certain class of 

observables, for example px ~dependent measurements or vetoes. 

The helicity operator formalism presented here applies identically 

to constructing SCETn operators. The collinear operators and 

matching coefficients are the same for both cases. 
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(A 2 ,1, A), with Ada small parameter. For example, 
for a jet of collinear particles in the H direction with total 
momentum pj, n-pj — 2 Ej corresponds to the large 
energy of the jet, while n-pj ~ rrij/Ej -C Ej, where mj 
is the jet mass, so A 2 ~ m'j/Ej <C 1. 

To construct the fields of the effective theory, the mo¬ 
mentum of ?t-collinear particles is written as 

77 M 

p M =p M + = n-p— +Pn± + , (14) 

where n ■ p ~ Q and p n ± ~ A Q are the large momentum 
components, while k ~ X 2 Q is a small residual momen¬ 
tum. Here, Q is the scale of the hard interaction, and 
the effective theory expansion is in powers of A. 

The SCET fields for n-collinear quarks and gluons, 
€n,p(x) and An : p(x), are labeled by the collinear direction 
n and their large momentum p. They are written in po¬ 
sition space with respect to the residual momentum and 
in momentum space with respect to the large momen¬ 
tum components. Derivatives acting on the fields pick 
out the residual momentum dependence, i<9 M ~ k ~ A 2 Q. 
The large label momentum is obtained from the label 
momentum operator V%, e.g. V%£ n ,p = P^ £ n ,p- If there 
are several fields, V n returns the sum of the label mo¬ 
menta of all n-collinear fields. For convenience, we define 
V n = n-Vm which picks out the large momentum compo¬ 
nent. Frequently, we will only keep the label n denoting 
the collinear direction, while the momentum labels are 
summed over (subject to momentum conservation) and 
are suppressed in our notation. 

Collinear operators are constructed out of products of 
fields and Wilson lines that are invariant under collinear 
gauge transformations [H, [HI- The smallest building- 
blocks are collinearly gauge-invariant quark and gluon 
fields, defined as 


Xn,u{x) = 6(W - Vn) Wl(x) £ n (x) 

= -U(cJ + V n )W^x)iD^.W n (x)] . (15) 


9 1 


With this definition of Xn,u>, we have w > 0 for an in¬ 
coming quark and w < 0 for an outgoing antiquark. For 
B n ,u±, w > 0 (w < 0) corresponds to an outgoing (in¬ 
coming) gluon. In Eq. (fl5l) 


x, we can treat Xn,u>( x) and x) like local quark and 
gluon fields. For later use, we give the expansion of the 
collinear gluon field 


rm — av _ 

°n,± ~ 


u 

n • p 


A-n,p “t“ 


(18) 


Here the ellipses denote terms in the expansion with more 
than 2 collinear gluon fields, which are not required for 
our matching calculations. 

In our case the effective theory contains several 
collinear sectors, ni,n 2 ,... (l49| | . where the collinear 
fields for a given sector n- 1 = (1 ,ni) describe a jet in 
the direction Hi, and we also define = (1,— Hi). 
A fixed-order QCD amplitude with N colored legs is 
then matched onto operators in SCET with N different 
collinear fields. The different collinear directions have to 
be well separated, which means 


rii-rij A 2 for i ^ j. (19) 


The infrared singularities associated with collinear or soft 
limits of legs in QCD are entirely described by the La- 
grangian and dynamics of SCET itself, so the QCD am¬ 
plitudes are only used to describe the hard kinematics 
away from infrared singular limits. 

Two different rq and n' with ni-n[ ~ A 2 both describe 
the same jet and corresponding collinear physics. Thus, 
each collinear sector can be labeled by any member of 
a set of equivalent vectors, {rij | , which are related by 
reparametrization invariance [1501 ]. The simplest way to 
perform the matching is to choose rii such that the large 
label momentum is given by 



( 20 ) 


with = 0. 

In general, operators will have sums over distinct 
equivalence classes, {rii}, and matrix elements select a 
representative vector to describe particles in a particular 
collinear direction. For many leading power applications 
there is only a single collinear field in each sector, and we 
may simply set the large label momentum of that build¬ 
ing block field to that of the external parton using the 
following simple relation, 


iD n± = K± + 9 A n± - 

is the collinear covariant derivative and 


W n (x) = 


^2 ex p(-=- n ■ A n (x 

-^ r n 


L perms 


(16) 


(17) 


is a Wilson line of n-collinear gluons in label momentum 
space. The label operators V n in Eqs. m and (ED only 
act inside the square brackets. W n (x) sums up arbitrary 
emissions of n-collinear gluons from an n-collinear quark 
or gluon, which are 0(1) in the power counting. Since 
W n (x) is localized with respect to the residual position 


Jdp6(p-p) f(p) = f(ni -py) , (21) 

where p is collinear with the Fth jet. Here the tildes on 
the integration measure and delta function ensure that 
the integration over equivalence classes is properly im- 
plementedEI Because of this, at leading power, the issue 
of equivalence classes can largely be ignored. 


3 The precise definition of this delta function and measure are 
<5(pi -p) = S{ ni y, p <S(wi - f ~n-p), 
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Particles that exchange large momentum of 0(Q) be¬ 
tween different jets are off shell by 0{rii ■ n ; Q 2 ). They 
are integrated out by matching QCD onto SCET. Before 
and after the hard interaction the jets described by the 
different collinear sectors evolve independently from each 
other, with only soft radiation between the jets. The cor¬ 
responding soft degrees of freedom are described in the 
effective theory by soft quark and gluon fields, q s (x) and 
A s (x ), which only have residual soft momentum depen¬ 
dence idv ~ X 2 Q. They couple to the collinear sectors 
via the soft covariant derivative 

iD^W+gAS, (24) 

acting on the collinear fields. At leading power in A, n- 
collinear particles only couple to the n-A s component of 
soft gluons, so the leading-power n-collinear Lagrangian 
only depends on n-D s . For example, for n-collinear 
quarks 0, 59] 

Cn = 1n(in-D s + gn-A n + iJp nl _W n =- • 

(25) 

The leading-power n-collinear Lagrangian for gluons is 
given in Ref. [60j . 


III. SCET OPERATOR BASIS 

In this section, we describe in detail how to construct 
a basis of helicity and color operators in SCET, which 
greatly simplifies the construction of a complete opera¬ 
tor basis and also facilitates the matching process. Usu¬ 
ally, a basis of SCET operators obeying the symmetries 
of the problem is constructed from the fields Xn,wi 
as well as Lorentz and color structures. This process be¬ 
comes quite laborious due to the large number of struc¬ 
tures which appear for higher multiplicity processes, and 
the reduction to a minimal basis of operators quickly be¬ 
comes nontrivial. Instead, we work with a basis of op¬ 
erators with definite helicity structure constructed from 
scalar SCET building blocks, which, as we will show, 
has several advantages. First, this simplifies the con¬ 
struction of the operator basis, because each independent 
helicity configuration gives rise to an independent helic¬ 
ity operator. In this way, we automatically obtain the 
minimal number of independent operators as far as their 



( 22 ) 


where 


^{ n i 


1 n,i-p = 0{ A 2 ), 

0 otherwise. 


(23) 


The Kronecker delta is nonzero if the collinear momentum p is in 
the {rii} equivalence class, i.e. p is close enough to be considered 
as collinear with the zth jet. The sum in the second line of 
Eq. J22J runs over the different equivalence classes. 


Lorentz structure is concerned. Second, operators with 
distinct helicity structures do not mix under renormal¬ 
ization group evolution, as will be discussed in detail in 
Sec. IVIIII The reason is that distinct jets can only ex¬ 
change soft gluons in SCET, which at leading order in 
the power counting means they can transfer color but 
not spin [see Eq. Therefore, the only nontrivial as¬ 

pect of the operator basis is the color degrees of freedom. 
The different color structures mix under renormalization 
group evolution, but their mixing only depends on the 
color representations and not on the specific helicity con¬ 
figuration. 


A. Helicity Fields 


We start by defining quark and gluon fields of definite 
helicity, out of which we can build operators with a defi¬ 
nite helicity structure. To simplify our discussion we will 
take all momenta and polarization vectors as outgoing, 
and label all fields and operators by their outgoing he¬ 
licity and momenta. Crossing symmetry, and crossing 
relations are discussed in Sec. HUB 

We define a gluon field of definite lrelicitjH 


B 


a 

z± 




(26) 


where a is an adjoint color index. For nf = (1,0,0,1), 
we have 


e±(rii, m) 


x/2 


(0,1, =pi, 0), 


(27) 


in which case 




(28) 


For an external gluon with outgoing polarization vec¬ 
tor e(p, k ) and outgoing momentum p in the ?Zj-collinear 
direction, the contraction with the field Bf ± contributes 


— e- 


TM 


(ni,rii) £^_.{p,k) 


p L . _ 

- -— Hi ■ £(p, k ) , 

rii-p i 


(29) 


where we have used the expansion of the collinear gluon 
field given in Eq. m- Since (jii,rii) is perpendicular 
to both Hi and fq, we can drop the _Lj labels in brackets. 
A convenient choice for the reference vector is to take 
k = fit, for which the second term in brackets vanishes. 
Equation (l29l) then becomes 


-e T (ni,ni) • e(p,n»), (30) 


which is equal to 0 or 1 depending on the helicity of 
e{jp, fii). Adopting this choice, the tree-level Feynman 


4 The label ± on B± refers to helicity and should not be confused 
with light-cone components. 











rules for an outgoing gluon with polarization ± (so e = 
e±), momentum p (with p° >0), and color a are 


Jff = - 


\^2 8 I {lli , Tlj ) Xi— 7/v.Xj- 


(37) 


(ffi (p) | ^i± 10> = S ab 5(pi — p), 
(g“(p)\B b i± \ 0 ) = 0 . 


(31) 


J -\J idj [TiiTij] 

where = fii -fti from Eq. as well as a scalar current 


Note that B b i± = B b ± ( 0), so we do not get a phase from 
the residual momentum. Similarly, for an incoming gluon 
with incoming polarization (e = e=p, so e* = £±), in¬ 
coming momentum — p (with p° < 0), and color a, we 
have 

(0\B b ± \g a T (-p))=6 ab ~S(p l -p), 

(°\B b i±\g±{-p)) = 0. (32) 

We define quark fields with definite helicitjH as 


XT+Xj- 


J&P _ 

' ,0 s/UiUj [riiTlj] 


= 


Xi 


-x5+ 


lj0 {riirij) 


■y“ = 1 ^ 75 y" 

Ai± r> Arij.-Wi ? 


= v“ 1 Z 75 
Az± Ani,Wj 


(33) 


2 i /vizr. ^ 

where a and d are fundamental and antifundamental 
color indices respectively. 

For external quarks with n,-collinear momentum p , the 
fields contribute factors of the form 

—2 -y U (W = — b±) = b±)n 4 , (34) 

where in the last equality, we have defined a shorthand 
notation |p±) ni for the SCET projected spinor. The 
spinor | p±) ni is proportional to |n±); see Eq. (1A30I) . 

The tree-level Feynman rules for incoming (p° < 0) and 
outgoing (p° > 0) quarks with helicity +/— and color a 
are then given by 


(38) 


In Eqs. (E3 and (l38l) the flavor labels of the quarks have 
not been made explicit, but in general the two quark 
fields in a current can have different flavors (for example 
in W production). Since we are using a basis of physical 
polarization states it is not necessary to introduce more 
complicated Dirac structures. For example, pseudovector 
and pseudoscalar currents, which are usually introduced 
using 7 ' 5 , are incorporated through the relative coeffi¬ 
cients of operators involving J + , J_ or Jo, Jo- As we 
shall see, this greatly simplifies the construction of the 
operator basis in the effective theory. 

At leading power, there is a single collinear field in each 
collincar sector, so we can choose n^ = p'lIp\ to represent 
the equivalence class {rip, so that pP = -n ■ Pi nf which 
gives 


b±)n 4 = |p±) = 


•Py±> = 


m -p 


m±). (39) 


Since we always work at leading power in this paper, we 


(o \xi+\q+(-p)) = s^hPi-p)\(-Pi)+)ru , 
(o \xi-\q-(-p)) = 5 0oi 5{pi-p) 1 {-Pi)-)m , 


will always make this choice to simplify the nr; 
With this choice, the tree-level Feynman rules 
fermion currents are 

(q+(p)\Xi+ 0) = 6 a P 5(pi - p) ni (pi+ 1, 

(q~ (p) Xi- 0) = 5 aP 5(Pi - p) ni (pi- 1, 

(35) 

(q?(pi)e 2 (P2)\J^ |0> 

= <S Q1 * ~S(Pi-Pi)hP2-P2), 

and similarly for antiquarks 


(q- 1 (pi) q+*( P 2 ) \ J\ 2 - 2 0) 

(0 xf+|9-(-p)> = S a ^S(pi -p) ni {{~Pi)+\ , 


= S a ^ S (hs a j(pr-pi)5(fb-p2), 

<0 xf_ \q+(-p)) = S a P 6(pi - p) ni {{~Pi)-\ , 


(q+ 1 (pi)q‘+(p 2 )\J^o 2 \°) 

{qZ(p)\xi+ o) = $ Pa 5(pi — p) \pi+) ni , 


= S a ^ 5^ 5(pi-pi)5(p2-P2), 

(<7?(p) xf- 0) = $ 0a S(pi-p ) | Pi-)ru ■ 

(36) 

{q a P(pi)fj( P 2 )\(H^ |o) 

The corresponding Feynman rules with the helicity of the 
external (anti)quark flipped vanish. 

= <5“ l/?1 <5 /32 “ 2 <5 (pi - Pi)S(p 2 ~ P 2 ) ■ 


To avoid the explicit spinors in Eqs. (1551) and 
and exploit the fact that fermions come in pairs, we also 
define fermionic vector currents of definite helicity 

VSetfanj) Xi+ 7pxf + 


T^P _ 

J ij+ ~ 


\J OJ'i UJj ( ) 


(40) 


The simplicity of these Feynman rules arises due to 
the unconventional normalization of the operators in 
Eqs. and (l38l) . This normalization has been cho¬ 
sen to simplify the matching of QCD amplitudes onto 
SCET operators, as will be seen in Sec. hyi 
W e will also make use of leptonic versions of the above 
currents. These are defined analogously, 


5 Technically speaking chirality, although we work in a limit where 
all external quarks can be treated as massless. 


Jij+ — 


n j) 1/ •- 7/tQ • 

\J (-dj 
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J _ _ y/^E+jTli, Tlj) £j± 'y f j£j± (41) 

sj [jljTlj] 

Unlike the collinear quark field %, the leptonic field i 
does not carry color and so does not contain a strong- 
interaction Wilson line. 

All couplings in the SM, except to the Higgs boson, 
preserve chirality. This limits the need for the scalar 
current, especially when considering only massless ex¬ 
ternal quarks. In the SM the scalar current can arise 
through explicit couplings to the Higgs, in which case, 
even though we still treat the external quarks as massless, 
the Wilson coefficient for the scalar operator will contain 
the quark Yukawa coupling. This is relevant for exam¬ 
ple for Hbb processes. The scalar current can also arise 
through off-diagonal CKM-matrix elements connecting 
two massless external quarks through a massive quark 
appearing in a loop. This can occur in multiple vector 
boson production, or from electroweak loop corrections, 
neither of which will be discussed in this paper. When 
constructing an operator basis in Sec. IIIIB1 we ignore 
the scalar current, as it is not relevant for the examples 
that we will treat in this paper. However, it should be 
clear that the construction of the basis in Sec. ME can 
be trivially generalized to incorporate the scalar current 
if needed. 


B. Helicity Operator Basis 


Using the definitions for the gluon and quark helicity 
fields in Eqs. (® and (ED, we can construct operators 
for a given number of external partons with definite helic- 
ities and color. (As discussed at the end of the previous 
section, for the processes we consider in this paper we do 
not require the scalar current J 5 .) In the general case 
with CKM-matrix elements, we must allow for the two 
quark flavors within a single current to be different. The 
situation is simplified in QCD processes, where one can 
restrict to currents carrying a single flavor label. 

For an external state with n particles of definite helic- 
ities ±, colors a*, oti, oij, and flavors /, /', ..., a complete 
basis of operators is given by 


O 


a\d2--cx.i 


±±-..(± : ...±) (Pi 1PI1 • ■ ■ 1 Pi— 1 ) Pi- ■ ■ ■ • Pn—Xi Pn) 


_ c f2 a 1 K2 a 2 T&i—ioti ja n — iOt n 

— o D 1± a 2± ■ ■ ■ J■ ■ ■ -Jff n -l t n± ■ 


(42) 


For example, / = q indicates that both quark fields in the 
current have flavor q. When it is necessary to distinguish 
different flavors with the same current, for example when 
we consider processes involving W bosons in Sec. ED we 
use a label / = ud such as Judi 2 -- For simplicity, we will 
also often suppress the dependence of the operator on the 
label momenta 2 ^. For the operator subscripts, we always 
put the helicity labels of the gluons first and those of the 
quark currents in brackets, with the labels for quark cur¬ 
rents with different flavor labels / and /' separated by a 


semicolon, as in Eq. (l42l) . The ± helicity labels of the in¬ 
dividual gluon fields and quark currents can all vary inde¬ 
pendently. Operators with nonzero matching coefficients 
are restricted to the color-conserving subspace. We will 
discuss the construction of the color basis in Sec. mm 
The symmetry factor S in Eq. ED is included to sim¬ 
plify the matching. It is given by 


S = 


I! nprii ! ’ 


(43) 


where nf denotes the number of fields of type i = 
g, u, u,d,d,... with helicity ±. We also use 


n = ^2( n t+ n i )> ( 44 ) 

i 


to denote the total number of fields in the operator. Each 
Bi counts as one field, and each J has two fields. 

For each set of external particles of definite helicities, 
colors, and flavors, there is only one independent oper¬ 
ator, since the physical external states have been com¬ 
pletely specified. All Feynman diagrams contributing to 
this specific external state will be included in the Wilson 
coefficient of that specific operator. For the case of pure 
QCD, quarks always appear in pairs of the same flavor 
and same chirality, and therefore can be assembled into 
quark currents labeled by a single flavor. In this case, to 
keep track of the minimal number of independent oper¬ 
ators, we can simply order the helicity labels, and only 
consider operators of the form 


0 +■■(.—) 



and analogously for any additional quark currents with 
different quark flavors |j 

With the operator basis constructed, for a given n- 
parton process we can match hard scattering processes in 
QCD onto the leading-power hard-scattering Lagrangian 


C 


hard 


« n 

j\[dvi C£;^)(pi,... ,p n ) 0 +(pi, • 


> Pn)i 
(46) 


In the general case with off-diagonal CKM-matrix elements, 
there is some more freedom in the choice of the operator ba¬ 
sis, because quarks of the same flavor do not necessarily ap¬ 
pear in pairs. However, it is still true that only a single oper¬ 
ator is needed for a specific external state. For example, for 
external quarks u_, d+, s+, c ~, one could either use the op¬ 
erators Jus—Jed— ■> °r the operators J C s—Jud— (where the color 
structures have been suppressed). Since different helicity combi¬ 
nations are possible, a single flavor assignment does not suffice to 
construct a complete helicity basis, and one must sum over a ba¬ 
sis of flavor assignments. As an example explicitly demonstrating 
this, we will consider the case of pp —> W-\- jets in Sec. eh 
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where a sum over all color indices is implicit. Lorentz in¬ 
variance implies that the Wilson coefficients only depend 
on Lorentz invariant combinations of the momenta. This 
hard Lagrangian is used in conjunction with the collinear 
and soft Lagrangians that describe the dynamics of the 
soft and collinear modes; see for example Eq. (l25l) . 

We emphasize that Eq. (USD provides a complete ba¬ 
sis in SCET for well-separated jets and additional non- 
hadronic particles at leading power. We will discuss in 
more detail in Sec. m the matching and regularization 
schemes, and demonstrate that no evanescent operators 
are generated for this case. At subleading power, the 
SCET operators would involve additional derivative op¬ 
erators, soft fields, or multiple SCET building blocks 
from the same collinear sector. 


C. Example with a Z-Boson Exchange 


Here S12 = (pi +P2) 2 , Pz is the ratio of the Z and photon 
propagators, 


Pz{s) 


s 

s - m| + i T z m z ’ 


and the couplings Vl,r to the Z boson are 


1 — 2 sin 2 9w e 2 sin 2 6w 
sin(2 Qw) ’ Vr sin(2 6w) 


(49) 


(50) 


Note that the presence of the spinor factors in Eq. (l48l) 
occur due to our normalization conventions for the cur¬ 
rents. 

Now, consider calculating the scattering amplitude in 
the effective theory, for example for the case when both 
electrons have positive helicity. The matrix element in 
the effective theory gives 


It is important to note that all kinematic dependence of 
the hard process, for example, its angular distributions, is 
encoded in the Wilson coefficients. Since the Wilson co¬ 
efficients can (in principle) carry an arbitrary kinematic 
dependence, our choice of helicity basis imposes no re¬ 
striction on the possible structure or mediating particles 
of the hard interaction. For example, the spin of an in¬ 
termediate particle may modify the angular distribution 
of the decay products, and hence the Wilson coefficients, 
but this can always be described by the same basis of 
helicity operators. 

As a simple example to demonstrate this point we con¬ 
sider e + e _ —> e + e~ at tree level. This process can pro¬ 
ceed through either an off-shell 7 or Z boson. Because 
the SM couplings to both of these particles preserve chi¬ 
rality, a basis of operators for this process is given by 

0 (++) = J Je 12+ j e 34+ , 

0(+~) = Je 12+ Je 34- , 

0( -) = J -4 12- -4 34- j (47) 

where the leptonic current is defined in Eq. (1411) . The fact 
that this is a complete basis relies only on the fact that 
the couplings preserve chirality, and is independent of 
e.g. the possible number of polarizations of the mediating 
Z or 7 . 

We now consider the calculation of the Wilson coeffi¬ 
cients for the matching to these operators (the matching 
procedure is discussed in detail in Sec. m- At tree level, 
the Wilson coefficients are easily calculated, giving 

C(++) = ^ e2 [l + V e R V e R Pz(si 2 )\ ————- + (1 -O' 3), 

s 12 

C (+ _) = -e 2 [l + «P z ( Sl2 )]ffiflM, 

s 12 

C(__) = -e 2 [1 + vlvtP z (s 12 )] 2[24](13) + (1 o 3). 

S 12 

( 48 ) 


(e + (pi) et (p 2 ) e + ( p 3 ) et {pa) | i£hard 10) 

r. n 

;(pi)et(p 2 )e+(p 3 )et(p 4 )| j \\_JpiC ++ 0 


= 1 e 


2=1 


= -ie 2 [1 + v e R v e R P z ( Sl2 )\ 2[13](24) + (1 o 3), (51) 


s 12 


using the Feynman rules of Eq. (l40l) . The effective theory 
therefore reproduces the full theory scattering amplitude. 
The same is true of the other helicity configurations, so 
the familiar angular distributions for e + e _ —> e + e _ , as 
well as the different couplings of the Z to left- and right- 
handed particles, are entirely encoded in the Wilson co¬ 
efficients. 


D. Color Basis 


In addition to working with a basis of operators with 
definite helicity, we can also choose a color basis that fa¬ 
cilitates the matching. When constructing a basis of op¬ 
erators in SCET, we are free to choose an arbitrary color 
basis. With respect to color, we can think of Eq. TO as 
having a separate Wilson coefficient for each color con¬ 
figuration. For specific processes the color structure of 
the Wilson coefficients can be further decomposed as 

°+-(»-) = zL G +-(—)W- = 1 c +"(-) ■ 

k 

(52) 

Here, T ai "'“ n is a row vector whose entries are 

suitable color structures that together provide a complete 
basis for all allowed color structures, but which do not 
necessarily all have to be independent. In other words, 
the elements of T ai span the color-conserving sub¬ 
space of the full color space spanned by {ai • • • a n }, and 
C is a vector in this subspace. Throughout this paper 
we will refer to the elements of T ai as a color ba¬ 
sis, although they will generically be overconrplete, since 
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this allows for simpler choices of color structures. As dis¬ 
cussed below, due to the overcompleteness of the bases, 
some care will be required for their consistent usage. 
Using Eq. (El, we can rewrite Eq. (® as 


Uhard = J JJdpj 

._)(Pl) • ■ ■ ,Pn)C + -(—){P\, ■ 

■ ■ ,Pn), 



(53) 

where O t is a conjugate vector defined by 


°T<-) 

/Qfll ■■■CXn 'T'ai ■ ■ -OLn 
- °+..(.-) 1 

(54) 


While the form C® 1 (“jy 0°£ ]in Eq. (l46l) is more con¬ 
venient to discuss the matching and the symmetry prop¬ 
erties of operators and Wilson coefficients, the alternative 
form in Eq. (1531) is more convenient to discuss the mixing 
of the color structures under renormalization. 

For low multiplicities of colored particles it can be con¬ 
venient to use orthogonal color bases, e.g., the singlet- 
octet basis for qqq'q' is orthogonal. However, using or¬ 
thogonal bases becomes increasingly difficult for higher 
multiplicity processes, and the color bases used for many 
fixed- order calc ulations are not orthogonal. (See e.g. 
Refs. EMEU for a discussion of the use of orthogo¬ 
nal bases for SU(N).) The use of a nonorthogonal color 
basis implies that when written in component form in a 
particular basis, the conjugate (A of the vector C is not 
just given by the naive complex conjugate transpose of 
the components of the vector. Instead, we have 


_ ^jCLi'-Oin j* rj-}a\-"CX. ri (J *T rj-\ 


(55) 


which follows from its color basis independent definition 
in terms of Wilson lines [see e.g. Ref. j98( or Eq. (12671) ]. 
Here we have suppressed the dependence of S on soft 
momenta. The action of the identity on an element of 
the color space is defined by 

i 17 ) , (60) 

and its matrix representation in any color basis is given 
by 1 = diag(l, 1, ■ ■ • , 1) . In the literature, see e.g. 
Refs. [82L 88l l89l 1 1 34l fl53[ | . often a different convention is 
used, where the T matrix is absorbed into the definition 
of the soft function. In this convention, the soft func¬ 
tion becomes explicitly basis dependent and is not the 
same as the basis-independent color-space identity. One 
should be careful to not identify the two. 

As an example to demonstrate our notation for the 
color basis, consider the process ggqq. A convenient 
choice for a complete basis of color structures is 

faba p = (jr a T b ) a p , (T b T a ) a p , tr[T°T & ] S a? ) 

( ( T a T b ) a 0 \ T 

= ( T b T a ) aB ■ (61) 

ytr[T a T b ] 5 a pJ 

For cases with many color structures we will write T as 
the transpose of a column vector as above. The transpose 
in this case only refers to the vector itself, not to the 
individual color structures. The color-sum matrix for 
this particular basis is 


where 

T= (T' ai "'“") t T ai (56) 

ai,...,Q! n 

is the matrix of color sums for the chosen basis. If the 
basis is orthogonal (orthonormal), then T is a diagonal 
matrix (identity matrix). Note that Eq. (15^1) implies that 
by definition T* t = T. 

Similar to Eq. E3b for an abstract matrix X in color 
space, the components of its Hermitian conjugate X' 
when written in a particular basis are given in terms of 
the components of X as 

Xt = f-ix* T f. (57) 


A proper treatment of the nonorthogonality of the 
color basis is also important in the factorization theo¬ 
rem of Eq. ©• Here, the color indices of the Wilson 
coefficients are contracted with the soft function as 


•••OL n j * CjCLi--(Xnbi •••/3 n 


c^s K c 

C* t TS k C. 


(58) 


At tree level, the soft function is simply the color-space 
identity 

(59) 


_ _ ^rj-taba/3^ rj-iaba^ 


c f n 

2 


/ 2C f 2 C F - C A 27> \ 
2 C F - C A 26V 27> 

y 2 TV 2 T f 2 T f N J 


(62) 


Our conventions for color factors are given in Ap¬ 
pendix [A2] Explicit expressions for T for the bases used 
in this paper are given in Appendix [Fl for up to five par- 
tons. 

Depending on the application, different choices of color 
basis can be used. For example, in fixed-order QCD cal¬ 
culations, color ordering j5|-[l] is often used to organize 
color information and simplify the sing ularit y structure of 
amplitudes, while the color flow basis tl54( | is often used 
to interface with Monte Carlo generators. For a brief 
review of the color decomposition of QCD amplitudes, 
see Appendix IA 31 Choosing a corresponding color basis 
in SCET has the advantage that the Wilson coefficients 
are given directly by the finite parts of the color-stripped 
helieity amplitudes, as defined in Eq. ESI)- which can 
be efficiently calculated using unitarity methods. In this 
paper we will use color bases corresponding to the color 
decompositions of the QCD amplitudes when giving ex¬ 
plicit results for the matching coefficients, although we 
emphasize that an arbitrary basis can be chosen depend¬ 
ing on the application. 


S K = 1, 
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Finally, note that the color structures appearing in the 
decomposition of a QCD amplitude up to a given loop 
order may not form a complete basis. The color basis in 
SCET must be complete even if the matching coefficients 
of some color structures are zero to a given loop order, 
since all structures can in principle mix under renormal¬ 
ization group evolution, as will be discussed in Sec. IVIIII 
In this case, we always choose a complete basis in SCET 
such that the color structures appearing in the ampli¬ 
tudes to some fixed order are contained as a subset. 


can infer that the Wilson coefficients must satisfy 

X^±)(PI’P2-,P3,Pa) = - C W 2 w(Pl^ 2 ;P4,P3) • 

( 66 ) 

In the color basis of Eq. ED, we can write this as 
C\ 1 \ 2 (±)(pi,p2-,P3,Pi) = VC Xl X 2 ( T )(p 1,P2!P4,P3) , 

~ / 0 -1 o\ 

with V = I —1 0 0 I . (67) 

V° 0 -v 


E. Parity and Charge Conjugation 

Under charge conjugation, the fields transform as 

f '1 V2(l rpa f '1 Y? a rpa 

^ — ~°i± 1 / 3 a J 

(63) 

The minus sign on the right-hand side of the second equa¬ 
tion comes from anticommutation of the fermion fields. 

Under parity, the fields transform as 

PB? ± (p i ,x)P = e ±2i ^B? T (pr,x F ), 

= (64) 

where we have made the dependence on pi and x explicit, 
and the parity-transformed vectors are pf = Wjfij/2, 
The cf> ni are real phases, whose exact defi¬ 
nition is given in Appendix IA II The phases appearing 
in the parity transformation of the helicity operators ex¬ 
actly cancel the phases appearing in the corresponding 
helicity amplitude under a parity transformation. This 
overall phase is determined by the little group scaling 
(see Appendix I A 1 1 for a brief review). 

Using the transformations of the helicity fields under 
parity and charge conjugation in Eqs. ED and ED, 
it is straightforward to determine how these discrete 
symmetries act on the helicity operators. Parity and 
charge conjugation invariance of QCD implies that the 
effective Lagrangian in Eq. (® must also be invariant. 
(For amplitudes involving electroweak interactions, par¬ 
ity and charge conjugation invariance are explicitly vi¬ 
olated. This is treated by extracting parity and charge 
violating couplings from the operators and amplitudes. 
See Sec. [VI] for a discussion.) This then allows one to de¬ 
rive corresponding relations for the Wilson coefficients. 

To illustrate this with a nontrivial example we consider 
the ggqq process. The operators transform under charge 
conjugation as 

GO a x b i f 2 0 {±) (p 1: p 2 -,p 3 ,p^f ab ^ C 
= CSB a 1Xi B b 2X2 J^ ± T ab ^C 
= -° b XlX2 ( T )(Piip2;p4,P3)T aba P , (65) 

where Ay 2 denote the gluon helicities, and T aba P is as 
given in Eq. ED- From the invariance of Eq. (1461) we 


Now consider the behavior under parity. For concrete¬ 
ness we consider the case of positive helicity gluons. The 
operators transform as 

P 0“ 6 + %)(pi,P2;P3,p 4 ) P = P \B a 1+ B b 2+ J^ ± P (68) 
= e i ( 2 ^u 2 ^ ± (^3-^ 4 ))0“^ )( pf- 

The invariance of Eq. (1461) under parity then implies that 
the Wilson coefficients satisfy 


C++(±)(P 1,P2;P3,P4) 

= d__ c=F) (p f,pf;pf,pDe- i( 2 ^ + 2 ^ ±( ^3-^ 4 )) 


= C—(t)(pi>P 2 !P 3 ,P 4 ) 


(■■}«-[■■] 


(69) 


Here we have introduced the notation (..) o [..] to indi¬ 
cate that all angle and square spinors have been switched 
in the Wilson coefficient. The fact that the phase appear¬ 
ing in the parity transformation of the operator exactly 
matches the phase arising from evaluating the Wilson 
coefficient with parity related momenta is guaranteed by 
little group scaling, and will therefore occur generically. 
See Eqs. (1A24I) and (IA25[) and the surrounding discussion 
for a review. 

Below we will use charge conjugation to reduce the 
number of Wilson coefficients for which we have to carry 
out the matching explicitly. We will use parity only when 
it helps to avoid substantial repetitions in the matching. 


F. Crossing Symmetry 

Our basis is automatically crossing symmetric, since 
the gluon fields £>,;± can absorb or emit a gluon and 
the quark current Jij± can destroy or produce a quark- 
antiquark pair, or destroy and create a quark or anti¬ 
quark. We will first illustrate how to use crossing sym¬ 
metry in an example and then describe how to technically 
have crossing symmetric Wilson coefficients. 

We will again consider the process ggqq as an example. 
Due to our outgoing conventions, the default Wilson co¬ 
efficient is for the unphysical processes with all outgoing 
particles: 

0 -4 gl(pi)g b _{p2)q+{p3)q-{pi) ■ Cf^ +) (pi,p 2 ;p3,p4 ,), 

(70) 
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where we picked one specific helicity configuration for 
definiteness. Crossing a particle from the final state to 
the initial state Hips its helicity, changes the sign of its 
momentum, and changes it to its antiparticle. In addi¬ 
tion we get a minus sign for each crossed fermion, though 
in practice these can be ignored as they do not modify 
the cross section. This allows one to obtain the Wilson 
coefficient for any crossing. For example, for the follow¬ 
ing possible crossings, the Wilson coefficients are given 
by 

g+(pi)g b -(p 2 ) -t q+(p3)q-(p4 ): C b ^ +) (-p 2 ,-pi;P3,P4), 
g%(pi)q+(p2) -> g b +(p3)q+(p 4 ) : ~ c +-(+)(P3, ~p^p^ ~p^ > 

g+{pi)q-{p2) -> g+{p3)q-(p4) ■ -C b ^ +) (p 3 , -pi; -p 2 ,P4) , 

q+(pi)q-{p2) -> g+(p3)g b -(p4) ■■ Cf?“ +) (P3,P4\-P2,-Pi) • 

(71) 

Since the signs of momenta change when crossing par¬ 
ticles between the final and initial state, care is required 
in taking the proper branch cuts to maintain crossing 
symmetry for the Wilson coefficients. In terms of the 
Lorentz invariants 

Sij = (pi+Pj) 2 (72) 

this amounts to the choice of branch cut defined by s l: j —> 

Sij + iO. In particular, we write all logarithms as 

L ij = ln ( - ^l ~ i0 ) = ln (^) ~ i 7 r 0 ( s b) • ( 73 ) 


channel n is given by 

H K ({Pi}) = <V(..a b )({ Pi}) ^..(.^^({Pi}) - ( 76 ) 

{>.} 

where {pi} = {pi,p 2 , • • For unpolarized experiments 
we simply sum over all helicity operators, so i7 K ({pi}) 
with its sum over helicities in Eq. (1761) appears as a mul¬ 
tiplicative factor. It is important to note that the color 
indices of the Wilson coefficients are not contracted with 
each other, rather they are contracted with the color 
indices of the soft function through the trace seen in 
Eq. (@. 

As an explicit example to demonstrate the treatment 
of both color and helicity indices, we consider the contri¬ 
bution of the ggqq partonic channel to the pp —> 2 jets 
process. In this case, the Wilson coefficients are given by 
Cx 1 \ 2 (\ 3 ), where Ai,A 2 denote the helicities of the glu¬ 
ons, A 3 denotes the helicity of the quark current, and re¬ 
call that the vector denotes the possible color structures, 
which were given explicitly for this case in Eq. dim) . The 
hard function for this partonic channel is then given by 

HggqqdPi}) ~ y^fi\iA 2 (A 3 )({Pi}) ^AiA 2 (A 3 )^^ i J'^ 

{Ai} 

= ^++(+)^++(+) + ^+-(+)^+-(+) + 
^-+(+)^!f(+) + c —(+)C f l_(+)+ 

C ! ++ (-)Ci +( _ ) + C ! + _ ( _ ) (7|_ ( _ ) + 

c ! - + (-)Cl +( _ ) + c__ ( _ ) cL ( _ ) . ( 77 ) 


For spinors, crossing symmetry is obtained by defining 
the conjugate spinors (p±| as was done in Eq. ([5]), result¬ 
ing in the following relation 

(P±| = sgn(p°) |p±). (74) 

The additional minus sign for negative p° is included 
to use the same branch (of the square root inside the 
spinors) for both spinors and conjugate spinors, i.e., for 
p° > 0 we have 

I(-P)±) =i|p±>, 

<(-p)±I = ~(-i)(P±l = i(P±l • (75) 


Here, explicit expressions are only needed for 

C + _l_( + ), (7 + _( + ) and C _( + ). One can obtain C _(_( + ) 

using Bose symmetry simply by interchanging the 
gluons, 

C ab +f +) (pi,P2]P3,P4) = C b ™f +) {p 2 ,pi]p3,pi ), (78) 

or equivalently, 

(?_+(+)(pi,p 2 ;p 3 ,P 4 ) = VC + _ {+) (p 2 ,Pi;P 3 ,P 4 ), 

- ( 01 °\ 

with V = 1 0 0 . (79) 

V° 0 V 


In this way all spinor identities are automatically valid 
for both positive and negative momenta, which makes it 
easy to use crossing symmetry. 


G. Hard Function 

In the factorized expression for the cross section given 
in Eq. m, the dependence on the underlying hard Born 
process appears through the hard function H K . In terms 
of the Wilson coefficients of the operator basis in the ef¬ 
fective theory, the hard function for a particular partonic 


As explained in Sec. IVII A~ 2 l the remaining C\i a 2 (-) can 
be obtained from the expressions for the other Wilson 
coefficients by charge conjugation. 

In Eq. (1771) . the Wilson coefficients are vectors in the 
color basis of Eq. m and thus the hard function is a 
matrix in this basis. As discussed in Sec. IIII Dl the tree- 
level soft function is the color-space identity, i.e., 


^,( 0 ) bib 2 f 3 i (32 a\a 2 OL\OL 2 
D ggqq 


5 biai 8 b2a2 Sp l0ll Sp 2 a 2 = 1 • (80) 


With the color trace in Eq. m this amounts to contract¬ 
ing the color indices of the Wilson coefficients. In the 
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color basis of Eq. this simply becomes 


c(0) 

°ggqg 


= 1 = 



(81) 


The tree-level soft function also has dependence on mo¬ 
menta depending on the measurement being made, which 
are not shown here. 

To demonstrate a complete calculation of the cross sec¬ 
tion using the factorization theorem of Eq. © together 
with the hard functions computed using the helicity op¬ 
erator formalism, it is instructive to see how the leading- 
order cross section is reproduced from Eq. ©. We con¬ 
sider the simple case of H + 0 jets in the rrit —> oo limit. 
For this channel, there is a unique color structure S ai a 2 , 
and using the results of Sec. IV Al and Appendix IB II the 
lowest order Wilson coefficients are given by 

<?++(»,»;« = ( 8 2 ) 

<?—(pi.feips) = 

C+-(pi,p2;P3) = C-+(p 1 ,p 2 ;p3) = 0 , (84) 

where v = (>/2 Gf) ^ = 246GeV. Note that these are 
simply the helicity amplitudes for the process, as will be 
shown more generally in Sec. IIVI Analytically continu¬ 
ing to physical momenta, squaring, and summing over 
helicities, the tree-level hard function is given by 


tt( 0) oio 2 b\b 2 
H ggH 


( Pl,P2',P3 ) 


Qs S 12 

37 tv 2 
2 2 

ajsi2 c 
18tt 2 v 2 


2 

2 ^0102 
a 1 a 2 s bl b 2 


Sb 1 b 2 


(85) 


Note that only 2 of the 4 helicity configurations con¬ 
tribute, hence the factor of 2 . 

The tree-level gluon beam functions are given by the 
gluon PDFs. Since there are no jets in the final state, 
there are no jet functions. The tree-level soft function is 
the identity in color spac^H 

^( 0 ) hih 2 01O2 _ ^biQijb 2 a2 (gyj 

The leading-order cross section is then given by 


1 1 

2 E 2 m [2(TV 2 - l )] 2 


dxi dx 2 
Xi x 2 


fg{xi)fg{ X 2 ) 


7 Since there is only one color structure, the tree-level soft function 
is normally defined as 

s$ = S ^ blb2 S b ^s b ^ = 1 . (86) 

Here we do not absorb numerical prefactors into our soft func¬ 
tions, because this is not useful for processes with more final-state 
partons. 


X J~^jz e ( p ^ 6 (P3~ m H) 

x (2n) 4 5 4 (^XiE cm 1 -^-+X2E cni 1 -Y~p 3 S j 


tt(0) aia 2 bib 2 

x n ggH 

a2 s m H f 
576 nv 2 E 2 m J 


(Pl,P2lP3)Sif blb2 0102 




m H 

E C m 



( 88 ) 


The l/(2i? 2 m ) factor is the flux factor and for each of the 
incoming gluons we get a 1/[2(iV 2 — 1)] from averaging 
over its spin and color. This is followed by integrals over 
the gluon PDFs, f g , and the Higgs phase space, where 
we have restricted to the production of an on-shell Higgs. 
The final expression in Eq. (l 88 l) agrees with the standard 
result, where the first factor is the Born cross section. 

We now briefly discuss our choice of normalization. 
The currents in Eq. (1371) were normalized such that the 
Wilson coefficients are simply given by the finite part of 
the QCD helicity amplitudes (see Eq. HMD and Sec. EYD. 
This is distinct from the normalization typically used for 
SCET operators, e.g. XiTXji which is chosen to facilitate 
the matching to QCD operators. We now show that the 
extra factors in Eq. (1371) arrange themselves to produce 
the standard normalization for the jet function (or beam 
function). Starting from the current and its conjugate, 


t“/3 ( r7<5 \f 
J ij±\ J ij±) 

'J2e -p (//,;, Tl j ) Xt±7gXj± 


(89) 


\/2e±(u-i, n-j) Xj±lvXii 


\J Mi ujj (riiT \ n jE) ajj) ( n jE 


= 4 


5 ia 5& s £±(n i ,n j )et! F (n i ,n j ) r rfi if 


N N 
if. 


2ni-rij |uii uij \ 

i: 


tr 




(- Ti \ ( - TO 

x ^Xi± ~2 Xi± J ^Xj± “ 2 ” Xj± 

2610505 (2ivkET^ i=b ^ Xi± )(2AU^ ± l Xi± ) ’ 


where we have rearranged the expression in a factorized 
form using the SCET Fierz formula in spin 


1 ( 8)1 = 


ifi 

~2 2 “ 


ii 75 

2 


2 


2 J’ 

(90) 


which applies for the SCET projected spinors. In the 
last line of Eq. (l89l) we have dropped the color nonsin¬ 
glet terms and terms which vanish when averaging over 
helicities, which are indicated by ellipses. The delta func¬ 
tions in color space highlight that the jet function does 
not modify the color structure. The factor 1 /u>ij, which 
arises from the normalization of the helicity currents, is 
part of the standard definition of the jet function and en¬ 
sures that this operator has the correct mass dimension. 
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IV. MATCHING AND SCHEME DEPENDENCE 

In this section, we discuss the matching of QCD 
onto the SCET helicity operator basis introduced in the 
previous section. We start with a discussion of the 
matching for generic helicity operators in Sec. IIV Al In 
Sec. IIVBI we discuss in detail the subject of renormal¬ 
ization schemes, and the issue of converting between reg¬ 
ularization/renormalization schemes commonly used in 
spinor-helicity calculations, and those used in SCET. We 
also demonstrate that evanescent operators are not gen¬ 
erated in our basis. 


A. Generic Matching 

In this paper, we work to leading order in the power 
counting, which means we only require operators that 
contain exactly one field per collinear sector. That is, 
different n, in Eq. (H^l) are implicitly restricted to be¬ 
long to different equivalence classes, {rij} ^ {rij} for 
i ^ j. Operators with more than one field per collinear 
direction are power-suppressed compared to the respec¬ 
tive leading-order operators that have the same set of 
collinear directions and the minimal number of fields. 

At leading order, the Wilson coefficients can thus be 
determined by computing matrix elements of Eq. (H51) . 
with all external particles assigned well-separated mo¬ 
menta, so that they belong to separate collinear sectors. 
The only helicity operator that contributes in this case is 
the one that matches the set of external helicities, picking 
out the corresponding Wilson coefficient. Since we only 
have one external particle per collinear sector, we can 
simply choose n,; = Pi/p° in the matching calculation to 
represent the equivalence class {rii}. 

To compute the matrix element of £hard, we first note 
that the helicity operators are symmetric (modulo mi¬ 
nus signs from fermion anticommutation) under simulta¬ 
neously interchanging the label momenta and indices of 
identical fields, and the same is thus also true for their 
Wilson coefficients. For example, at tree level 

(g a + 1 {pi)gT(P2)\o b ^\o) tiee 

= \[s aibl 6 a ^~6(p 1 -p 1 )~6(P2-P2) 

+ 6 a ^ - P2 ) 5(P2 ~ Pi)] (91) 

so the tree-level matrix element of £hard gives 

(g+ {pi) 9+ i.P2) |£hard |0) tree (92) 

= \ [C??(P1,P2) + C? + 01 (P2,P1)] = C??(P1,P2) ■ 

By choosing m = Pi/p the label momenta pi on the 
right-hand side simply become pi = n ■ Pirii/2 = pt. 

Taking into account the symmetry factor in Eq. (1431) . 
one can easily see that this result generalizes to more 


than two gluons or quark currents with the same helic¬ 
ity. In the case of identical fermions, the various terms 
in the operator matrix element have relative minus signs 
due to fermion anticommutation which precisely match 
the (anti)symmetry properties of the Wilson coefficients. 
Hence, the tree-level matrix element of £i iar d is equal to 
the Wilson coefficient that corresponds to the configura¬ 
tion of external particles, 

( 9l92 ' ' ' Qn — lQn |£hard | 

= Cl^Z) n - lSn (Pl,P2, ■ • ■ ,Pn-l,Pn) ■ (93) 

Here and below, gt = g±(pi) stands for a gluon with 
helicity ±, momentum p,, color a,;, and analogously for 
(anti)quarks. From Eq. (1931) we obtain the generic tree- 
level matching equation 

(pi,...,Pn) = -L4 tree (Si ■•■?«), (94) 

where _4 tree denotes the tree-level QCD helicity ampli¬ 
tude. Intuitively, since all external particles are ener¬ 
getic and well separated, we are away from any soft or 
collinear limits and so all propagators in the QCD tree- 
level diagram are far off shell and can be shrunk to a 
point. Hence, the tree-level diagram simply becomes the 
Wilson coefficient in SCET. 

The above discussion can be extended to higher orders 
in perturbation theory. In pure dimensional regulariza¬ 
tion (where e is used to simultaneously regulate UV and 
IR divergences) all bare loop graphs in SCET are scale¬ 
less and vanish. Here the UV and IR divergences pre¬ 
cisely cancel each other, and the bare matrix elements 
are given by their tree-level expressions, Eq. (IM1) . In¬ 
cluding the counterterm do(euv) due to operator renor¬ 
malization removes the UV divergences and leaves the IR 
divergences. Schematically, the renormalized loop ampli¬ 
tude computed in SCET using £hard is 

-dsCET = /’((d t ) tree + (d t ) l00P )i(7 = [1 + do (hr)] id, 

(95) 

where we used that the loop contribution is a pure coun¬ 
terterm and thus proportional to the tree-level expres¬ 
sion. In general, the counterterm So is a matrix in color 
space, as we will see explicitly in Sec. IVIIII and Ap¬ 
pendix [G] By construction, the 1/e IR divergences in the 
effective theory, C do(em), have to exactly match those of 
the full theory. Therefore, beyond tree level the matching 
coefficients in MS are given by the infrared-finite part of 
the renormalized full-theory amplitude, M r en, computed 
in pure dimensional regularization. The IR-finite part is 
obtained by multiplying A le n by SCET MS renormaliza¬ 
tion factors, which cancel the full theory 1/hr poles. De¬ 
composing the renormalized QCD amplitude in a color 
basis so that -4//,)"“" = T ai '" an A r en(gi • • • q n ), the all¬ 
orders form of Eq. (1M1) becomes 

C/|_A) (ftl ) * * * ) Pn ) lA fin {gi ' ' ' Qn) 


(96) 
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= -iT ttl '" Q! "Z c 1 A-en(gl •'■;») 

^ 7^9 / ^ ^^9 / ^ 

£ -A 

The SCET renormalization factors Zc , and are 
discussed in Sec. IVIII Al At one-loop order this corre¬ 
sponds to taking (—L4“ e 1 n ''“ n ) and simply dropping the 
1/eiR terms. In Sec. IIVB1 we will discuss in more detail 
the use of different renormalization schemes to compute 

*^ren(gl * * * Qn) ■ 

If the same color decomposition is used for the QCD 
amplitude as for the Wilson coefficients in Eq. (l52l) . we 
can immediately read off the coefficients C in this color 
basis from Eq. (1961) . As an example, consider for simplic¬ 
ity the leading color n gluon amplitude, which has the 
color decomposition (see Appendix I A 3D 

An(gi---g n ) = igr 2 E tr[r°^) ■ • -T°'W] 

CreSn/Zn 

. <r( fl )), (97) 

i 

where the first sum runs over all permutations a of n 
objects (S n ) excluding cyclic permutations (Z n ). The 
An' 1 are the color-ordered or partial amplitudes at i loops. 
Each is separately gauge invariant and only depends on 
the external momenta and helicities (pj±) = (i±). If we 
choose 

Tg 1 ' an = tr ■ ■ ■ T a ”kW] ? (98) 

as the color basis in Eq. (53, where a k is the kth per¬ 
mutation in S n /Z n , then the Wilson coefficients in this 
color basis are given directly by 

c 'A 1 -A 2 (Pi ) ---.Pn) = 

9 n s~ 2 E A% n (a k ( l Al ),..., a k (n x *)), (99) 

i 

where the subscript “Hn” denotes the IR-fmite part of 
the helicity amplitude, as defined in Eq. (IMP . This is 
easily extended beyond leading color, given a valid choice 
of subleading color basis. Our basis therefore achieves 
seamless matching from QCD helicity amplitudes onto 
SCET operators. 

B. Renormalization Schemes 

In this section we discuss in more detail the issue of 
renormalization/regularization schemes in QCD and in 
SCET. In particular, the construction of a basis of helic¬ 
ity operators discussed in Sec. Illll relied heavily on mass¬ 
less quarks and gluons having two helicity states, which 
is a feature specific to 4 dimensions. We clarify this issue 
here and discuss the conversion between various schemes. 

In dimensional regularization, divergences are regular¬ 
ized by analytically continuing the particle momenta to d 


dimensions. In a general scheme, the helicities of quarks 
and gluons live in d 9 s , d q s dimensional spaces respectively. 
We shall here restrict ourselves to schemes where quarks 
have two helicities, but d 9 s is analytically continued. This 
is true of most commonly used regularization schemes, 
but is not necessary [1551 ] . Different schemes within di¬ 
mensional regularization differ in their treatment of d 9 s for 
internal (unobserved) and external (observed) particles. 
In the conventional di men sional regularization (CDR),’t 
Hooft- Veltm an (HV) |l56l | , and four-dimensional helicity 
(FDH) [l57l . Il58j ] schemes the internal/external polariza¬ 
tions are treated in d/d (CDR), d/4 (HV), 4/4 (FDH) 
dimensions. 

For helicity-based computations, the FDH scheme has 
the advantage of having all helicities defined in 4 dimen¬ 
sions, where the spinor-helicity formalism applies, as well 
as preserving supersymmetry. Indeed, most of the recent 
one-loop computations of helicity amplitudes utilize on- 
shell methods and therefore employ the FDH scheme. 
However, most existing calculations of SCET matrix ele¬ 
ments (jet, beam, and soft functions) use d-dimensional 
internal gluons, corresponding to the CDR/HV schemes^ 
As we will discuss below, CDR and HV are identical for 
matching onto SCET. 

Although the FDH scheme is convenient for helicity 
amplitude com putations, it leads to sub tleties beyond 
NLO fjfil [163. As explained in Ref. jl62| , this discrep¬ 
ancy arises due to the different number of dimensions 
for the momenta in the loop integral and the spin space, 
leading to components of the gluon Held whose couplings 
to quarks are not protected by gauge invariance and re¬ 
quire separate renormalization. Nevertheless, it has been 
shown that FDH is a consistent regularization scheme 
to NNLO [l60l ]. The presence of these extra degrees of 
freedom in the FDH scheme is quite inconvenient in the 
formal construction of SCET, especially when working 
to subleading power. Because of this fact, and because 
most SCET calculations are performed in CDR/HV, our 
discussion of SCET schemes will focus on regularization 
schemes where the dimension of the gluon field and the 
momentum space are analytically continued in the same 
manner. We will also discuss how full-theory helicity am¬ 
plitudes in the FDH scheme are converted to CDR/HV 
for the purposes of matching to SCET. 

We will now describe how helicity amplitudes in the 
FDH scheme can be converted to CDR/HV. To get a 
finite correction from the O(e) part of the gluon polar¬ 
ization requires a factor from either ultraviolet (UV) or 
infrared (IR) 1/e divergences. Although the regulariza¬ 
tion of UV and IR divergences is coupled in pure dimen¬ 
sional regularization schemes by use of a common e, they 


Recently while this paper was being finalized, a calculation of 
the inclusive jet and soft functions in both FD H and dimensional 
reduction (DRED) [l59l| appeared in Ref. |160| . The conclusions 
of this section agree with their study of the regularization scheme 
dependence of QCD amplitudes. 
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can in principle be separately regulated, and we discuss 
their role in the scheme conversion separately below. 

When matching to SCET, the UV regulators in the 
full and effective theory need not be equal. Indeed, 
the effective theory does not reproduce the UV of the 
full theory. In massless QCD, scheme dependence due 
to the UV divergences only affects the coupling con¬ 
stant through virtual (internal) gluons. Therefore, the 
CDR and HV schemes have the same standard MS 
coupling, a s (/i), while FDH has a different coupling, 
a™ H (g). The conversion between these couplings is 
achieved by a pert urbatively calculable shift, known to 
two loops [53, ll58l il63| 


a s {n) = a s {n) 


1 + 


Ca a 8 {ij) 
3 47t 


+ (y c a - 2C F T F n,f S J 


47T 


( 100 ) 


This replacement rule for the coupling captures the ef¬ 
fect of the scheme choice from UV divergences. One can 
therefore perform a matching calculation, treating a s in 
the full and effective theories as independent parameters 
that can be defined in different schemes. A conversion 
between schemes can then be used to ensure that the 
matching coefficients are written entirely in terms of a s 
defined in one scheme, for example using Eq. (1 1001) . The 
issue of UV regularization is therefore simple to handle 
in the matching. 

The structure of 1/e 2 and 1/e IR divergences in one- 
loop QCD amplitudes is well known, and allows one 
to determine their effect on converting amplitudes from 
FDH to CDR/HV. For a QCD amplitude involving n q 
(anti)quarks and n g gluons the FDH and HV one-loop 
amplitudes A ^ are related by [H, Il55( | 

A 1 J = 4 1 A-g(fc F + ^)>im, (101) 

where denotes the tree-level amplitude, and the pre¬ 
cise scheme of the a s entering here is a two-loop effect. At 
one loop, the FDH scheme can therefore be consistently 
used when calculating full-theory helicity amplitudes and 
results can easily be converted to HV with Eqs. (1 1001) and 
(I101D for use in SCET Wilson coefficients. 

We will now compare CDR and HV schemes for SCET 
calculations and the construction of the operator basis. 
In the HV scheme, all external polarizations are 4 di¬ 
mensional, so that one can use a basis of helicity oper¬ 
ators, as was constructed in Sec. m However, in CDR 
external polarizations are d dimensional, with the limit 
d —> 4 taken. In particular, this implies that one must 
work with d— 2 gluon polarizations at intermediate steps, 
potentially allowing for the presence of evanescent oper¬ 
ators corresponding to operators involving the additional 
components of the gluon field, so-called e-helicities. How¬ 
ever, we will now argue that there is no real distinction 
between the two schemes, and that one does not need to 
consider evanescent operators in SCET at leading power. 


First consider the Wilson coefficients and matching. 
In the case of CDR, the operator basis must be extended 
to include operators involving the e-helicities. However, 
their presence does not affect the matching coefficients for 
operators with physical helicities, since they do not con¬ 
tribute at tree level and all loop corrections are scaleless 
and vanish. Additionally, in Sec. IVIIII we will discuss 
the fact that the SCET renormalization of the opera¬ 
tors is spin independent at leading power, and therefore 
there is no mixing under renormalization group evolu¬ 
tion between the physical and evanescent operators. For 
the beam and jet functions, azimuthal symmetry implies 
that the difference between a field with 2 or 2 — 2e po¬ 
larizations is simply an overall factor of 1 — e and thus 
can be easily taken into account. The independence of 
the soft function to the differences in the CDR/HV reg¬ 
ularization schemes follows from the insensitivity of the 
soft emission to the polarization of the radiating parton, 
which is made manifest by the SCET Lagrangian and the 
fact that the soft function can be written as a matrix ele¬ 
ment of Wilson lines. Thus there is no difference between 
CDR and HV and the helicity operator basis suffices. 


V. HIGGS + JETS 

In this section, we consider the production of an on- 
shell Higgs + jets. We give the helicity operator basis and 
matching relations for H- 1-0,1, 2 jets, and the correspond¬ 
ing helicity amplitudes are collected in Appendix Q3] 


A. H + 0 Jets 

The ggH and qqH processes contribute to the H + 0 
jets process. For qqH, the scalar current in Eq. (1551) is 
required, and the helicity operator basis is given by 

of = j°£ 0 h 3 , 

Of = (J')f 0 H3, ( 102 ) 

with the unique color structure 

T aB = (5 aP ) . (103) 

These operators are relevant when considering Higgs de¬ 
cays to massive quarks, for example H —> bb. However, 
we will not consider this case further since for Higgs pro¬ 
duction the bbH and ttH contributions are much smaller 
than the dominant gluon-fusion hard scattering process. 
For ggH , the basis of helicity operators is given by 

Of + = l -B" + BfH :i , 

o a _b_ = 1 B a_ B b_ Hs ( 104 ) 

The operator O _|_is not allowed by angular momentum 

conservation. Similar helicity operators, extended to in¬ 
clude the decay of the Higgs, were used in Ref. [l64l |. 
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There is again a unique color structure for this process, 

fab = _ ( 10 5) 

Writing the QCD helicity amplitudes as 

A(g l92 H 3 ) =iS a ^ A(l,2; 3 H ), 


(106) 


the Wilson coefficients for ggH are given by 


C++(pi,P2',P3) = A fin (l + ,2 + ;3 ff ), 
C__(p 1 ,p 2 ;p 3 ) = AU I", 2"; 3 H ). (107) 


The subscript “fin” in Eq. (11071) denotes the IR-finite 
part of the helicity amplitudes, as discussed in Sec. [TV] 
Note that the two amplitudes appearing in Eq. (11071) are 
related by parity. The results for the gluon amplitudes 
up to NNLO are given in Appendix IB II They corre¬ 
spond to the usual gluon-fusion process, where the Higgs 
couples to a (top) quark loop at leading order. The LO 
amplitude including the dependence on the mass of the 
quark running in the loop is well known. The NLO am¬ 
plitudes are also known including the full quark-mass 
dependence Jni5Hl69| . while the NNLO fl7(m72l l and 
N 3 LO 0, [nT amplitudes are known in an expansion in 
m H /m t . 


B. H + 1 Jet 

The gqqH and gggH processes contribute to the H+l 
jet process. For gqq , the basis of helicity operators is 
given by 


n“ a ^ _ I3 a T a P IT, 

— °1+ J 23+ ) 

— B a T a P H, 

Ly _(+) ~ °1- J 23+ > 

/-> a &P _ tra ra/3 tt 

0 + (_) - £>1+ J-23- I 

0““f) = B a x _ j°£_ H a . (108) 

Note that we consider only QCD corrections to the ggH 
process, so the qq pair is described by J t j± ■ For ggg , the 
helicity operator basis is 

0?Z + = ±B a 1+ B b 2+ B c 3+ H 4 , 

0 ^_= l -B a 1+ B b 2+ B c 3 _H A , 

0--+ = B"_ B b _ B' i+ H 4 , 

O ab l_ = i B a x _ B\_ B c 3 _ H 4 . (109) 

For both cases the color space is one dimensional and we 
use the respective color structures as basis elements 

faa£ = ( T a_) , fabc = ^jabcj _ ( 110 ) 


In principle, there could be another independent color 
structure, d abc , for gggH. The gggH operators transform 
under charge conjugation as 

col b Xx 3 (pi,P2,h;P4)T abc c 

= ~O cb Xx 3 (P^P2,PS-,P4)f abc . (Ill) 

Charge conjugation invariance of QCD thus leads to 

Cx%Xa(Pl,P2,h)P4,) = -CaiA 2 A 3 (Pi>P2,P3;P4) , (H2) 

which implies that the d abc color structure cannot arise 
to all orders in perturbation theory, so it suffices to con¬ 
sider i f abc as in Eq. (jllOjl . This also means that the 
d abc color structure cannot be generated by mixing un¬ 
der renormalization group evolution, which will be seen 
explicitly in Eq. (12591) . 

Using Eq. (I110D . we write the QCD helicity amplitudes 
as 


-4(513233^4) = i (i/“ l0a “ s ) A(l, 2,3; A H ), 
A( 9l q 2 q 3 H 4 ) = i2£ a3 24(l;2„3 9 ;4 ff ) . (113) 

The Wilson coefficients for gqqH are then given by 

C+(+){PBP2,P3]Pi) = A fin (l+;2+,3^;4 ff ), 
C_ i+) (pBp 2 ,p 3 ;p 4 ) = A fin (l“;2+,3^;4 ff ), 

C+(-)(pi;p2,P3;p4) = C'+(+)(p 1 ;P3,P2;P4), 
C-(-){PBP2iP3\P±) = C-(+)(pBp 3 ,p 2 \p 4 ) , (114) 


where the last two coefficients follow from charge conju¬ 
gation invariance. The Wilson coefficients for gggH are 
given by 


c +++ 

<?++_ 

c__ + 


(Pi,P2,P3',P4) = A fln (l+,2+,3 + ;4 ff ), 
{Pi,P2,P 3 ;P4) = A fln ( 1 + ,2 + ,3~;4 h ) , 


(Pi,P2,P3',Pa) = C ++ _{p 1 ,p 2 ,p 3 ;p 4 ) 




C -(pi, P2 , P3 ; Pi) = C+++ {pi ,P2,P3 \Pa) 


<••>«•[••] 

(115) 


where the last two relations follow from parity invariance. 
As before, the subscript “fin” in Eqs. (11141) and (11151) de¬ 
notes the finite part of the IR divergent amplitudes. The 
NLO helicity amplitudes were calculated in Ref. j32|, and 
are given in Appendix IB 21 and the NNLO helicity am¬ 
plitudes were calculated in Ref. 173]- Both calculations 
were performed in the m t —> oo limit. At NLO, the first 
corrections in m 2 H /m 2 were obtained in Ref. 174]. 


C. H + 2 Jets 

For H + 2 jets, the qqq'q'H , qqqqH , ggqqH, and 
ggggH processes contribute, each of which we discuss 
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in turn. Again, we consider only QCD corrections to the 
ggH process, so qq pairs are described by the helicity 
currents Jij±■ The LO helicity amplitudes for H + 2 jets 
in the m t -» oo limit were calculated in Refs. ^3, 1175f 
and are collected in Appendix IB 31 for each channel. The 
LO amplitudes including the mt dependence were calcu¬ 
lated in jl76l ] (but explicit expressions for ggggH were 
not given due to their length). The NLO helicity ampli¬ 
tudes were computed in Refs. EH HI EMUS- 


The QCD helicity amplitudes for qq q'q'H can be color 
decomposed in the basis of Eq. (11181) as 


-^(91929394-^5) = 2 lZ> 


+ 


1 

N 


^aiQ 4 ^a3a2^(^5 2q, 3q ', 4 q> ,5//) 

*^aia 2 ^a3a 4 71(lq, 2^; 3gq 5//) , 

(119) 


1. qq q'q'H and qqqqH 


For the case of distinct quark flavors, qq q'q'H, the 
helicity basis consists of four independent operators, 


O 

O 

o 

o 


a(3'i8 

(+;+) 

(+;-) 

(-;+) 

afl'yS 


t«/3 

J q 12 + 

r«/3 
J q 12 + 

T a/3 
J q 12 - 

ja/3 
J q 12 - 


Jq> 34+ i 

4 S . 34-^, 

^34+^5, 


(116) 


where the additional labels on the quark currents indi¬ 
cate the quark flavors. For the case of identical quark 
flavors, qq qqH , the basis only has three independent he¬ 
licity operators, 


n ot0j5 _ 1 ja/3 7-76 tt 
Ly (++) — 4 J 12 + ^34+ n 5 , 

n a0jS _ ra/3 j0S jj 
Ly (H—) — '-' 12 + ^34- 11 5 ; 

Of^ = (117) 

since both quark currents have the same flavor. In both 
cases we use the color basis 

T"™ = 2T F (S aS 5 lB , 8 ajS <J 7 j) . (118) 


where we have included a factor of 1 /N for convenience. 
The amplitude vanishes when the quark and antiquark 
of the same flavor have the same helicity, in accordance 
with the fact that the operators of Eq. (11161) provide a 
complete basis of helicity operators. For identical quark 
flavors, qq qqH, the amplitudes can be obtained from the 
qq q'q'H amplitudes using the relation 


- 4 ( 91929394 ^ 5 ) = ■A(qiq 2 q' 3 q' A H 5 ) - A(qiq4q' 3 q' 2 H 5 ) . 

( 120 ) 

The Wilson coefficients for qq q'q' H are then given by 


3+ 


;+)(Pl,P2;P3,P4;P5) 


C'(+;-)(Pl,P2;P3,P4;P5) 


C(_ ;+ )(pi,p 2 ;p 3 ,p 4 ;p5) 

C' ( _.._)(pi,p 2 ;P3,P4;P5) 


( A fin (l+,2r;3+,4-;5 ff ) \ 
^^( 1 + 2553 +, 4^:5^’ 
( A fin (l+,2,-;3-„4+;5 ff ) \ 
Vii? fi n(l+,2 9 -; 3 -„4+;5 ff ); 5 

C (+; _)(p 2 ,pi;p 4 ,P3;P5), 
C'( + . + )(p 2 ,pi;p4,p 3 ;P5), 

( 121 ) 


and for qq qqH they are given in terms of the amplitudes 
Afi n and B Rn for qq q'q’H by 


( 7(++)(pi J P2;P3,P4;p5) 


C'(+- ) (pi,p 2 ;p 3 ,p 4 ;p5) 


( ^4fin(l+, 2 g ; 3+ 4 g ; 5 H ) - ^B fi „(l+, 4 ? -; 3+ 2 g ; 5 H )\ 
\jjBfi n (l+,2g-,3+Ag-,5 H )-A Rn (l+,4g-,3+,2g-,5 H )J ’ 

( A fin (l+2r ; 3-,4+;5 ff ) \ 

\±BU 1+ 2r;3-,4+;5 ff ) y / ’ 


C( -)(P1,P2;P3,P4;P5) = C ! (++) (p 2 ,pi;p4,P3;P5) • 


( 122 ) 


The relations for C(_ ; ±) and C(_) follow from charge 

conjugation invariance. Note that there is no exchange 

term for C(_|_), since the amplitude vanishes when the 

quark and antiquark of the same flavor have the same 
helicity (both + or both —). Also, recall that the sym¬ 
metry factors of 1/4 in Eq. (11171) already take care of 
the interchange of identical (anti)quarks, so there are no 


additional symmetry factors needed for C( ++ y Explicit 
expressions for the required amplitudes at tree level are 
given in Appendix IB 3 al 
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2 . ggqqH 


3■ ggggH 


For ggqqH , the helicity basis consists of a total of six 
independent operators, 

n abaP _ 1 na nb ra/3 tt 
°++(+) — 9 °1+ °2+ J 34+ i 

t+aba0 _ na nb raft tt 

U +-(+) — °1+ C 2- J 34+ , 


o 


ab a(3 


i . 


= -B a x _B\_ &H S , 


—(+) “ 2 


r^abotfi _ * na nb jap tt 

u ++(~) - 2 Cl + " 2 + J 34- ' Hs ’ 

O ab ^_ ) = lB a 1 _B b 2 _J^-H 5 . 


1 


(123) 


We use the color basis already given in Eq. m, 

faba p = ((T“T% , (T b T a ) a p , tr [T a T b ] S a0 ) . (124) 

Using Eq. (11241) . the color decomposition of the QCD 
helicity amplitudes into partial amplitudes is 

-^(9152 9394^5) 

= i E [ r ' {1)r ’ P) ] « 3 a 4 A ( ct (1 )’ f7 ( 2 ); 3 9 ,4g; 5 H ) 

creS 2 

+ i tr[T ai T a2 ] 5 Q 3 a 4 5(1, 2 ; 3„ 4 5 ; 5h) ■ (125) 

The 1? amplitudes vanish at tree level. From Eq. (11251) 
we can read off the Wilson coefficients, 

/A fin (l+,2-;3+,4^;5 H p 
C+_^+ ) {p 1 ,p 2 ;p 3 ,p 4 -,p 5 ) = A fin (2-, l+;3+,4^;5 ff ) 

\Bfin(l + ,2 _ ;3+,4^;5_f/) y 

/A fin (l+,2+;3+,4,j-;5 ff )^ 
C ++ ( +) (pi,P2;P3,P4;P5) = A fin (2+,l+;3+,4,j ; 5 ff ) 

Vl? fin (l+,2+;3+,4r;5 ff ) y 

/A fin (l-,2-;3+ 4 ? r ;5ff )\ 

C —( + )(pi,p 2 ;P3,P4;P5) = Afin(2", l~;3+,4^;5 ff ) . 

\B fi „(l- 2-;3+4r;5*)/ 
(126) 

The Wilson coefficients of the last three operators in 
Eq. (11231) are obtained by charge conjugation as discussed 
in Sec. IIII El Under charge conjugation, the operators 
transform as 

c °K±) (Pi, P 2 ; P3 , Pi ; p 5 ) T ab ^ c 

= -O b x a i ^ iT) (p 1 ,p 2 ;p A ,p 3 ;p 5 )f ab ^, (127) 
so charge conjugation invariance of QCD implies 
C'a 1 a 2 (-)(pi,P2;p3,P4;p5) = UC'a 1 a 2 (+)(pi,P2;p4,P3;p5) 


0-10 
with V = | — 1 0 0 

0 0 - 1 , 


(128) 


Explicit expressions for the required amplitudes at tree 
level are given in Appendix IB 3 bl 


For ggggH, the helicity basis consists of five indepen¬ 
dent operators, 


O 


1 


abed _ 

H—I—I—b ^l+^2+^3+^4+ 


o 


abed 


= -B a 1+ B b 2+ B c 3+ B a 4+ H 5 , 
1 


+++. = ^B a 1+ B D 2+ B c 3+ B a 4 _H 5 , 


O 


abed 
++- 


= - B a l+ B° 2+ B c 3 ^B a 4 _ H 5 , 


1 


O abc _ d _+ = B a 4 _B b _B c 3 _B d 4+ H 5 . 
O ab f__ = ^ Bl_B b _Bl_Bi_ H 5 . 


We use the basis of color structures 


(129) 


rj-iabcd _ 


2-2 T f 


/tr [abed] + tr [dc 6 a] Y 
tr [aedb] + tr \bdca\ 
tr [adbc] + tr [cbda] 
2 tr[afe] tr[cd\ 

2 tr [ac] tr[d&] 

V 2 tr[ad] tr[ 6 c] ) 


(130) 


where we have used the shorthand notation 

tr[o 6 ] = tr[T a T b ], tr[a 6 cd] = tr[T a T b T c T d ]. (131) 

Note that the three independent color structures with a 
minus sign instead of the plus sign in the first three lines 
in Eq. (11301) can be eliminated using charge conjugation 
invariance, see Sec. IVII A 31 

The color decomposition of the QCD helicity ampli¬ 
tudes into partial amplitudes using the color basis in 
Eq. (11301) is 


^(#l£2#3£4#5) = -^jT 


^ ^ tr[a a (i) a (J (2) g<t(3) &<t(4) ] 
L <t£*S4 / Z 4 

X A(er(l),cr(2),cr(3),cr(4);5 ff ) 

T ) ( tr[a 0 -(i)a CT ( 2 )]tr[a (T ( 3 )a (T ( 4 )] 

X u(cr(l),cr(2),cr(3),cr(4);5 ff ) 

(132) 

where the B amplitudes vanish at tree level. From 
Eq. (11321) we obtain the Wilson coefficients, 


/2A fin (l+,2+,3-,4-;5 ff )\ 

2A fin (l+,3-,4-,2+;5 ff ) 

A ( .x 2A fin (l+,4-,2+,3-;5 ff ) 

C++ — (P1,P2,P 3 ,P4;P5) = 5 g n (l + , 2 + , 3 _ ,4 _ ;5#) 

S fin (l+,3-,4-,2+;5 H ) 

Vu fin (l+,4-,2+,3-;5 ff )/ 
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C+++-(pi,P2,P3,P4-,P5) 


C++++(P\,P2,P3,P4-,P5) 


/2A fin (l+,2+,3+,4-;5 ff )\ 
2A fin (l+,3+,4-,2+;5 ff ) 
2A fin (l+,4-,2+,3+;5 ff ) 
S fin (l+, 2+,3+,4-;5 h ) 
B fln (l+,3+,4-,2+;5 ff ) 
VB fln (l+,4-,2+,3+;5 ff )/ 

/2A fin (l+, 2+ 3+ 4+;5 ff )\ 
2A fin (l+,3+, 4+,2+;5 ff ) 
2A fin (l+,4+,2+,3+;5 ff ) 
R fln (l+,2+,3+,4+;5 ff ) 
R fln (l+,3+,4+,2+;5 ff ) 

V Sfi n (l+, 4 + , 2+, 3 + ; 5h) ) 


c_ 

1 

1 

+ 

■ ■ 1 P 5 ) 

— ^+++-(pii • 

■ ■;P5) 

c _ 

— (pi,. 

■ •;P5) 

= (Pi > • 

• ■ ;P5) 


(■■>«■[••] ’ 

<•■>++[■■] ' 
(133) 


The last two coefficients follow from parity invariance. 
The factors of two in the first three entries of the coef¬ 
ficients come from combining the two traces in the first 
three entries in Eq. (ESP using charge conjugation invari¬ 
ance. Because of the cyclic symmetry of the traces, the 
partial amplitudes are invariant under the corresponding 
cyclic permutations of their first four arguments, which 
means that most of the amplitudes in Eq. (11331) are not 
independent. Explicit expressions for the necessary am¬ 
plitudes at tree level are given in Appendix IB 3 cl 


VI. VECTOR BOSON + JETS 

In this section, we give the helicity operator basis and 
the corresponding matching for the production of a 7 , 
Z, or W vector boson in association with up to two jets. 
The corresponding helicity amplitudes are collected in 
Appendix O 

We work at tree level in the electroweak coupling and 
consider only QCD corrections, so any external qq pairs 
are described by the helicity vector currents Jij± in 
Eq. @7|). We always include the subsequent leptonic de¬ 
cays 7 /Z —*• tl , W ± —> vl/tv. In the following, for 7 jZ 
processes, i stands for any charged lepton or neutrino fla¬ 
vor, and q stands for any quark flavor. For W processes, 
we use £ to denote any charged lepton flavor and v the 
corresponding neutrino flavor. Similarly, we use u and 
d to denote any up-type or down-type quark flavor (i.e. 
not necessarily first generation quarks only). 

The operators in the helicity bases satisfy the trans¬ 
formation properties under C and P as discussed in 
Sec. IIII El However, the weak couplings in the ampli¬ 
tudes explicitly violate C and P. Therefore, to utilize the 
C and P transformations of the operators and minimize 
the number of required amplitudes and Wilson coeffi¬ 
cients, it is useful to separate the weak couplings from 
the amplitudes. 

We define P z and Pw as the ratios of the Z and W 


propagators to the photon propagator, 

p z,w(s) =- 2 --' ( 134 ) 

s ~ m z,w + 11 z,wmz,w 


The left- and right-handed couplings vl,r of a particle 
to the Z boson are, as usual, 


*l = 


2 Tj - 2 Q i sin 2 6 W 

sin(2f?w) 


2 Q i sin 2 0 W 
sin (2 6w) 


(135) 


where TJ is the third component of weak isospin, Q 1 is 
the electromagnetic charge in units of \e\, and 9\y is the 
weak mixing angle. 

The "f/Z amplitudes can then be decomposed as 


A(- ■ ■ i£) 

= e 2 | [QV + vl R vl R P z (s a )\A q (- • • it) 

n f 




l 'L,R ' 


- p z{su) Av(---ii) 


U L,R 


sin(2$vv) 


Pz(Stf) A a (’ 


(136) 


Here, A q corresponds to the usual contribution where the 
vector boson couples directly to the external quark line 
with flavor q. (There is one such contribution for each 
external qq pair, and this contribution is absent for pure 
gluonic amplitudes like gggZ.) For A v , the 7 /Z couples 
to an internal quark loop through a vector current and 
the sum runs over all considered internal quark flavors. 
For A a , the Z boson couples to an internal quark loop 
through the axial-vector current. This means that when 
using parity and charge conjugation we have to include 
an additional relative minus sign for this contribution. 
We have also made the assumption in Eq. (I136D that all 
quarks, except for the top, are massless. Since A a van¬ 
ishes when summed over a massless isodoublet, this has 
the consequence that only the b, t isodoublet contributes 
to A a , hence the lack of sum over flavors. We have made 
this simplification following the one-loop calculation of 
Ref. fl4j |. which calculated the amplitude in an expan¬ 
sion in 1 /to 2 , assuming all other kinematic invariants to 
be smaller than the top mass. From the point of view 
of constructing a basis these assumptions are trivial to 
relax. 

The amplitudes can be written as 


e 2 V j 

A{- ■ ■ t~v + ) = - y l - Pw{siv) A q {■ ■ ■ £~D + ), 

2 sin 6w 

A- ■ ■ v-t*) = f V k Pw(s vi ) A q (--- v-~t), (137) 
2 sin 9 W 

where V u d is the appropriate CKM-matrix element. The 
A q amplitudes are the same in Eqs. (1 1361) and (11371) . 
since all electroweak couplings have been extracted, but 
we have explicitly included the helicity labels (not to be 
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mistaken as charge labels) to emphasize that these are 
the only possible helicities. The analogs of A v and A„ 
do not exist for W production. 

We note again that Eqs. (11361) and (11371) hold at tree 
level in the electroweak coupling, which is what we con¬ 
sider in this paper. At this level, the leptons always 
couple to the vector boson through the currents [see 
Eq. JAH] 


(±) = (pfFlVIWT) • (138) 

This allows us to obtain the Wilson coefficients for oppo¬ 
site lepton helicities simply by interchanging the lepton 
momenta. 


A. V + 0 Jets 

For 'y/Z + 0 jets, the partonic process is qqll- and the 
basis of helicity operators is 

34±, 

0^. ±) =J%_Jt 34±- (139) 

In principle, the process ggtl is allowed through the ax¬ 
ial anomaly, but its contribution vanishes because in the 
matching calculation the gluons are taken to be on shell, 
and we neglect lepton masses. 

For W T + 0 jets, the partonic processes are vMv and 
duvl, respectively. Since the W only couples to left- 
handed fields, the helicity basis simplifies to 

n &P _ r<5/3 t 

^(W~) — j ud 12- 34- i 

°tw +) = Jlu 12-^34-. (140) 

Here, we have explicitly written out the flavor structure 
of the currents. However, we use the shorthand subscript 
(W^) on the operators and Wilson coefficients, since we 
will not focus any further on the flavor structure. In 
an explicit calculation, one must of course sum over all 
relevant flavor combinations. 

The unique color structure for V + 0 jets is 

T a ~ 0 = (6 a p) , (141) 

and extracting it from the amplitudes, we have 

Aq,v : a(qiq2^3^4:') 1 (lg, , 3f, 4^) . (142) 

Here, A v and A a first appear at two loops. In addition, 
A a is proportional to the top and bottom mass splitting 
due to isodoublet cancellations. It drops out when both 
top and bottom are treated as massless (e.g., when the 
matching scale is much larger than the top mass). 

We use the same electroweak decomposition as in 
Eqs. (I136D and (11371) to write the Wilson coefficients. For 
'y/Z + 0 jets, we have 


= e 2 | [QV + vivlPz(s 3 4)\ C q(Xg;Xe) (...) 

n f • 

+E [qv+< ^^ Pz(s 34 )] c v(Xq . M) {. ..) 


2=1 

+ ^fc) Pz(S34)(?a(A9;A^)( ■■■ ) } , (143) 

where the weak couplings are determined by the helicity 
labels of the quark and lepton currents, 

v+ = v e R , v e _ = v e L , v\ = v q R , v q _ = v q L . (144) 

For W + 0 jets, we simply have 


C(W-)(Pl,P2-,P3,p4) 

C(W+){PUP2-,P3,P4) 


e 2 Vud 

2 sin 2 9w 
2 sin 2 9w 


PW(S 3 4) C,(_;_)(. ..), 
PW(S 34 ) C ? (_;_)(. ..). 

(145) 


In all cases, the momentum arguments on the right-hand 
side are the same as on the left-hand side. Note that the 
C'q(_._) coefficient is the same in all cases. The Wilson 
coefficients are given by 

C X ( + -+)(P1,P2;P3,P4 ) = A Xt fin(lg ,2j ;3|,4p , 

C x (+.-)(Pl,P2-,P3,P4:) = ^(+;+)(Pl,P2;P4,P3) , 
Cq,v(—\±){Pli P2i P3i P4) ^q,v(-\-;4z){P2j Pi: P3i P 4 ) j 

Ca(-;±)(Pl,P2;P3,P4) = -C a ( +; ±)(p2,Pi;P3,P4) , 

(146) 

where x = q,v,a and as discussed in Sec. m the sub¬ 
script “fin” denotes the IR-finite part of the helicity am¬ 
plitudes. The second relation follows from Eq. (I138D . 
The last two relations follow from charge conjugation in¬ 
variance. At tree level and one loop only C q receives a 
nonvanishing contribution. The A q amplitude is given in 
Appendix 1C II 


B. V + 1 Jet 


1. gqqV 


For 7 /Z + 1 jet, the partonic process is gqqii, and the 
basis of helicity operators is 

< 5/3 _ fia jdt(3 j 

^+(+ ; ±) — ^1+ Jq 23+ ^45± 5 

^+(- ; ±) — ®l+ Jq 23- ^45± 7 

O a _l f ;±) = J ? t 3+ -W , 

O a _^. ±) =Bt_J%_J t 4B±. (147) 

Fop W T + 1 jet, the partonic processes are gudtv and 
gduui , respectively, and the helicity operator basis is 


n a a P 
u ± (W~) 


r>a ra.fi j 
°1± J ud 23- J iv45- ) 


C( Xq - Xl ){pi,P2\P3,P4) 
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0“^ +) =^ ± j|f 23 _J^45-. (148) 

The unique color structure for gqqV is 

faa p = (T^) , (149) 

and extracting it from each of the amplitudes, we have 

A x ( gi q 2 q 3 hh) = iZ£ aa ^(!; 2 „ 3 ? -; 4 ,, 5 ,-), (150) 

where the subscript x stands for one of q , v, a. 

As for V + 0 jets, we write the Wilson coefficients using 
the electroweak decomposition in Eqs. (I136D and (11371) . 
For 7 /Z + 1 jet, we have 

<?A(A,;A/)(Pi;P2,P3;P4,P5) 

= e 2 { [Q l Q q + v{y Xq Pz(s 45 )]C gHXq;Xe) (...) 

n f 






v'l + Vr 


i =1 


Pz(s 45 ) C vX ( Xq;Xt )(...) 


+ 


—^-yPz( S 45 )C aA(Ag;A£) (...) 


(151) 


where the weak couplings are determined by the helicity 
labels of the quark and lepton currents, 


vy = v 


R > 


V_ = V 


L > 


v 9 + =vl, vl = v q L . (152) 


For W + 1 jet, we have 
_ e 2 F (t) 

Cx(WT)(. • •) = 0 . 2™? iV(«4s) C 9A( _ ; _)(. . .) . (153) 

2 sin % 

The Wilson coefficients are given by 


2. gggV 

The partonic process gggll first appears at one loop, 
and thus contributes only at relative 0(a 3 ) to 7 /Z+l jet. 
Nevertheless, for the sake of completeness (and curiosity) 
we briefly discuss it here. The helicity operator basis is 

1 


OTl^ = ^ B a 1+ B b 2+ B c z+ J u5 ± 


-\abc 

W++(±) - 31 ^1+ 

i B'!', B b m_ J t 


++-(±) = 2 l ~’ 1 + U2 + ° 3 - J£45± ' 
1 


Ofl_ (±) = -B! + B b _B^_ ,W 


O a _ b l_ {±) = | B\_ B\_ B%_ J e45 ± • (156) 

The color space is two dimensional. We use the basis 
T abc = (i f abc , d abc ^j , (157) 

in terms of which we can write the gggll amplitudes as 

A v ( 9l g 2 53 44) = id aia2a3 A V {1, 2 ,3; 4,, 5*), 
A a (gig 2 g 3 hl 5 ) = i (if m3 )A a (l, 2 ,3; 4*, 5g ). (158) 

We will justify shortly that to all orders, only a single 
color structure appears for each of A Vl A a ■ This process 
can only occur via a closed quark loop, so there is no 
A q contribution. The gggV operators transform under 
charge conjugation as 


C °Ai a 2 Aa (±) (Pi i .P 2 , P3; P4, P5) T abc C 

= Of 1 a X2X3(T) (p 1 ,p 2 ,p 3 ;P5,P4) T abc . 


(159) 


C x+ (+-+)(pi-,P2,p 3 -,P4,P5) = A X)fin (l + ;2+,3^;4|,5jT), 

C xX (+--){pi;P2,P 3 ;P4,P5) = ^xX(+;+){pl'>P2,P 3 ',P5,Pi) , 
C g ,v\(—,± )(Pi;P2,P3;P4,P5) = — (?g,w A (+;±)(Pi;P3)P2;P4)P5) 
C a \(--,±){Pl\P2,P 3 ;P4,P5) = C aX{+ . ±) {pi;p 3 ,P2-,P4,P5) ■ 

(154) 

The second relation follows from Eq. (11381) . and the last 
two relations follow from charge conjugation invariance. 
The Wilson coefficients with a negative helicity gluon fol¬ 
low from parity invariance, 


C q , v - (+ ; ±) (pi ; P2 , P 3 ; Pi , P 5 ) 

= C qtV+ (-; T )(pi; P2, p 3 ; p 4 , P5) 

Ca-(+-,±)(Pl',P2,P 3 -,P4,P5) 

= -C a+{ _. T) (p 1 ;p 2 ,p 3 ;p 4 ,p 5 ) 


(155) 


The helicity amp litudes for gqqll were calculated in 
Ref. 0, Il81| . We provide the tree-level and one- 

loop results in Appendix 1C 21 The two-loop amplitudes 
were computed in Refs. ll82Ul83j . 


Charge conjugation invariance of QCD thus leads to 

C v\ix 2 X 3 (±)(P^P2,P3;P4,P5) 

= CZ b XlX2X3(T )(pi,P2,P 3 ;P5,P4) 

= C^AiAaA 3 (±)(PuP2,Ps;P4,P5) , (160) 

where we used Eq. (11381) in the last line. This implies 
that to all orders in the strong coupling, only the fully 
symmetric color structure d abc can contribute to A v and 
C v . For C a the same relation holds but with an addi¬ 
tional minus sign on the right-hand side due to the weak 
axial-vector coupling in A a - This implies that for A a 
and C a only the fully antisymmetric color structure i f abc 
contributes, as given in Eq. (11581) . 

We decompose the gggli Wilson coefficients as 


C. 


Ai A 2 A3(A^) 


= e 


£[qV 


+ <- k -r JL PzM 


a 


vAi A2 As(A^) 


2=1 

V 


—j^PzMC aXlXaX3lXl) j, (161) 
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where 


v+ = v R , 


V_ = V L 


(162) 


and we have 


^vXi A2A 3 (+) (Pi 5 P2i P3i P4:i Pb) 

^ > fi„(l Al ,2 Aa ,3^;4+ 5jT), 

^-'aAiA 2 A 3 (+) (Pi7 P 2 i PS] P§) 

1 Ai 0^2 *^A 3 

0 

<A,aAiA 2 A 3 (-)(ihjP2,P3;P4,P5) 

P'v.aXi A 2 A 3 (+) (lT 5 -P2 7 P3i P5 7 P 4 ) ■ 


A a , Rn ( 1 A S2 A %3^;4+ 57 ) 


(163) 


For brevity, we have not written out the various gluon he- 
licity combinations. The one-loop amplitudes for gggZ 
were calculated in Ref. HP, and the two-loop ampli¬ 
tudes were computed in Ref. [185] ]. Since their contribu¬ 
tion is very small we do not repeat them here. 


C. V + 2 Jets 


Here we consider the processes q'q'qqV , qqqqV, and 
gg qq V. The ggggV process is allowed as well, but only 
arises at one loop, so we do not explicitly consider here. It 
can be treated similarly to gggV , but using the gggg color 
basis analogous to that for ggggH given in Eq. (I130D . 

The NLO helici ty a mplitudes for R + 2 jets were calcu¬ 
lated in Refs. fl4l Il86l] assuming that all kinematic scales 
are smaller than the top mass rrit and including the 1 / m\ 
corrections. We give the full expressions for the LO re¬ 
sults in Appendix 1C 31 Since the NLO results are rather 
long, we do not repeat them, but we show how to convert 
the results of Refs, [H, ll86j | to our notation. 


1. q'q'qqV and qqqqV 
For q'q'qqil, the helicity operator basis is 


Of-%) = \ J %~ J %~ • (165) 

For W + 2 jets, the corresponding partonic processes are 
qqudiD and qqduizi , and the helicity operator basis is 


^\a/3-y8 to;/3 7-7 8 j 

U (±]W~) ~ J q 12± J ud3A- J tv 56- > 

n a0i5 _ tS/3 t7<5 t 

U (±\W+) ~ J q 12+ J du 34- J ^56- • 

(166) 

We use the color basis 


f^ = 2 TF ^ a - s 6 J - 0 ,5 a pS 1 s)- 

(167) 

For distinct quark flavors, the color decomposition of the 
amplitudes in this basis is 

Ax(qW 2 Q3Q^^6) 

— 2T F i ^Qia 4 ^a 3 a 2 -^xip-q' j 2 ^/ 5 3^, 4^, 5^, 6 j?) 

(168) 

H - 2TV4 S 0 il 0 i 2 ^ 0 : 30:4 (lq/, 2g/, 3^, 4gr, 5^, 6 ^). 

For identical quark flavors the amplitudes can 
tained from the distinct flavor amplitudes using 

be ob- 

Ac id! 92930444) = A x (^2Q3Qihh) 

- A x (q[ ^935244) , 

(169) 


where it is to be understood that the electroweak cou¬ 
plings of q' must also be replaced by those of q. 

Writing the Wilson coefficients in the decomposition 
in Eqs. (1 13611 and (113711 , we have for the q'q'qqti channel 


C(\ q , ; \ q; \ e )(pi,P2;P3,PVP5,P6) 

= e 2 | [Q e Q q + v{ ( v q xq Pz{s m )\ 

X PA P3,PVP5,P6) 

+ [Q e Q q ' +v{ e vi q ,Pz(s 56)] 

X A,; A 9 ,; A,) (P3, P4; Pi, P2; P5, P6) 




e v L 


l -Pz(s 56 ) C v ( A ,;A,;Ai)(-- ■) 


»=1 


sin( 2 $vv) 


Pz{s 5e ) 



(170) 


O 

o 

o 

o 


afi'yS 
(+; + ;±) 
ci.fi'yS 
(+:-;±) 
a/3 7 <5 
(-;+;±) 
a. (3^ 8 

>±) 


_ T a 0 t7<5 t 

— J q' 12+ J q 34+ 56± ) 

_ T&P r7<5 t 

— J q' 12+ J q 34- 56+ 7 

_ r“/3 r7<5 7 

— J q' 12- J q 34+ 56+ 7 

_ -rap j 

— ^q' 12- Jg34- 56+ ■ 


with the weak couplings 


U R 7 


v_ = v L , 


R ! 


U L ■ 


(171) 


(164) 


The same decomposition is used for the case of identi¬ 
cal flavors, qqqqlL For the channels, qqudiv and 
qq du vt. we have 


For identical quark flavors, qqqqii, the basis reduces to 


n a 0 7<5 
W++;±) 


u (+-;±) 


_ 1 r a(3 T-yS J 

~ ^ J gl2+ J q 34+ J i56± 7 

_ T&P 77<5 T 
— J q 12+ J q3A- ^ 56+ 7 


_ e 2 y ( t) 

C(A q ;WT)i- ■■)= . 2 d Pw(s 56 ) C^A, ■ .) ■ 

z sin uw 

( 172 ) 

The coefficients for q'q'qq V are given by 
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~ ,. I A XiSn (l+,2-,;3+ 4^5+ 6J) 

C x (+-, + ;+){Pl,P2-,P3,Pr,P5,P6) = J_ R /? + 1- 0 + A~ \+ e-\ I > 

' jy n xfi n \ l q '> Z q' I' 5 ? i 4 ? i°£ i ) J 


A A ^,fin(lJ,2r,;3-,4;5+ 6J) 

C®(+; —;+) (Pi? P2? P3i P4i Pbi P6J ( R „ n + 9 — . q — /j + . t^+ ft —' 


j£s,fin(i+,,2^;3-,4±;5+6r); ’ 
C x{+ . ± ._ ) (p 1 ,p 2 ;p 3 ,P4;P5,P6) = C x(+ , ± . +) (p 1 ,p 2 ;P3,P4-,P6,P5) , 
Cg iV (_. + . ± )(pi,p 2 ;P3,P4;P5,P6) = -^g,u(+;-i±)(P 2 >Pi;^4.P3;P5,P6) , 

Ca(-, + -,±)(PliP2\P3iP4\P5iPG) = (£+._. + ) (P2,Pl',P4,P3',P5,P6) , 
C q , v (---±){pi,p 2 - l P3,P4- l P5,P6) = -Cg, l ,(+;+ ; ±)(P2,Pi;P4,P3;P5,P6) , 
Ca(---±){pl,P2\P3,P4\P5,P&) = C a ^ + . + .+)(p 2 ,pi;p4,P3;P5,P6) , (173) 

and for qqqqV they are given in terms of the amplitudes A x + n and B x fi n for q'q'qqV by 

A - x M x , fin (l+ 2=-; 3+ 4g-;5/,6-)-£B*, fi „(l+ 4r ; 3+ 2^;5+ 6-) 

^(++; + )(pi,P 2 ;P 3 ,P 4 ;P 5 ,P 6 ) = Uo n+V%+W 6 -'i_ 4 « n+ 4 -.4 2 -- 5 + r-i 

\ JV - D 3?,finV 1 (j » 5 °q- 5 1 5 ) ^ 1 x,finV 1 g 5 5 °q 1 1 5 / 

( 4 —;+) (.Pi ? ^ 2 ? P3? P4? P 6 j 1 _ 1 _ d „ /1 + 9 —. o — /i + . c+ £? — \ ) 5 

\ JV 1 5 5 J/ 

C x (+±.-){pi,p 2 ',P3,P4\P5,P6) = C x ( +±:+ ){pi,p 2 - l P3,P4-,P6,P5) , 

Cq,v{ -;±)(Pl,P2;P3iP4;P5!P6) = — ^g,t)(++;±)(P2iPi;P4,P3;P5iP6) j 

C a( ■±)(P1,P2;P3,P4;P5,P6) = C a ( ++ . ±) {p2,Pi;P4,P3',P5,P6) • (174) 


The various relations for the coefficients with flipped he- 
licities follow from Eq. (113811 and charge conjugation in¬ 
variance. The tree-level helicity amplitudes are given in 
Appendix IC3al 


A gg qq V 

For ggqqll , the helicity operator basis consists of 12 
independent operators, 


soabap _ * foa fob rap j 

H ++ ( +; ±) — 2 W+^2+ J q3i+ J t56± , 

n aba0 _ 1 „a K b j&P j 
H++(- ; ±) ~ 2 W+^2+ J q 34- ^56+ > 

/")ab aft _ foa fob raft j 

H-|—(+;±) — °1+ c 2- Jq 34+ 56+ i 

/oab ap _ foa fob raP j 

H-|—(-; + ) — ° 1+ i5 2 - ^q 34— 3^56+ i 

Q a b °tfi 

/+ a b 3/3 _ 1 foa fob jap j 

11-(-; + ) — o ^1- ^2- Jq 34- 3^56+ ■ 


1 


0.0 cup _ * foa rob rap j 

-(+; + )“ 2 ^2- J q 34+ 56+ 


(175) 


For the corresponding partonic processes are 

gg ud iv and gg du vi , and the helicity operator basis re¬ 
duces to six independent operators, 


O' 


ab ot.fi 
++(W 


_ N = Ib^bw j: 


OcP 


1+ u 2+ J ud. 34- 56- > 


Jfi i 


r\aba.p _ Ra K?f> T“^ 7_ 

'-'+_ (W-) — D l+ °2- ^ud 34— J £u56- > 


soob ap _ - m to 

^-(r+-) — 9 °1- °2- J ud34- u tvb3- > 


T a P 


o 


ab ap 


= l B ( !\ B f A J'l 


Jti 
Jvt 56— 


++ (W+) 2 2+ °du 34— 

0 Qba L,T = Bf + ^_ jf 34 _ ^se¬ 


ct 


X +-(W+) 

ab otfi 
— (W+) 




(176) 


We use the color basis 

f abap = (( T a T b ) a p , (r b T“)^ , tr[T“T b ] S a p) , (177) 
and the amplitudes are color-decomposed as 
A x {gig 2 q 3 q4hh) 

= i E [T^ 1) T“^)] Q3 . 4 A x (a( 1 ), ( 7(2);3 g ,4 ? ;5,,6,-) 

CT€lS2 

+ i tr[T Q 1 T a2 ] ^ 3fi4 B x (l, 2; 3„ 4*; 5/, 6 *). (178) 

Writing the Wilson coefficients in the decomposition 
in Eqs. (113611 and (11371) . we have for the ggqqil channel 

Cx 1 x 2 (\ q -\ 1 ){P1,P2\P3 1 P4\P3 1 P6) 

= e 2 f [Q e Q q + v{y x Pz(s 56 )]Cq XlX2{Xq ; Xe) (...) 


£[^ 


v \, 


-Bz(s 56) l^uAi A2(A g ;A^) (• • ■) 
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u \i 


s\n{29 w ) 
with the weak couplings 


Pz(s 56 ) C aXl \ 2 (\ q -\ e ){■ • ■) f ) 


(179) 


U R ! 


vi_ = V 


L > 


v 9 +=v%, vl = v q L . (180) 


in Refs. fl87l \m for 
qqqq and in Refs. |l93i ll94l] for gggg. 

1. qqq'q' and qqqq 


Refs. fl89l4T93 for qqq’q’, 


In the case of distinct quark flavors, qq q’q’, the helicity 
basis consists of four independent operators, 


For the W T channels gg ud tv and gg du vt, we have 

_ e 2y(t) 

C Xl X 2 (W*)(- • ■ ) = ~ . 2a P w(s 56 ) C qXl x 2 (---)(- ■ •) ■ 
z sin uw 

(181) 

The coefficients for gg qq V are then given by 

C xXl x 2 { + - + ){pi,P2',P3,P4',Pb,P&) 

(A xM ( l Al ,2 A2 ;3+,4r ; 5+6r)\ 

= A a , ifin (2 A =,l^;3+4r;5+,6 £ -) , 

\^,fin(l Al , 2 A -; 3+, 4^; 5+, 6^)/ 

Ai A 2 (+;-) (Pi, P2! P3, P4; p.5, p 6 ) 

= C xXlX2 ( + . +) (p 1 ,p 2 ;p3,P4;P6,P5) ■ (182) 

The remaining Wilson coefficients are obtained by charge 
conjugation invariance as follows, 

Cq,v Xl X 2 (-;±)(PliP2',P3,P4',P5,P6) 

= VC qtVXlX2 ( + -±)(pi,ft2]PA,PZ]P5,Ps) , 
C a \ 1 \ 2 (--±)(Pl,P2-,P3,P4-,P5:P6) 

= -VC aXlX2{+ . ±) (p 1 ,p 2 -,P4,P3-,P5,P6) , 

- f 010 \ 

with 7= 1 0 0. (183) 

V° 0 V 


n a/3^S _ toiP t7<S 
'-'(+;+) ~ J q 12+ J q' 34+ > 

n a/3jS _ toi/3 t7<5 

~ J q 12+ J q' 34- S 

n a/3jS _ toi/3 -rjS 
^(-;+) ~ J q 12- J q' 34+ > 

0B = 4-J?34-- (184) 

For identical quark flavors, qq qq, the helicity basis only 
has three independent operators, 

n a0^s _ ]_ r ap t7 S 
(++) — 4 Jl2 + Js4 + ’ 

_ t&( 3 j'yS 

u (-\—) — J 12+ J 34- 5 

(185) 

Here we have not made the flavor label explicit, since 
both quark currents have the same flavor. In both cases 
we use the color basis 

T a ^~ S = 2 T f (5 aS , 6 a p 5 yi ) . (186) 

The QCD helicity amplitudes for qqq’q’ can be color- 
decomposed in the basis of Eq. (11861) as 


* 4 ( 171 ^ 2 ^ 3 ^ 4 ) — 2Tpi d aia4 5 a3Ct2 A(\q, 2q‘, 3 q > , Aqf (187) 


+ jy d aia2 5oi 3 a 4 B{tq, 2qJ 3qr', 4 ? /) 


The tree-level helicity amplitudes are given in Ap¬ 
pendix 1C 3 bl 

VII. pp -+ JETS 

In this section, we give the operator basis and matching 
relations for pp —> 2,3 jets. We consider only the QCD 
contributions, so that quarks only appear in same-flavor 
quark-antiquark pairs with the same chirality, and so are 
described by the currents Jy ±. The helicity amplitudes 
for each channel are given in Appendix [D] 


A. pp —f 2 Jets 

For pp —> 2 jets, the partonic channels qqq'q', qqqq, 
qqgg, and gggg contribute. We will discuss each in turn. 
The one-loop helicity amplitudes for all partonic channels 
were first calculated in Ref. (23|. The tree-level and one- 
loop results are given in Appendix ID II The two-loop 
amplitudes have also been calculated, and can be found 


where we have included a factor of 1 /N for convenience. 
The amplitude vanishes in the case that the quark and 
antiquark of the same flavor have the same helicity. This 
is equivalent to the fact that the operators of Eq. (11841) 
provide a complete basis of helicity operators. For iden¬ 
tical quark flavors, the QCD amplitudes can be written 
in terms of the amplitudes for the distinct flavor case as 


-4(<7 i<729394) = *4(91929374) - *4(91947372) • ( 188) 

The Wilson coefficients for qq q’q' are then given by 

A=r ( Afi n (l^, 2q ; 3^1, Aq,) 

p 4) = ( iBan( !+ i2 - ; l Ji l,, ) 

At ~ ~ ~ \ l *4fin(l^) 2q ; 3 ,, 4^1) 

= (^B, n( l+, 2 -;3-,4+ ) 

C(--+){pi,P 2 -,P 3 ,Pi) = <?(+;-)(P2,Pi;P4,P3) , 

C(-,-){pi,P 2 ]P 3 ,P 4 ) = C {+ . +) {p 2 ,Pl-,P 4 ,P 3 ) , (189) 

and for qqqq they are given in terms of the amplitudes 
Afi n and B Rn for qqq'q’ by 
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C(++)(Pl,P2',P3,P4) 


C(+-){pi,P2;P3,P'l) 


C(- )(Pl,P 2 ;P 3 ,P 4 ) 


( A fin ( 1+ ; 3+ 4g ) - ^£ fin (l+ 4*; 3+ 2f )\ 

V^B fin (l+ 2r ; 3+ 4^) - A fin ( 1+ 4* ; 3+ 2J)J ’ 

/ A fin (l+ 2^;3-,4+) \ 

V^B fin (l+ 2r;3- ) 4+) y /’ 

C(++)(P2,Pi;P4,P3) • 


(190) 


The relations for C(_ ; ±) and C7_) follow from charge 

conjugation invariance. The Wilson coefficient C(_|_) is 

equal to C( +; _), since the amplitude vanishes when the 
quark and antiquark of the same flavor have the same he- 
licity (both + or both —), so there is no exchange term. 
The subscript “fin” in Eqs. (11891) and (11901) denotes the 
IR-finite part of the helicity amplitudes as discussed in 
Sec. El see Eq. ©. Recall that the symmetry factors 
of 1/4 in Eq. (I185D already take care of the interchange 
of identical (anti)quarks, so there are no additional sym¬ 
metry factors needed for CV_|_ + ). Explicit expressions for 
all required partial amplitudes at tree level and one loop 
are given in Appendix ID 1 al 


2- ggqq 

For ggqq, the helicity basis has a total of six indepen¬ 
dent operators, 


O' 


ab cx.fi 
++(+) 


— 1 K a K° T 

— „ Oil Dn I d- 


1+ ^2+ °34+ ) 


s^ab a/3 na yjb raft 

—(+) — °1+ °2- J 34+ > 


O' 


aba/3 _ 1 na ryb tOc/3 

-(+) — 9°1- °2- J 34+ > 


o 


ab otfi 
++(-: 


— - R a n° r 

— ^ 4 - *-> 9 . 4 - ? 


1+ ^2+ ° 34 - > 


o 


ab otfi _ 


o 


+-(-) 
ab otfi 


= £>?, B b j: 


1+ u 2- J 34- > 


1 


_ Yia job jcx.p 

-(-) — 2°!- °2— J 34- * 


(191) 


Note that the use of a helicity basis has made it easy to 
count the number of required operators. [f| For the color 
structure, we use the basis 


faba 0 = ( (T a T 6 )a - , (T b T a ) aB , tr[T“T b ] 6 a? ) . (192) 


9 This should be contrasted with the more complicated basis given 
in Eq. ( 126 ) of Ref. [l 95 l] which is built from fields Xn* and Bn^ 
and standard Dirac structures. It can be reduced to a minimal 
basis using identities such as O2 = — Oi, Os = 07 + 4^03 — AtO^ 
and 06 = O5 — 20 i + O(e) where t = -71.3/2, and then 

can be related to the basis used here. 


The color decomposition of the QCD helicity ampli¬ 
tudes into partial amplitudes using the color basis of 
Eq. (11921) is 

A(gig 2 9394 ) 

= i E [T“^)T“^)] Q3 _ 4 A(a(l), ( 7(2);3 g ,4 ? -) 

<t£S2 

+ itr[T ai T° 2 ] 5 a3 ai B( 1,2; 3„ 4,-), (193) 

from which we can read off the Wilson coefficients, 

/A fin (l+,2-;3+,4,r)\ 

C + _ {+) {p l ,p2\P3,P4) = -4fln(2“,l + ;3+,4 g ) , 

\i? fi „(l + ,2-;3+,4r)y 

/A fin (l+,2+;3+4£)\ 

C ++{+) (Pi,P2;P3,P4) = A fln (2+,l+;3+,4f) , 

\B Rn (l+,2+-3+,^)J 
(A &a ( l-,2-;3+,4 9 -)\ 

c —(+)(pi,P2;P3,P4) = A fin (2-,l _ ;3+,4 5 ) . (194) 

VB fi „(l-,2-;3+,4r)y 

The remaining coefficients follow from charge conjugation 
as discussed in Sec. UTTeI 

C\ 1 \ 2 (-)(P1,P2',P3,P4) = VC Xi x 2 ( + ){Pi,P2;P4,P3) , 

0 -1 0 \ 

-1 0 0 . (195) 

0 0-1/ 

At tree level, the partial amplitudes are well known, and 
only the first two entries in C + _(-m are nonzero. Explicit 
expressions for all amplitudes at tree level and one loop 
are given in Appendix ID 1 1)1 



3■ gggg 

For gggg, the helicity basis has five independent oper¬ 
ators, 


r^abcd. _ 1 yya fob yoc yoc1 

'-'++++ — 41 °l + D 2+ i5 3+ i5 4+ > 

O a + Y + _ = ± i B a 1+ B b 2+ 13 c 3+ BL, 
Off__ = ^ + B b + B c 3 _Bt, 
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O abcd 1 yoa fob yoc yod 

-+ — 0 l- 0 2- 0 3-°4+ 5 

O ab l d __ = ^ B°_B b _B^_B d _ . 

We use the color basis 


T 


abed 


l 

2 • 2Tp 


/tr [abed] + tr [dc6a]\ 
tr \acdb] + tr [bdea] 
tr [adbc] + tr [cbda] 
2tr[a6]tr[cd] 
2tr[ac]tr[d6] 

\ 2tr[ad]tr[&c] / 


T 


where we have used the shorthand notation 


(196) 


(197) 


tr[o6] = tr[T a T b ], tr [abed] = tv[T a T b T c T d } . (198) 
Under charge conjugation, the operators transform as 

C seabed rpabed r 1 s^idcba rpabed /-i r\r\\ 

( ^AiA 2 A 3 A 4 i U - Cy AiA 2 A 3 A 4 i * 

Thus, charge conjugation invariance of QCD leads to 

seabed _ r'ideba fOC\C\\ 

°AiA 2 A 3 A 4 — u AiA 2 A 3 A 4 • V ZUU J 

In principle, there are three more color structures with a 
minus sign instead of the plus sign in the first three lines 
in Eq. (119711 . Since charge conjugation is a symmetry 
of QCD, Eq. (I200D holds to all orders, so these addi¬ 
tional color structures cannot contribute. In particular, 
the color structures in Eq. (11971) cannot mix into these 
additional structures at any order. Hence, it is sufficient 
to consider the reduced basis in Eq. (11971) instead of the 
9 different color structures, which were used for example 
in Ref. [H]. Note that for N = 3 it is possible to further 
reduce the color basis by one using the relation 

tr[a6cd + deba] + tr [aedfr + bdea ] + tr[ad6c + cbda] 

= tr[a&]tr[cd] + tr[ac]tr[d&] + tr[ad]tr[bc]. (201) 

We refrain from doing so, since it makes the structure of 
the anomalous dimension matrix less visible, and because 
there are no such relations for N > 3. 

The color decomposition of the QCD amplitude into 
partial amplitudes using the color basis in Eq. (I197D is 


Agigigsg*) = -^r 


t r [ a er(l) a o-(2) a o-(3) a <r(4) ] 

ct£<S 4 j 

x xT(cr(1), cr(2), cr(3), cr(4)) 

T ^ ' tr[o cr (i)a 0 .(2)]tr[a 0 .(3)a 0 .(4)] 

cres 4 /z| 


X R(cr(l),a-(2),cr(3),cr(4)) 


( 202 ) 


from which we obtain the Wilson coefficients 


C++ - {pl,P2,P3,Pi) 


/2H fin (l+,2+,3-,4-)\ 
2H fin (l+,3-,4-,2+) 
2H fin (l+,4-,2+,3-) 
B fin (l + ,2+ 3-,4") 
H fin (l+,3-,4-,2+) 
VR fi n(l + ,4-,2+,3-)/ 


C+++-(Pl,P2,P3,P4) = 


C++++(Pi,P2,P3,P4) = 


/2H fin (l+,2+,3+,4-)\ 

2H fin (l+,3+,4-,2+) 

2H fin (l+,4-,2+,3+) 

R fin (l+,2+,3+,4-) 

R fin (l+,3+,4-,2+) 

VR fln (l+,4-,2+,3+)/ 

/2A fin (l+,2+, 3+,4+)\ 
2H fin (l+,3+,4+,2+) 
2H fin (l+,4+,2+,3+) 
B fln (l+,2+3+4+) 
B fln (l+,3+,4+,2+) 

V Rfin(l + , 4+, 2 + , 3 + ) / 


c- +{pl,P2,P3,P4) = C+++_(pi,P2,P3,P4) 

C - (Pl,p2,P3,P4) = C ++++ (p 1 ,P2,P3,P4) 


(••}«•[..] ’ 

<■■>+*[••] ' 
(203) 


The last two coefficients follow from parity invariance. 
The factors of two in the first three entries of the coef¬ 
ficients come from combining the two color structures in 
the first three entries in Eq. (I197D using charge conjuga¬ 
tion invariance in Eq. (12001) . 

The tree-level amplitudes are well known. At tree level, 
only the A amplitudes with two positive and two nega¬ 
tive helicity gluons are nonzero. Because the A ampli¬ 
tudes correspond to a single-trace color structure, which 
possesses a cyclic symmetry, the corresponding partial 
amplitudes are invariant under the corresponding cyclic 
permutations of their arguments. Explicit expressions 
for the required amplitudes at tree level and one loop are 
given in Appendix ID 1 cl 


B. pp -+ 3 Jets 


The four partonic channels gqqq'q', gqqqq, gggqq, 
and ggggg contribute to pp —> 3 jets, which we discuss 
in turn. The one-loop partial amplitudes for the differ¬ 
ent partonic channels were calculated in Refs. I 22 I [24l [25| . 
Tree-level results for the helicity amplitudes for each par¬ 
tonic process are given in Appendix ID 21 


1■ gqqq'q' and, gqqqq 

For the case of distinct quark flavors, gqqq'q the he¬ 
licity basis consists of eight independent operators, 


to/3 t7<5 

q 23+ J q' 45+ ’ 

tS/3 t7<5 

'q 23+ J q' 45- > 

tS/3 t7<5 
q 23— J q’ 45+ ’ 

1 - 3/3 t 7<5 


(204) 
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For identical quark flavors, gqqqq , the basis reduces to 
six independent helicity operators, 


In both cases we use the color basis 


O 


a afi-yS 
±(++) 
a Oifi-yS 

±(H—) 
a a/3'7<5 
±(—) 


1 


naafi'yS _ 


— — n a T a P t 

— - *->1 4- ^93-1- I 


a/3 j'yS 


1± e '23+ ^45+ > 


2 Tf (r“j 6rf, t; b 5 ah T“p 5 lh t;- 5 6 aP ) . 


_ K a T a P T^ S 
— 0 1± J 2 3+ J 45- 1 


/-\a a 0 ~/S _ kjo toi /3 ryS 

u ^' i — - D l± J 23- J 4:5- ’ 


(205) 


(206) 

The QCD helicity amplitudes for g qq q'q' can be color 
decomposed into partial amplitudes in the color basis of 
Eq. (12061) as 


-4(9192939495) = 27 >i -^20:5^0:40:3 4 ( 1 ; 2 q , 3 q \ 4 q >, 5 q >) + T a ' &3 ^2 0:5 A(l; 4 ^, 5 r ; 2 „ 3 ? ) 

T j^T a l a3 S ai a 5 B(l; 2 q , 3^; 4 q /, 5g/) + —T a ' &5 6 a2 a 3 B( 1; 4 9 /, 5^; 2 9 ,3g) 


(207) 


where we have used the symmetry qq o q'q', and inserted the factors of 1/iV for later convenience. The amplitude 
vanishes when the quark and antiquark of the same flavor have the same helicity (both + or both —), in accordance 
with the fact that the operators of Eq. (12041) provide a complete basis of helicity operators. For identical quark flavors, 
the amplitudes can be written in terms of the amplitudes for the distinct flavor case as 


-4(0192539455) = -4(9192939495) - A(giq 2 q' 5 q' 4 q3 ) • 
The Wilson coefficients for gqqq'q' are then given by 


<?+(+;+) (Pi;P2,P3;P4,P5) = 


f A Rn ( l+;2+,3r;4+,5r,) \ 

A Rn ( 1 +;4+ 5 ^ 52 + 35 ) 
iB fin ( l+;2+,3r;4+,5y) 
ViB fi „(l + ;4+,5-,;2+,3r)y 

/ 4 fin (l+; 2+,3r;4-„ 5+) \ 
A fin (l+;4-„5+;2+3r) 
iB fin (l+;2+3r;4-„5+) 
\^i? fi „(l + ;4- / ,5+;2+,3r)/ 


c+(+;-)(pi;p2,P3;p4,P5) = 

JY J -'nnv, J - , ^q ,^q , 

\jr B U l+;4-„5+; 

and for gqqqq they are given in terms of the amplitudes ^4fi n and 13fi n for gqqq'q' by 
C+( ++ )(pi;P2,P3;P4,P5) = 


Mfi n (l + ;2+,3 5 ;4+,5 ? ) - ^S fin (l+;2+,5^ ;4+,3 5 )\ 

4fi n (l + ; 4+ 5=-; 2+ 3f) - ^i?fin(l + ; 4+ 3, ; 2+ 5*) 

1 


c '+(+-)(9i;p2,p3;p4,p5) = 


5 > “g ) 1 “g 1 ’ ’9 ’ ' J <2 > “q ’ "9 

iS fin (l+; 2+ 3=-; 4+ 5f) - 4 fln (l+; 2+ ; 4+ 3f) 

ViB fi „(l + ;4+5r;2+3 5 -)-A fin (l+;4+ 3r;2+5^)7 

/ A fin (l+;2+,3 9 -;4-,5+) \ 

A fin (l+;4-,54;2+,3=-) 
iB fin (l+;2+ 3=-;4",5t) 
V^S fi „(l + ;4-,54;2+3r)y 


(208) 


(209) 


( 210 ) 


Charge conjugation invariance of QCD relates the Wilson 
coefficients, 


^A(-;±)(Pi;P2,P3;P4,P5) = VC A(+;=F ) (pi ; p 3 , p 2 ] p5, Pi) , 

C\( - )(Pl\p2,p3',p4,p5) = VC K++) (pi\P3,P2-,P3,Pa) , 

( 211 ) 


with 


V = 


( 0 -1 0 0 
-10 0 0 
0 0-10 
Vo 0 0-1; 


( 212 ) 


The remaining Wilson coefficients for a negative helicity 
gluon follow from parity invariance, 

C'-(+;±)(Pi;P2,p 3 ;P4,P5) 
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= C , +C-ST)(Pi;P2,P3;P4,P5) 

c '-(++)(pi;p2,P3;p4,P5) 

= C+( - )(pl-,P2,P3\P4,P5) 


(■■>«•[•■] 


<••>«•[■■] 


(213) 


Explicit expressions for all required partial amplitudes at 
tree level are given in Appendix ID 2 al 


2 . ggg gg 

For 33333 , we have a basis of eight independent helic- 
ity operators, 

ciabca.fi _ f rya fib fic T&fi 

Ly +++(±) — 31 °1+ D 2+ °3+ J 45± J 

ciabcafi _ f fia vyb tic fafi 

°++-(±) — 2 ° 1+ ° 2 + ° 3 ~ J 45± ’ 

O abcafi _ f fia fib tic ja.fi 

-+ (±) — 2 °1- °2- °3+ J 45± i 

0'^% } = B b _ B c 3 _ J 4 f ± , (214) 

and we use the color basis 


from which we can read off the Wilson coefficients, 


C++=F(+) (pi > • • • i Pi > Ps) — 


Mfin( 1+ 2+ 3=F;4+ 5f)\ 
A fin (2+, 3*,l+;4+ 5r) 
A fin (3^,l+,2+;4+52-) 
A fin (3T,2+,l+;4+5 ? -) 
A fin (l+, 3*,2+; 4+5*) 
A fin (2+,l+,3+;4+ 5,-) . 

B fin (3+, l+,2+;4+ 5*) 
B fin (l+,2+,3+;4+ 5*) 
B fin (2+,3T,l+;4+5 ? -) 
C' fin (l+,2+,3+;4+ 5*) 
\C' fin (3+2+,l+;4+ 5*)/ 

(217) 


Charge conjugation invariance of QCD relates the coeffi¬ 
cients with opposite quark helicities, 


Cx 1 x 2 \ 3 (-)(P1,P2,P 3 ;P4,P5) = VC Xl \ 2 \ 3 (+){Pi,P2,P3;P5,P4), 

\ 


/03x3 13x3 
13x3 03x3 


with V = 


43x3 


0 1 
1 0 J 


(218) 


where l raxn denotes the n-dimensional identity matrix 
and the empty entries are all zero. The remaining coeffi¬ 
cients follow from parity invariance 


rj-tabc ot(3 


( [T a T b T c ] a p \ 
[T b T c T a } a p 
[T c T a T b ] aB 
[T c T b T a } a0 
[T a T c T b } aB 
[T b T a T c } a ~ p 
tr[T c T a ]T b p 
tr [T a T b ]T^ 
tr[T b T c ]T“, 
tr[T a T b T c ] a 6 a p 

\tr [T c T b T a }S a ~J 


(215) 


The color decomposition of the QCD helicity ampli¬ 
tudes into partial amplitudes using Eq. (12151) is 


C -+ (±)(Pl,P2,P3;P4,P5) 

= C ++ -( T )(P1,P2,P 3 ;P4,P5) 

c - (±){Pl,P2,P3-,P4,P5) 

= C +++{T) (p 1 ,p 2 ,P3;P4,P5) 


(•■>+>[■•] 




(219) 


At tree level, the partial amplitudes are well known, and 
only the A amplitudes are nonzero. Furthermore, the 
partial amplitudes with all negative or all positive he¬ 
licity gluons vanish. Combining the charge and par¬ 
ity relations of Eqs. (12181) and (12191) . there are only 
three independent amplitudes at tree level, which we 
take to be A(l+ 2 +, 3"; 4+ 5 ? - ), A{ 2 +, 3", l+;4+, 5 ? -), 
and A(3 _ , 1 + , 2 + ; 4+, 5^ ). These amplitudes are given 
in Appendix ID 2 bl 


-4(313233 34®) = 1 E [T a ^T a ^T >-(»] ^ 

<tG S3 

x ^4( cr (l)) cr ( 2 ), u(3); 4 g ,5*) 

+ i tr[T a ^T a ^]T^ 

/ Z 2 

x%(l), ff (2), ff (3);4,,5 5 -) 

+ i ^2 tr 

crG S3/Z3 

xC(a(l), 4 ( 2 ),a( 3 ); 4 „ 5 ,-), (216) 


3■ ggggg 


For ggggg , the basis consists of six independent helicity 
operators, 


O 


abode 


o 


abode 

++++- 


o 


abode 
+++- 


= ^B1 + Bl + Bl + Bi + Bl +) 

= ^B a 1+ B b + B c 3+ B d 4+ Bt_ ) 

= 2^3 ]B a 1+ B b 2+ B c 3+ BtB e 5 _, 
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O a _ b f_ e ++ = ^ Bl_B b 2 _Bl_Bt + Bl + , 

0 “^-+ = \Bl_B\_Bl_Bi_Bl+ , 

0 ahcde__ = _L Bl_B\_Bl_B\_Bl_ . (220) 


^ ' tl [flcr(l) ^cr(2) ^tr(3) ] U[^cr(4) ^(t(5) ] 

erCSs / (Z3 XZ2) 


x ^(c r (l),cr(2),cr(3),cr(4),CT(5)) 

from which we obtain the Wilson coefficients 


(225) 


As before, we only need one operator for each number of 
positive and negative helicities. We use the color basis 


V 


abcde 


1 


2 • 2T f 


( tr [abode] 
tr [acdeb] 
tr [acbed] 
tr [abced] 
tr [abdec] 
tr [acbde] 
tr [adceb] 
tr [adcbe] 
tr [aebdc] 
tr [abdce] 
tr [aecbd] 
tr [acebd] 
(tr[ced] — 
(tr[a6e] — 
(tr [acd] — 
(tr[6ec] — 
(tr [adb] — 
(tr [ace] — 
(tr[6dc] — 
(tr[aed] — 
(tr[ac5] — 
\(tr[6ed] — 


— tr [edcba] \ 

— tr [bedca] 

— tr[de6ca] 

— tr [decba] 

— tr[ced6a] 

— tr[ed6ca] 

— tr \becda] 

— tr [ebcda] 

— tr[cd6ea] 

— tr [ecd&a] 

— tr [dbcea] 

— tr [dbeca] 
tr[dec])tr[a6] 
tr[e5a])tr[cd] 
tr[dca])tr[6e] 
tr[ce6])tr[ad] 
tr[6da])tr[ce] 
tr[eca])tr[6d] 
tr[cd5])tr[ae] 
tr[dea])tr[6c] 
tr[6ca])tr[de] 
tr[de6])tr[ac] / / 


( 221 ) 


where we have used the shorthand notation 


tr [a6 ■■■cd\= tr [T a T b ■ ■ ■ T c T d ] . 


( 222 ) 


A priori, there are twice as many color structures as in 
Eq. (12211) with a relative plus sign instead of a minus sign 
between the two traces. Under charge conjugation, the 
operators transform as 

r-i r\abcde rpabcde n 1 r^edcba (rfiabcde\ /'ooqA 

^ A 2 A 3 A 4 A 5 1 ^ — — 1 u 'AiA 2 A 3 A 4 A 5 

Therefore, charge conjugation invariance implies for the 
Wilson coefficients 


iabcde _ _ fiedcba 

AiA 2 AsA4A5 AiA 2 AsA4A5 ? 


(224) 


and hence these additional color structures cannot ap¬ 
pear at any order in perturbation theory, either through 
matching or renormalization group evolution. 

The color decomposition of the QCD amplitude into 
partial amplitudes using the color basis of Eq. (122 ID is 

-4(5i52ff35455) 


1 


^ ' tl [u^i) a<j(2 ) ^<t( 3 )^<t( 4) ^er(5)] 
GS 5 /Z 5 

X A(cr(l),cr(2),cr(3),cr(4),cr(5)) 


C+-H- {pi, ■ ■ • ,Pd) — 2 


/A fin (l+,2+,3+ 4-,5")\ 
A fin (l+,3+,4-,5-,2+) 
A fin (l+,3+,2+,5-,4-) 
A fin (l+,2+,3+,5-,4-) 
A fin (l+,2+,4-,5-,3+) 
A fin (l+,3+,2+,4-,5-) 
A fin (l+,4-,3+,5-,2+) 
A fin (l+,4-,3+,2+,5-) 
A fin (l+,5-,2+,4-,3+) 
A fin (l+,2+,4-,3+,5-) 
A fin (l+,5-,3+,2+,4-) 
A fin (l+,3+,5-,2+,4-) 
B fin (3+,5-,4-,l+,2+) 
B fin (l+,2+,5-,3+,4-) 
B fin (l+,3+,4-,2+,5-) 
B Rd ( 2 + ,5" 3+1+4") 
S fin (l+,4+,2-,3+,5-) 
S fin (l+,3+,5-,2+,4-) 
B Rn ( 2+ 4-,3+,l+,5-) 
S fin (l+,5-,4-,2+,3+) 
S fin (l+,3+,2+,4-,5-) 
VB fi „(2+,5-,4-,l+,3+)/ 

C -±(pi, ...,p 5 ) = C+++-kf(pi,... ,p 5 ) 


c -++(pi, ■ ■ • ,p 5 ) = C+++ -(Pl> • ■ • ,P5) 


<••>++[..] 

<->«•[••] ■ 
(226) 


For brevity, we have not written out the coefficients 

C+++_|_and C+++++. They have exactly the same 

structure as C+-H_with the replacements 4“ —► 4+ 

and 4 _ ,5~ —»■ 4 + ,5 + , respectively, in the arguments of 
the helicity amplitudes. The remaining Wilson coeffi¬ 
cients are given by parity invariance as shown. The over¬ 
all factor of two comes from combining the two color 
structures in Eq. (12211) . which are related by charge con¬ 
jugation. 

At tree level, all the B amplitudes vanish, as do all 
the amplitudes in C++++± and C _=p. By the par¬ 

ity relations given in Eq. (12261) . only the A amplitudes 
in C+-h _are then required for the tree-level match¬ 

ing. Since these amplitudes correspond to single trace 
color structures, which posses a cyclic symmetry, the re¬ 
quired partial amplitudes are invariant under the corre¬ 
sponding cyclic permutations of their arguments. There¬ 
fore, at tree level, there are only two independent am¬ 
plitudes, which we take to be Afi n (l + , 2 + , 3 + , 4 _ , 5“) 
and Ag n (l + , 2 + , 4 _ , 3 + , 5 _ ). These are given in Ap¬ 
pendix ID 2 cl Simplifications also occur at one loop, since 
the B amplitudes can be expressed in terms of sums of 
permutations of the A amplitudes Si- 
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VIII. RENORMALIZATION GROUP 
EVOLUTION 

In this section, we discuss the renormalization group 
evolution (RGE) of the Wilson coefficients. We start with 
a general discussion and give the solution of the RGE to 
all orders in perturbation theory. For completeness, we 
also explicitly derive the (known) anomalous dimension 
at one loop. To discuss the RGE, it is convenient to 
consider the operators O t in Eq. (El), which are vec¬ 
tors in color space. Lastly, we give explicit results, in 
a manifestly crossing symmetric form, for the relevant 
color mixing matrices for the color bases we have used 
in the previous sections. Since the operators’ renormal¬ 
ization is independent of their hclicity structure, we drop 
all hclicity labels throughout this section for notational 
simplicity. 


A. General Discussion 

The renormalization of the hard scattering in SCET 
can either be carried out as operator renormalization, 
where the relation between bare and renormalized ma¬ 
trix elements is (Ot)bare = z~ nq/2 Z^ ng/2 (O^Zo, or 

with coefficient renormalization where (Cd} bare G bare _ 
Z” ? / 2 Z n ^^ 2 (Ot)b are ZcC Ten . The relationship between 

the two is Zc = Zq 1 . Here Z^ and Z A are the wave- 
function renormalizations of the SCET collinear quark 
and gluon fields and A n , defined in Sec. IIIBI and 


Next consider the renormalization group equations, 
working to all orders in a s . As usual, the fj, indepen¬ 
dence of the bare operator implies the renormalization 
group equation for the Wilson coefficient 




d d{fi) 

d/i 


= 7c(/r) C(n), 


(230) 


where the anomalous dimension matrix is defined as 


7cO) 




- d 

. d In /r 


ZM 


(231) 


The solution of the RGE in Eq. (I230D can be written as 


C(ri = U( f x 0 ^)C(no), (232) 


with the evolution matrix 

■ An /i 


U(n o,/x) = V exp 


/In 


din p! jc(v') 


(233) 


Here, V denotes path ordering along increasing p,, and 
H > fiQ. The path ordering is necessary since 7 c{n) is a 
matrix in color space. 

The anomalous dimension matrix has the general form 

7 c(A = TcuspM//)] MA) + 7[“«(m)] . (234) 

where T cusp is the cusp anomalous dimension and A(p, 2 ) 
is a process-dependent mixing matrix in color space, 
which does not depend on a s . Its /i dependence is given 
by 


n 9 = n t + n 9 > n q =n+ + n q (227) 

are the total number of quark and gluon helicity fields in 
the operator (recall that there are two quark fields in each 
of the fermionic helicity currents). The UV divergences 
for (Ot^bare are given in terms of a local product (as 
opposed to a convolution over label momenta), since we 
are working at leading power where the operators contain 
a single field per collinear sector. 

Let us consider more explicitly how the renormaliza¬ 
tion works at one loop. The counterterm Feynman rule 
at this order is 

(O t ) tree ( Zl q/2 Z\l g/2 Z c - l) • (228) 

At one loop, the UV divergences of (CT) bare are p ro _ 
portional to the tree-level matrix element as (Ch) tree D , 
where D is a matrix in color space, which denotes the 
1/e 2 and 1/e UV divergences (with /r defined in the MS 
scheme) of the bare matrix element. The counterterm 
has to cancel these UV divergences so 


A (fi 2 ) = 1 (rigC A + n q C F ) ln (^) + MA) , (235) 

which will be demonstrated explicitly at one loop in 
Sec. lVIIIBl We can then perform the integral in Eq. (12331) 
by using the running of the coupling, da s (/r)/dln/r = 
(3(a s ), to switch variables from ln/x to a s . We find 


t/(/r 0 ,//) = 


(236) 


x Va s exp 


»7r(Mo,M) A(//q) + # 7 (mo.m) 


where V as now denotes path ordering along decreasing 
a s , with a s {n) < a s (n 0 ), and 




Kv{lio,A = 

/ da. 


J a s ( F o) 



77T (Mo, At) = 

/ da, 


Ja B {ii 0 ) 


roie(n) 


/ da. 


J a s (iio) 


bcusp (c^s ) 

P{0L S ) 

Tcusp (^s ) 

P(0is) 

7 («s) 

/3(a s ) ’ 


da( 


s(Mo) 


1 

jsR) ’ 


(237) 


(6^ ee (z; q/2 Z n / 12 Z C - l) = -(O t ) tree TI I (229) 


Up to two loops, the noncusp piece j(a^] in Eq. (12341) is 
proportional to the identity operator jl96l ll97| 


which fixes Zc at one loop. 


j(a s ) = {n q Ac+n g ~i 9 c ) 1 . 


(238) 
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FIG. 1. (a) Collinear one-loop diagrams, (b) Soft one-loop 
diagrams connecting two fields i and j in the operator. 


In this case, the evolution factor simplifies to 


f7(/U 0 ,/t) = e -( n a G A+nqC F )Kr(lJ.o,v)+K~,(iJ.o,v) 


x exp 


Vr(fio,li) A(n 2 0 ) 


(239) 


where Sij = 2 pi ■ pj. 

First, we recall the wave function renormalization con¬ 
stants. In Feynman gauge at one loop, 


Z £ = 1 -S > + •")■ 

QJ 1 

Za = 1 + --03o — 2 Ca + ■■■): 

47t e 


(241) 


where /? 0 = 11/3 CU — 4/3 TpUf is the one-loop beta 
function coefficient [see Eq. dEU)], and rif is the number 
of considered quark flavors. Here and below, the ellipses 
denote possible UV-finite terms, which are irrelevant for 
our discussion here. (Using the on-shell scheme for wave 
function renormalization, the Z, L contain UV-finite pieces, 
see Appendix [Gi l 


The collinear diagrams in Fig. 1(a) contribute 


i q c = n = 

n = 


ol s Cf f 2 

4-7T Ve 
cxsCa / 2 
4-7T Ve 2 


(2 2 2 

(— H- 

Ve 2 e e 


2 2 r 

- Li± + 1 

e e 

1 2 r 

- Li _l + ■ 

e e 


•)(O t ) t 


(242) 


Starting at three loops the noncusp anomalous dimension 
is not color diagonal, and starts to depend on a conformal 
cross ratio built from factors of Pi ■ p, fujl. (For earlier 
work beyond two loops see Refs. |199[ - |207| . The result 
of Ref. jl98j | implies that the conjectured all-order dipole 
color structure in Refs. |202i . l203j ] is violated.) 

The evolution factors Kr(fio, l*), and r/r(lio, 2) are uni¬ 
versal. Explicit expressions for the integrals in Eq. (12371) 
to NNLL order, together with the required coefficients 
for r cusp and the (3 function to three loops, are given for 
reference in Appendix [E] 


B. One-loop Anomalous Dimension 


The anomalous dimension 7 c(li) is process dependent. 
In this subsection, we derive its general form at one loop. 
The anomalous dimension of the operators is determined 
from the UV divergences in the effective theory. The rel¬ 
evant one-loop diagrams in SCET are shown in Fig. [1] 
In pure dimensional regularization the UV and IR di¬ 
vergences cancel such that the bare results for the loop 
diagrams vanish. To extract the UV divergences, we reg¬ 
ulate the IR divergences by taking the external particles 
off shell with pf = p? ± ^ 0 . 

Since all fields in the operators correspond to dis¬ 
tinct collinear directions, the collinear loop diagrams in 
Fig. 1(a) only involve one external line at a time. Differ¬ 


ent external lines can only interact through the exchange 
of a soft gluon, shown by the diagrams in Fig. |l(b)| 
When expressing our results, we use the notation [see 
Eq. 473D] 


L = ln (-yf) 


l »= in (~(/~ 10 


(240) 


where I l c denotes the result of the diagram for an external 
leg of type i, either quark or gluon . 

The soft diagrams in Fig. |l(b)| differ from each other 
only in their color structure. The result of the diagram 
connecting particles i and j (with i ^ j) is given by 


IV = 


4-7T 


2 2 2 

— Lii Li .L- La l + ■ 

e e e 




tree 2 c 2 c 
l j : 


(243) 


where t £ and are matrices in color space. From 
Eqs. (I242D and (I243D we see explicitly that the opera¬ 
tors only mix with respect to the color structure, with 
no mixing between operators with distinct helicities. 

The action of the matrix on the color space is to 
insert a generator acting on the color index of the ith 
particle, i.e., 


rrr;)-'‘- 

(;r 

(t t?)"‘ ai "' 


- T , T'-L- 

OLiPi 


_ _ rp-’-pi-- rp 


PiOii ? 


= i/' 


aicbi 


rp---b 


(244) 


for quarks, antiquarks, and gluons, respectively. Our is 
identical to what is usually denoted as Tj in the notation 
of Refs. f208l.[20i. 

To give an explicit example, consider gg qq. Then, for 
quark i = 3 and antiquark j =4 we have 

Qf £c £c _ Qa 1 a 2 a 3 a i (rj, £c ^01020304 

= T^ 3 , (245) 


while for gluon i = 1 and quark j = 3, 

q] fc fc _ Qa 1 a 2 a 3 a 4 ^j,a 1 cb 1 j,c _ rjib 4 a 2 ft 3 a 4 

13 J ± a 3 p 3 ± 


(246) 


Plugging in the explicit basis in Eq. (11921) and using the 
relations in Appendix IA 21 we can rewrite the resulting 
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color structures above in terms of the basis in Eq. (11921) , 
which yields 


(C F -\C a 0 0 \ 

ili% = - 0 C F -\C A 0 , 

\ Tp T f C f ) 

(\C A 0 T f \ 

i c J § = - 0 0 -7> I . (247) 

\ 0 -T f 0 ) 

The other combinations are computed analogously. 

In general, one can easily see that for i = j 

T 0l '" a " t? t? = Ci f ai '" a " , (248) 

where Cj = Cf for quarks and Ci = C A for gluons. 
By construction, the color basis T ai '"“ n conserves color, 
because each index corresponds to an external particle. 
Since measures the color charge of the itli particle, 
color conservation implies 


n 

fo! ' 01 " = 0 . (249) 

i =1 

As a simple example, consider gqq for which T al02 “ 3 = 
Tala 3 - In this case, Eq. (12491) gives 

• raicbi rpbi _j_ rpC _ rpa\ rpa\ rpC 

1 J 1 02fi 3 + i a 2 fe J ftfi3 1 a 2 P3 ± 

= (if a ^T bl + [T c , T ai ]) a2Si3 = 0 . (250) 

The total bare one-loop matrix element is given 
by summing Eq. (12421) for each external particle and 
Eq. (I243D for each pair of distinct particles. The infrared 
logarithms Li± have to drop out in the sum of all UV- 
divergent contributions. To see that this is indeed the 
case, we can use Eq. (12481) to rewrite the collinear contri¬ 
butions. Then, the sum of all Li± terms is proportional 
to 


(d f ) tr ee Li jc H + ^(L l± + L j± ) % % 

i i<j 

= (o t ) tree (E L ^i %+E %) 

i i^tj 

= (o t ) tree (E L ^i) (E *i) = 0 - ( 251 ) 

* j 

where in the last step we used Eq. (I249D . For the same 
reason the 1/e 2 poles in the soft diagrams cancel against 
half of the 1/e 2 poles in the collinear diagrams. The 
remaining UV-divergent part of the matrix element is 
given by 


(dt)treej) = (Qt)t 


Os_ 

4-7T 


Mhc) 


+ n q C F (^ + ^) A( M 2 ) 

where the color mixing matrix is given by 

^(^ 2 ) = - E % L b' • 

i<j 


(252) 


(253) 


Combining this result with the identities in Eqs. (12481) 
and (12491) . we can easily check that the fi dependence of 


A(/i 2 ) is as in Eq. (I235D : 


A( M 2 ) - A(^o) = -2 E*?«i In(^) =E^ ln (^) 

i<3 ^ i ^ 

= l(n g C A + n q C F ) ln(^) . (254) 

We can now compute the anomalous dimension of the 
operators. From Eqs. (I229D and (12521) . we find at one 
loop 


Zq = 1 — D — 1 


L }(Z A - 1) + ^(Z £ -1) 


(255) 


which using Eq. (|231|) yields the one-loop anomalous di¬ 
mension 


7c(m) = [ 4 ^(m 2 ) ^ Hn g /3o + n q 3 C F )]. (256) 


The coefficient of 4 in front of A(// 2 ) is the one-loop cusp 
anomalous dimension coefficient [see Eq. (|E2I) ]. The re¬ 
maining terms determine the noncusp 7 (a s ) in Eq. (12341) 
at one loop, 


7( a s) = — ( n g/3o + n q 3 Cf) 1 • (257) 

47r 


C. Mixing Matrices 

In this section, we give explicit expressions for the mix¬ 
ing matrices for the color bases used in Secs. El IVI1 and 
ivm For simplicity, we only give explicit expressions for 
up to four partons, but allow for additional colorless par¬ 
ticles, such as a Higgs or vector boson. The matrices are 
straightforward to evaluate using the color relations in 
Appendix IA 21 but become rather lengthy for more than 
four partons, due to the large number of allowed color 
structures, and are more easily evaluated in an auto¬ 
mated way (see for example Ref. mi). For convenience, 
we introduce the following shorthand notation for sums 
and differences of logarithms Lij , 

Lij-ki-... = Lij 4 - Lki + ... , 

Lij..../(kl....) = ( Lij ....) {Lki. ), 

with L.^ = ln(— Sij/n 2 — iO) as defined in Eq. (1731) . 


(258) 
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1. Pure Gluon Mixing Matrices 


For gg and ggg in the bases used in Eq. (11051) and Eqs. (II101) and (11571) . we have 


^■ggiP 2 ) ~ Ca L\2 , A ggg ([i 2 ) — -Ca Fi2 i3-23 


1 0 
0 1 


(259) 


For gggg in the basis used in Eqs. (11301) and (11971) . we have 


A ggggil i ) 


( 5^4^12-14-23-34 0 0 

0 5CA-b12.13-24.34 0 

0 0 5CA-b13.14.23.24 

TfLi 2.34/(13-24) TfL\2- 34/(14-23) 0 CaL\2-3A 

0 TfLi 3 . 24/(14-23) 2>b 13 .24/(12-34) 0 

\TfL 1 4.23/(13-24) 0 7 1 pb 14 .23/(12-34) 


2TpL 14 . 23/(13.24) 0 

2TpL 13 . 24/(14-23) 2TpL 12 . 34/(14-23) 


2 TfL 12 . 34 /(i 3 . 24 )\ 

0 


0 


0 


2TpLi4. 23/(12-34) 2Tf-L 13 . 24/(12-34) 


0 

CAbi 3 .24 

0 


CaL±4 •23 / 


(260) 


For our color bases formed from multitrace color structures, the structure of the mixing matrices is simple. Since the 
mixing matrices are determined by single gluon exchange, cyclicity is maintained, and all that can occur in the mixing 
is that a single trace splits into two or two traces recombine into one. For example, the color structure tr[T a T b T c T d ] 
can only mix with 

tr [T a T b T c T d ], tr[T a T b ] tr[T c T d ], and tr[T d T a ] tr[T b T c ]. (261) 


Therefore, although the mixing matrices quickly get large as the number of color structures grows, their structure 
remains relatively simple. (An alternati ve a pproach to the organization of the anomalous dimensions for a large 
number of partons has been given in Ref. m .) For the dijet case, i.e., in the absence of additional colorless particles, 
the kinematics simplifies to 


S = Sl2 = S34 j t = Si 3 = S24 , U = S\4 = S23 , 
and these matrices were given in Ref. Il53| . which also gives their eigenvectors and eigenvalues. 


(262) 


2. Mixing Matrices Involving qq Pairs 

For qq and gqq in the bases used in Eq. (11411) and Eqs. (11 101) and (11491) . we have 

^qq(p 2 ) = Cf bl2 , A gqg(tl 2 ) = — CaL\2-13 + (2Cj? — Ca)L 2 3 
For qqq'q' in the basis used in Eqs. (11181) . (11861) . and (I167D . we have 
A - (,, 2 \ - A t- (Cf LlA-23 + (Cf ~ \Ca) b 12 -34/(13-24) 

^qqqqKP ) ~ L ±qq q'q' \P ) ~ I rp j 

\ Iff 12-34/(13-24) 


TpL 

14-23/(13-24) 

C F L 

12-34 + {C f -\C a )L 

14-23/(13-24) 


For ggqq in the basis used in Eqs. (11241) . (11921) . and (11771) . we have 


(263) 


(264) 


A -F- 2 
gg qq 


J C'A bl2-13-24 + {Cf — 5 Ca) L 3 4 
(P 2 ) =[ 0 

T F L 

12-34/(14-23) 


0 T F L 

13-24/(14-23) 

\Ca bi2-i4-23 + (Cf — \Ca) b 34 TpL 14-23/(13-24) 
Tf bi2-34/(13-24) Ca bi2 + Cf b 34 y 


(265) 


Again, these simplify in the dijet case, for which they 
were given along with their eigenvectors and eigenvalues 
in Ref. fUl. 


D. Soft Function Evolution 

In this section, we review the renormalization group 
evolution of the soft function, focusing on our use of the 
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color basis notation of Sec. lIII Dl for nonorthogonal bases. 
We will consider the particular case of the IV-jettiness 
event shape [95], which allows for a definition of exclu¬ 
sive iV-jet production with a factorization theorem of the 
form of Eq. ©. 

The color mixing matrices of the previous section are 
in general complex valued for physical kinematics. For a 
physical channel, some of the appearing Sij are positive, 
giving rise to imaginary terms from the logarithms, as in 
Eq. (l73l) . Since the cross section is real, these imaginary 
terms generated by the renormalization group evolution 
must drop out of the final result. We start by describing 
the properties of the soft function that ensure that this 
is the case. 

Recall that the hard function H K for a particular par- 
tonic channel k has its color indices contracted with those 
of the soft function. Explicitly, 


tr (H K S K ) = 


(266) 


E [Wt 


Pn 

(-A„) 


Qbl ■■■finO'l •••Qin S^ a l'" a n 


{Ai} 


The soft function is defined as a vacuum matrix element 
of a product of soft Wilson lines Y as 


S K (M,{m}) 


T^(W) S(M-M) T Y({m}) 


0 

(267) 


where T({rq}) is a product of soft Wilson lines in the 
rii directions. It is a matrix in color space, and Y 1 is 
its Hermitian conjugate. Here T and T denote time or¬ 
dering and antitime ordering respectively. The matri¬ 
ces Y and Y 1 are multiplied with each other, i.e. one 
of the color indices of the corresponding Wilson lines 
are contracted, and the external indices correspond to 
b\ ■ ■ ■ 0 n and a± ■ ■ ■ a n , respectively. Thus, for example 
yty = 5 aibl ■ ■ ■ _ The dependence of the soft func¬ 

tion on the particular measurement, as well as the details 
of the jet algorithm, are encoded in the measurement 
function M, whose precise form is not relevant for the 
current discussion. 

From the definition of the soft function in Eq. (12671) 
we see that it is Hermitian, namely (S£ 1 '"^"“ 1 "'“ n )* = 
gai—j n abstract notation, this means Sj. = 
S K: which implies that the product C'S K C appearing in 
the cross section is real, so imaginary terms that appear 
in the Wilson coefficients due to renormalization group 
evolution drop out in the final cross section. 

While this argument is trivial in a basis independent 
form, it is important to emphasize that in a nonorthog¬ 
onal basis it takes a slightly more complicated form. As 
discussed in Sec. MIDI in a specific nonorthogonal color 
basis, Eq. (12661) takes the form & S K C = C* t T S k C as 
in Eq. (foSj) , where the matrix T is defined in Eq. (157>1) . 
Similarly, the matrix representation of S K is not Hermi¬ 
tian with respect to the naive conjugate transpose of its 


components. Instead, the condition on the reality of the 
cross section is given by [see Eq. m] 

S K = Sl = f~ 1 S* T f. (268) 

The invariance of the cross section under the RGE 

fi-^-a N =0, (269) 

d/r 

implies relations between the anomalous dimensions of 
the SCET functions appearing in the factorization the¬ 
orem of Eq. ([2[) . In particular, it allows the anomalous 
dimension of the soft function to be determined from the 
anomalous dimensions of the Wilson coefficients, along 
with the anomalous dimensions of the beam and jet func¬ 
tions. The anomalous dimensions of the jet and beam 
functions are proportional to the color-space identity. 
The anomalous dimensions of the beam and jet func¬ 
tions appearing in the 7V-jettiness factorization theorem 
are equal to all orders in perturbation theory t212i j allow¬ 
ing us to use only the jet function anomalous dimension 
in the following discussion. Renormalization group con¬ 
sistency then implies that the contributions of the soft 
function anomalous dimension not proportional to the 
identity, including the color off-diagonal components, are 
completely determined by the anomalous dimensions of 
the Wilson coefficients. 

The soft function for N -jettiness can be written in the 
general form of Eq. (12671) , but with an explicit measure¬ 
ment function 


Sk (^a ? > • ■ * j kN 5 ) 


(270) 


= 0 


Ty f (K}) I] Sfc-ti) TFdnJ) 


Here 71 picks out the contribution to the JV-jettiness ob¬ 
servable from the momentum region i, whose precise def¬ 
inition can be found in Ref. j98[. The soft function for 
7V-jettiness was first presented to NLO in Re f. [98| . and 
more recently analyzed to NNLO in Ref. [213| . 

The all-orders structure of the renormalization group 
evolution for the soft function can be derived from 
Eq. (12691) . and is given by [H, ;98|] 


V ■^S K ({k i },iJ,) 


U dk i 


2 L 


7s({fci - /c'},m)S'«({A:'},/z) 


■S K {{k'i},ix) 7 s({^i - K},n) 


(271) 


The soft anomalous dimension 75 , and its conjugate 7 ^, 
are given in terms of the anomalous dimension 7 j of the 
jet function and the anomalous dimension of the Wilson 
coefficients, 7 c defined in Eqs. (I230D and (I231D . as 

i 
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~ 2 lc(VYi 5 ( k i) ■ 


(272) 


(Here, the Qi are related to the precise IV-jettiness def¬ 
inition, see Ref. [98].) The Hermitian conjugates of 7 c 


and 7 s above again refer to the abstract Hermitian conju¬ 
gate in color space. In a nonorthogonal color basis, they 
are given in terms of the complex conjugate transpose 
components according to Eq. m as 


rji —1 '~'*T rji 

7 c= T 7c T , 


-'"'-t rji —1 rji 

Ts = t 7 s T - 


(273) 


IX. CONCLUSIONS 


carefully treated, and expressions for the color sum ma¬ 
trix T are given for the used color bases for all processes 
considered in the paper. Explicit results are also given 
for the one-loop mixing matrices describing the renor¬ 
malization group evolution in color space for the case of 
pp —> up to 2 jets with an arbitrary number of uncolored 
external particles and in a manifestly crossing symmetric 
form. 

Combining the methods of this paper with known ex¬ 
pressions for jet, beam, and soft functions for particu¬ 
lar exclusive jet cross sections, or jet shapes/observables, 
should facilitate analytic resunnnation for a large number 
of processes for which fixed-order amplitudes are known, 
or are soon to be calculated. 


In this paper, we have presented a helicity operator ap¬ 
proach to SCET. Helicities are naturally defined with re¬ 
spect to the external lightlike reference vectors specifying 
the jet directions in the effective theory, eliminating the 
need to consider complicated Lorentz and gamma matrix 
structures in the operator basis. The helicity operators 
correspond directly to physical states of definite helicity 
and color, which when combined with color organization 
techniques, greatly simplifies the construction of a min¬ 
imal operator basis. Furthermore, the helicity operators 
are automatically crossing symmetric, and make mani¬ 
fest parity and charge conjugation symmetries, making 
it simple to determine relations amongst Wilson coeffi¬ 
cients. 

We demonstrated the utility of the helicity operator 
approach by explicitly constructing the basis valid to 
all orders in perturbation theory for a number of key 
processes at the LHC involving jets, and then determin¬ 
ing the matching coefficients. In particular we consid¬ 
ered pp —> H + 0,1 jets, pp W/Z /y + 0,1 jets, and 
pp —> 2 jets at next-to-leading order, and pp —> H + 2 
jets, pp —► W/Z /7 + 2 jets, and pp —> 3 jets at leading or¬ 
der. We also discussed the dependence of this matching 
on the regularization scheme, considering schemes with 
helicities in 4 and d, dimensions. An important and well- 
known simplification of the SCET approach is that when 
dimensional regularization is used for both IR and UV 
divergences, all loop graphs in the effective theory are 
scaleless, and thus vanish. As a result, the hard SCET 
Wilson coefficients in the MS scheme, determined from 
matching QCD to SCET, are given directly by the IR- 
finite parts of color-ordered helicity amplitudes, defined 
using Eq. (l96l) . The use of our helicity operator basis 
therefore makes it simple to combine analytic resumma¬ 
tion in SCET with fixed-order calculations of helicity am¬ 
plitudes. 

The all-orders structure for the renormalization group 
evolution of the helicity operator basis was discussed in 
detail. At leading power, distinct helicity structures do 
not mix, with renormalization group evolution causing 
mixing only in color space. This feature is made man¬ 
ifest at the level of the SCET Lagrangian due to the 
expansion in the soft and collinear limits. Subtleties as¬ 
sociated with the use of nonorthogonal color bases were 
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Appendix A: Spinor and Color Identities 


1. Spinor Algebra 


The overall phase of the spinors | p±) is not determined 
by the Dirac equation, ( |p±) = 0, and so can be chosen 
freely. In the Dirac representation, 


(1 °\ 

i i 

( 0 crA 


(0 -l) ’ 

■ 7 = | 

0 ) ' 

, 75 


° M 

1 Oj ’ 


(Al) 

and t aking n/ = ( 1 , 0 , 0 , 1 ), we have the standard solu¬ 
tions |l46j 



( Vp~ \ 


/ vT+e ^ \ 

1 


■ |p -> = 7! 


71 


— sfpVe lcf>p 


Vv/p+e 1 ^/ 


\ J 


(A2) 

where 

P ± =/tp 3 > exp(±i</>p) = p . (A3) 

VP p 
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For negative p° and p± we use the usual branch of the 
square root, such that for p° > 0 

|(-p)±) =i|p±). (A4) 

The conjugate spinors, (p±|, are defined as 

(p±| = sgn(p°) \p±) . (A5) 

The additional minus sign for negative p° is included to 
use the same branch of the square root for both types of 
spinors, i.e., for p° > 0 

{(-P)±I = -IF^F) = -(-i)(p±l = i<P±| • (A6) 

In this way all spinor identities are automatically valid 
for both positive and negative momenta, which makes 
it easy to use crossing symmetry. The additional signs 
only appear in relations which involve explicit complex 
conjugation. The most relevant is 

<p-|9+>* = sgn(A°) (q+\P~) • (A7) 

The spinor products are denoted by 

( pq } = (p-\q+), [pq] = ( p+\q-) • (A8) 

Similarly, for products involving additional gamma ma¬ 
trices, we write 

(pNM = (p-|7 M l q~) , [pNk) = {p+\Y\q+) , (A9) 
(p\k\q\ = {p—\lf\q —), [p\k\q) = (p+|#k+) > (A10) 

ip\qk\l) = {p-\0\l+} , \p\qk\l] = (p+|#|Z-), (All) 

etc. 

Some useful identities, that follow directly from the 


definition of the spinors, are 

{pq) = ~{qp), [pq] = ~[qp\ , (A12) 

[pNb) = (pNb] = 2p M • (A13) 

From the completeness relations 

\p±){p±\ = ^*fi, (A14) 

i> = IpKpI + \p)\p\, (Ai5) 

one finds 

{pq) [qp\ = i tr{ (1 - 75 W} = 2 P ' 9 ■ (A16) 

Combining this with Eq. (lA7l) . it follows that 

\{pq )I = IN]I = V\ 2 P-q\ ■ ( Al7 ) 

The completeness relation is also useful to reduce typical 
expressions like 

\p\q\k) = \pq\(qk), (A18) 

to spinor products. 


Charge conjugation invariance of the current, the Fierz 
identity and the Schouten identity are 

(pNM = feNb). (Ai9) 

\p\ln\q)[k\l li \ l ) = 2 \pk](lq ), (A20) 

( pq)(kl) = ( pk)(ql) + (pl)(kq) . (A21) 

Finally, momentum conservation X]”=i Pi = 0 implies 

n 

Y,\ji]{ik)=0. (A22) 

i=1 

From Eq. (1A2I) , we see that under parity the spinors 
transform as 

|p p ±)=±A 7 0 |pt), (A23) 

and therefore 

{p p g p )=-e i (^+^«)|N], (A24) 

[p p < 7 P ] = —e~ 1 ^ p+tl>q \pq) . (A25) 

When applying the above result to a helicity ampli¬ 
tude, the phases which appear ar e det ermined by the 
little group scaling (see e.g. Refs. [l46l Il47l . 12141 ] for a 
review). The little group is the subgroup of the Lorentz 
transformations that fixes a particular momentum. In 
terms of the spinor helicity variables, the action of the 
little group, which preserves the momentum vector p, is 
given by 

Ip)-Mp), \p\^z\p\- ( A26 ) 

In the case that the particle with momentum p has he¬ 
licity h 7 the corresponding helicity amplitude scales as 
z~ 2h under the little group scaling. This property of the 
helicity amplitudes then predicts the phases that appear 
in the amplitude under a parity transformation. 

The following completeness relation for the polariza¬ 
tion vectors is also useful 

g) ( e v(p> g))* = - 9 ^ + Vil<lv ^ v a v% • ( A27 ) 

A=± P ^ 

In SCET the collinear quark fields produce projected 
spinors 



The projected spinor trivially satisfies the relation 

*(^|p±>)=0, (-A29) 

so it is proportional to |n±). Working in the basis in 
Eq. (|A2[) , we have 
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_l_ e *Op <M g j n 


fl P]=V? 


a i{it> p 4> n ) cog I 


sin 


9r. 


- sin 


e p 

sm -f 


e p 

cos -f 


*}> 


(A30) 


Here 9 n ,(f > n , and 9 P ,4 > P , are the polar and azimuthal angle 
of the n and p vectors, respectively. In particular, we see 
that choosing = P^ /p°, which can always be done at 
leading power since there is a single particle per collinear 
sector, we have <f> p = <t> n , 9 p = 9 n , and the simple relation 


where we have chosen the standard normalization 


T F = 

2 


(A36) 


Throughout the text, and for the amplitudes in the ap¬ 
pendices, we have kept Tp arbitrary. This can be done 
using C F = T F (N 2 — 1 )/N, C A = 2T F N. The strong 
coupling constant, g Sl can be kept convention indepen¬ 
dent, by using g s -> g s /y/2T F . 

Some additional useful color identities are 


: \P±) = 


n ■ p 




= (C r 


(A 3 7) 

n±) . 

(A31) 

T a T b T c T a = T 2 § b Cl + ^ Cf _ 

T b T c ; 

(A 3 8) 


2. Color Algebra 

The generators of a general irreducible representa¬ 
tion r of SU(N) satisfy 

[t a r ,t b r \=if abc t c r , t a r t a r =C r 1, tr [t a r t b r ]=T r 8 ab , 

(A32) 

where f abc are completely antisymmetric, and C r is the 
quadratic Casimir of the representation r. The normal¬ 
ization T r is given by T r = C r d r /d , where d r is the di¬ 
mension of the representation and d the dimension of the 
Lie algebra. 

We denote the generators in the fundamental repre¬ 
sentation by tp = T a , and the overall normalization is 
fixed by choosing a specific value for T F . The adjoint 
representation is given by {t a A )bc = —i f abc , which implies 

facdfbcd = Ca 6 ab ( A33 ) 

We also define the symmetric structure constants as 

d abc = _L tr [ T a{ T b, r c|j _ ( A34 ) 

T f 

For the fundamental and adjoint representations we have 
d F = N, d A = d = N 2 — 1, and so 

N 2 - 1 

C F = ^ r , C a =N , (A 3 5) 

I 


where the second relation is equivalent to the complete¬ 
ness relation 

Kp T yS = t f (&aS ~ Sjg'j . (A39) 

We also have 

T b i jbac T c = -A T a t (A40) 

rpc • j?cad^jdberj^e rp 2 ^ ab Tt A rparpb (A41 j 


3. QCD Color Decompositions 

Here we briefly review a common color decomposition 
for QCD NLO amplitudes @-{8j]. The color bases used for 
the processes discussed in the text are specific examples 
of the decompositions given below, and were chosen to 
facilitate the extraction of the matching coefficients from 
the amplitudes literature. For a pedagogical introduc¬ 
tion to color d ecom positions in QCD amplitudes see for 
example Refs. [l46lll47| . 

For an n gluon process, a one-loop color decomposition 
in terms of fundamental generators T a is given by 


A n ( gi ---g„)=g :- 2 Y trir-w-r-w] 

<res„/z„ 

Ln/2J+1 

+ 97 Y Y trp 1 ®^ 1 ) ■ ■ • T a ‘ T (°- 1 )]tr[T a ' T < 0 ) • ■ ■ T a<T ( n )]A„ ;c (cr(l), ■ • -, a(n)) , (A42) 

c—3 <j£S n /S c — i,n—c+i 

where A n;1 , A^ ree are primitive amplitudes, which can be efficiently calculated using unitarity methods, and the A n . c 
are partial amplitudes which can be written as sums of permutations of the primitive amplitudes. The amplitudes 
appearing in this decomposition are separately gauge invariant. In this formula, S n is the permutation group on n 


AT e (a(l), ■ ■ ;a(n)) + g 2 s C A A n -i (<r(l), • • ■, a(n)) 
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elements, and Sij is the subgroup of Si+j which leaves the given trace structure invariant. At tree level, only the 
single trace color structure appears. 

In the case that additional noncolored particles are also present, an identical decomposition exists, since the color 
structure is unaffected. For example, for a process involving n gluons and a Higgs particle, the amplitude satisfies the 
same decomposition as in Eq. (IA42I) . but with the partial and primitive amplitudes in Eq. (IA42I) simply replaced by 
A(</>, ct( 1), • • -, cr(n)), where (j> denotes the Higgs particle |l77f . 

A similar decomposition exists for processes involving qq pairs. For example, the one-loop decomposition for a 
process with a qq pair and n — 2 gluons is given by [25] 


A„(q iq2 g 3 ...g n ) = <?r 2 Y, ( T<M3) • • • T<Mn) U- 

crG*S' n _2 

n— 3 

+ 9sJ2 E tr[r°'W ...T a ^](T a ^ ••■r°'W) aj gA n;c ( lff , 2 , ; <T(3 ) ,---,tr(n)) 

C=3 aes„_ 2 /z c -i 

+ 9. n E tr P^ 3) • ■ l,-,2,;<7(3),. --^(n)) 

aGSn-2/ Z n — 3 

+ g n s Y, trpT^M ■ • • T a ' M ]t a p A n;n-i (Iff, 2,; <r(3), ■ • a(n)) . (A43) 

O’GStl — 2 / Z n — 2 


^n ee ( 1 g ) 2 g ;CT(3),- • -,CT(n)) + CUA n; i (l,, 2 9 ; cr(3), • • -,a{n)) 


This decomposition is easily extended to the case of additional qq pairs. As with the gluon case, the same color 
decomposition also applies if additional uncolored particles are included in the amplitude. 

For more than five particles, the one-loop color decompositions given above do not give a complete basis of color 
structures beyond one loop, since color structures with more than two traces can appear. A complete basis of color 
structures is required for the SCET basis to guarantee a consistent RGE. A convenient basis of color structures for 
one-loop matching is then given by extending the one-loop decomposition to involve all higher trace structures. 


Appendix B: Helicity Amplitudes for Higgs + Jets 

In this appendix we give explicit results for the hard matching coefficients for H + 0,1,2 jets. We only explicitly 
consider gluon-fusion processes, where the Higgs couples to two gluons through a top-quark loop, and additional 
jets correspond to additional gluons, or quark antiquark pairs. When matching onto SCET we perform a one-step 
matching and directly match full QCD onto SCET, as was done for H + 0 jets in Ref. [qt| . Most QCD results are 
obtained in the limit of infinite top quark mass, by first integrating out the top quark and matching onto an effective 
ggH interaction, 


£hard= l ^HG%G^ a , (Bl) 

which is then used to compute the QCD amplitudes. Here v = (\/2G_f) - 1 / 2 = 246 GeV. From the point of view 
of the one-step matching from QCD onto SCET, using Eq. (IBID is just a convenient way to compute the full QCD 
amplitude in the mt 00 limit. In particular, the a s corrections to Ci in Eq. (IB 1 ll are included in the amplitudes 
below, and therefore also in the SCET Wilson coefficients. In this way, if higher-order corrections in 1 /mt or the 
exact mt dependence for a specific amplitude are known, they can easily be included in the QCD amplitudes and the 
corresponding SCET Wilson coefficients. We illustrate this for the case of H + 0 jets below. 

We separate the QCD amplitudes into their IR-divergent and IR-fmite parts 

A = Adi v T Afi n , 

B = I?div + Ba n , (B2) 

where Afi n , Ba n enter the matching coefficients in Sec. [V] For simplicity, we drop the subscript “fin” for those 
amplitudes that have no divergent parts, i.e. for Adi v = 0 we have Ag n = A. For the logarithms we use the notation 

L a = ln (-^r - i0 ) > L n/H = ln (~^r - i0 ) - ln ('^ - i0 ) ■ 
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1. H -\- 0 Jets 


We expand the amplitudes in powers of a s (fi) as 

2 T F a s {n) 


A = 


E 


37 TV 


(B3) 


n —0 


The amplitudes with opposite helicity gluons vanish to all orders because of angular momentum conservation, 

A(l ± , 2=F; 3^) =0, (B4) 

corresponding to the fact that the helicity operators for these helicity configurations were not included in the basis of 
Eq. (11041) . The lowest order helicity amplitudes including the full m t dependence are given by 


_ Sl2 t 12 3 p(0) / 

' Sl 4 

- fii F (0)l 

' Sl \) 

2 (12) ' 

^4mf7 

2 ' 

< 4 m 2 / 

_ «12 (12) p( 0 )f 

' Sl \) 

- fH F (o)l 

' Sl2 4 

2 [12] ' 


2 ' 

v. 4 m\ ) 


i<£_ j- 


where the function F^(z) is defined as 

F (°U Z ) = — - — i _ I / arcsin2 (\/^) - _ 0 < z < 1, 

2z 2z z ^ In 2 [—i(^/z + \Jz — 1)] , z>\. 

For simplicity, we have extracted the (irrelevant) overall phases 

i$ ++H __ ilHl $—h _ 


(B5) 


(B6) 


(B7) 


( 12 ) [ 12 ] 

Since the two helicity amplitudes for ggH cannot interfere and are equal to each other by parity up to an overall 
phase, their higher-order corrections are the same as for the spin-summed ggH form factor. The divergent part of the 
NLO amplitudes is given by 


41;( 1± ’ 2± ;3tt) = ^ (0) (l ± 1 2 ± ;3 ff ) 


-^Ca + -(2C a L 12 ~/3 0 ) 


(B8) 


The IR-finite parts entering the matching coefficients in Eq. (1 1071) at NLO are 


4n(l ± .2 ± ;3 ff ) = 4l( 0 )(l ± ,2 ± ;3 ff ) 




FW{z) = C a (5--z 

+ C F ( 


4m 

1289 _ 2 155 _ 3 5385047 

4725 ^ _ 1134" ~ 65488500 


*) 


„ 307 25813 9 3055907 

—3 +- z 4-- z -\- - 

90 18900 3969000 


659504801 4 n 
1309770000 ) K ’ 


(B9) 


The full analytic expression for F^(z) is very long, so we only give the result expanded in 2 . Since the additional m* 
dependence coming from F^\z) is small and the expansion converges quickly, the expanded result is fully sufficient 
for on-shell studies of Higgs production. The IR-finite parts at NNLO are [93] 

4 n (1*, 2*; 3 H ) = A<°> (1±, 2 ± ; 3 ff )| \c 2 a L{ 2 + \c A p 0 L\ 2 + C A [c A (~ + y ) - §A> - ^ (1) (y|) 


3 

Ci(f - 2 Cs ) + - y) - + ^'(H) 

F< 2 >(.-) = (7 C\ + UCaCf - 6 C F ffo) ln(-42 - iO) + C \(-^ + E- + ^ - 44C S ) 

+ CAOr( .f _ £ + 44C3) + + ^ + ^)- IC.T F 

27 o /41 \ 4 

+ —Cp + CfP o — 12(3 J — -CfTf + 0{z). 


(B10) 


Here we only give the leading ter ms i n the m t —> 00 limit. The first few higher-order terms in z in F^ 2 \z) can be 
obtained from the results of Refs. Ezlllzi. 
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2. H + 1 Jet 

The amplitudes for H + 1 jet were calculated in Ref. [32] in the m t —» oo limit. Reference ^2j uses Tp = 1 and 
g s T a /y /2 for the ggg coupling. Thus, we can convert to our conventions by replacing T a —» -\/2T a , and identifying 
1/N = Ca — %Cf and AT = Ca in the results of Ref. [32|. We expand the amplitudes in powers of a s (n) as 

A= 2T^ sM ^ A(n ,^y (Bll) 

n —0 


a. gggH 

The tree-level amplitudes entering the matching coefficient C++± in Eq. (11151) are 


A ( 0 ) (l + ,2 + ,3 + ;4 H ) = ~^= 


m 


V% (12)(23)(31) yj 2 1 S 12 513^23 


H g i'J > +++n 


A( 0) (l + ,2+ 3-;4 ff ) = 


v/2 [13][23] Y / 2|si2Si3S23 

where we have extracted the (irrelevant) overall phases 
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vW a/I s 13| \/1S231 
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(12) (31) (23) ’ 

The divergent parts of the one-loop amplitudes are 


[12] \/\sr2\ \/|ai31 \Z|523| 
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(B12) 


(B13) 


(B14) 


f(s 12 ,s 13 ,s 23 ,m 2 H ,g) = - C A 


12 + -^13 + ^ 23 ) + Li 2 /hLi 3 /h + L12/HL23/H + Li 3 /hL 23 /h 
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-3 Cp. 


The finite parts of the gggH amplitudes, which enter the matching coefficient C++± at one loop are 

4^(1+ 2+3+;4 ff ) = A( 0 )(l+,2 + ,3+;4 ff ){/( Sl 2, S i3,S23,m^, M ) + \{C A - 2T F n f ) SuSl3 + ' Sl2 ; 23 + 513523 

l 6 m H 

4n(l + ,2 + ,3-;4 ff ) = ^°)(l + ,2+,3-;4 ff )|/( Sl 2,s 1 3,S23,m^, / r) + l(C A - 2 T F n f ) j , (B15) 

where we have extracted the common function 
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(B16) 


b. gqqH 

The tree-level amplitudes entering the matching coefficient C±(+) in Eq. (11141) are 


A(°)(l+;2+3r;4 H ) = -4= ^ 


v/2 [23] \/2| S 23 


Sl2 e i*+(+)n ) 
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The divergent parts of the one-loop amplitudes are 

4il(l ± ; 2 ^3r;4 H )=^ (0) (l ± ;2+ 3~-4 h ){-L(C a + 2C f ) + 1 -[Ca(L 12 + L 13 -L 23 )+C f (2L 23 -3)-^ 


/3o 
2 J 


(B19) 


The finite parts of the gqqH amplitudes, which enter the matching coefficient C±(+) at one loop are 
A«(l+;2+ 3r;4*) = A™ (1+; 2+ 3^ ; 4 H )| 3 (s 12) s 13 , s 23 , m 2 H , g) + (C F - C A ) , 

4n(l _ ;2+,3r ; 4 ff ) = A(°)(l-;2+ 3f ; 4 H ){< ? (s 12 , s 13 , s 23 , m 2 ,, + (C F - CU) gj , 


(B20) 


where we have extracted the common function 
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(B21) 


3. H + 2 Jets 


The full set of tree-level hclicity amplitudes for H + 2 jets in the mt —> 00 limit were calculated in Ref. [3lj, and all 
amplitudes below are taken from there. We expand the amplitudes A, B , in the decomposition of Eq. (11191) . Eq. (11241) . 
and Eq. (11301) . as 


A = 

B = 


2 T F a s {g) 
3nv 

2 T F a s (g) 
3nv 


[9.m] 2 X> w (4t)"' 

n =0 

b.M]’i> w (^)”. 

n =0 


(B22) 


For simplicity, we only give explicit results for the tree-level amplitudes in this appendix. To reduce the length of 
expressions, we use the kinematic variables Sijk defined by 

Sijk — (Pi + Pj + Pk ) 2 = Sij + Sik + Sjk ■ (B23) 


The H + 2 jets process is nonplanar, which means that we cannot remove all the relative phases in the amplitudes. 
It is therefore most convenient to keep all expressions in spinor helicity notation. We will explicitly demonstrate an 
example of the phases which appear in Eqs. (IB28I) and (IB29I) . 


a. qq q'q'H and qq qqH 


The tree-level amplitudes entering the Wilson coefficients (7( +; ±) and C{+±) in Eqs. (1 121 D and (I122D are 
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(B24) 


b- ggqqH 


The tree-level amplitudes entering the Wilson coefficients (7 + -(+)) C' + +(+), and C__i + \ in Eq. (11261) are 


A(°)(l+,2-;3+,4r ; 5 ff ) = 


(24) 3 
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^ ( 0 ) (2-,l+;3+ 4r ; 5 ff ) 
A(°)(l+,2+;3+ 4r ;5j? ) 
AW(l-,2-;3+ 4r;5 ff ) 
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(B25) 


In these expressions we have eliminated the Higgs momentum, p$, using momentum conservation, so that all momenta 
appearing in the above expressions are lightlike. We have also used an extended spinor-helicity sandwich, defined by 
[i\j + k\l) = [ i\j\l) + [i\k\l) to simplify notation. 

All the B amplitudes vanish at tree level, 


B (0) (l+, 2"; 3+ 4r ; 5 H ) = H (0) ( 1+, 2+; 3+ 4r ; 5 H ) = 2“; 3+ 4r ; 5 H ) = 0 . 


(B26) 


c- ggggH 


The tree-level amplitudes entering the Wilson coefficients C + _|_, C++h _, and C++++ in Eq. (11331) are 

A(°)(1 + ,2 + ,3 + ,4+;5h) = ~ 2Mr 
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(B27) 


To illustrate the relative phases that appear in these amplitudes, we can rewrite the amplitude 
#>(1 + , 2 + , 3 _ ,4 _ ; 5 h) in terms of the Lorentz invariants Sij 


A ( 0 ) (l+,2 + ,3-,4-;5 ff ) = 2e ^++— 1 


5 12 


5 34 
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(B28) 


with 


_ 2/3 ar I i\/^[~ Sl2S 34 + s 13 s 24 + S 14 S 23 ~ Kv^ + 2v / Sl3\/'S23 s 14)] 1 

\ —Si 2 S34(- v /Si3 — i^/sis) + (si3S 2 4 — Si4S 23 + i\/a)(- v /si3 + iy / s 23 ) / ’ 

a = 16( P 3 P 4 ) 2 = 4 S 13 S 14 S 23 S 24 - (S 12 S 34 - Si 3 s 24 - si 4 s 23 ) 2 > 0 , 

/3 = sgn (e^pap^p^pl ). (B29) 


The branch cut of the square root is given by the usual prescription, y/s/J = Sij + iO = iy/JsiTT if s v < 0. Our 
convention for the antisymmetric Levi-Civita tensor is €0123 = —1. For this process we can choose a frame where all 
but one of the momenta pi through P 4 lie in a plane (with p 3 determined by momentum conservation). The phase tp 
is needed to determine the momentum of the nonplanar momentum and the sign [3 resolves which side of the plane 
this particle is on, which is not captured by the (because they are symmetric with respect to a reflection about 
the plane). We note the simplicity of the spinor-helicity expression as compared with the explicit expression for the 
phases. 


Appendix C: Helicity Amplitudes for Vector Boson + Jets 

In this appendix we give all required partial amplitudes for the vector boson + jets processes discussed in Sec. IVII 
For each of the amplitudes A q v ^ a . B q , va defined in Sec. ED we split the amplitude into its IR-divergent and IR-finite 
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parts, 


X — Xdiv + Xfi n , 

where X stands for any of A qiVia and B q Vta . For the logarithms we use the notation 


Lij = ln(-^ - io) , L ij/k i = Lij - L u = ln( - iO ) - ln( -M - : 




(Cl) 


(C2) 


as 


1. V + 0 Jets 

In this section we give the amplitudes A qtVta for V + 0 jets. For each partonic channel, we expand the amplitudes 

Mm)' 


X - . 


(C3) 


n —0 


where X stands for any of R 9i „ j0 . The tree-level and one-loop helicity amplitudes entering the matching coefficient in 
Eq. (11461) are given by 
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(C4) 


2. V + 1 Jet 

In this section we give the amplitudes A q ^ v ^ a for V + 1 jets. Each amplitude is expanded as 

'Mm)' 


X = 


»wS^(T)' 


n =0 


(C5) 


where X stands for any of A g ^ >a . The tree-level and one-loop helicity amplitudes for V + 1 jets were calculated in 
Refs, pm, US EH. We use the results given in Ref. 0 , which uses Tp = 1 and g s T a /V 2 for the qqg coupling. We 
can thus convert to our conventions by replacing T a —»• \/2T a , and identifying 1/X = CU — 2 Cf and X = Ca ■ The 
one-loop amplitudes are given in the FDH scheme in Ref. 0 , which we convert to the HV scheme using Eqs. (11001) 
and (11011) . 

The tree-level amplitudes entering the matching coefficient <7 x+ (+;+) i 11 Eq. (11541) is given by 


M)(i+;2+,3r;4+,5r) = ^2v / 2 


<35} 2 
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The divergent part of the one-loop helicity amplitude is given by 
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The finite parts entering the matching coefficients at one loop are 
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| ~Y~ - ^13 + 31/13 - 7 + — ^ + ( Cf y-j (-7/23 + 37/45 - 8 + — j 
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The contributions from virtual top quark loops are calculated in an expansion in 1 /rrit to order 1/m 2 in Ref. [l4j |. 
hence the divergent behavior of as m* —> 0. To reduce the length of the expressions, we have used the commonly 
defined functions 


L o(r) = 


In r 
1 — r 


Li(r) = 


Lp(r) + 1 
1 — r 


Ls_i(ri,r 2 ) = Li 2 (l — n) + Li 2 (l — r 2 ) + lnri lnr 2 — 
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(C9) 


The proper branch cut of logarithms follows from the prescriptions s^ 
dilogarithm follows from that of the logarithm through the identity 


Sij + iO. The proper branch cut of the 


Im [Li 2 (1 — r)] = — ln(l — r) Im [ In r] . 


(CIO) 


3. V + 2 Jets 

In this section, we give the amplitudes A qiVa , B qtVta for V + 2 jets. Each amplitude is expanded as 

v=b.wi 2 ften) 

n= 0 

where X stands for any of A qiV ^ a or B q , v ^ a . We also define the kinematic variables Syfc as 

Sijk = (Pi + Pj + Pk) 2 = Sij + Sik + Sjk ■ (C12) 

The one-loop helicity amplitudes for q'q'qqV and qqqqV were calculated in Ref. [l86l |. The one-loop helicity ampli¬ 
tudes for ggqqV were calculated in Ref. )14j j. which also gives compact expressions for the four-quark amplitudes, 
which we use here. The contributions from virtual top quark loops are calculated in an expansion in 1/mt to order 

l/? 7 i 2 in Ref. El- 

Reference El uses Tp = 1 and g s T a /\/2 for the qqg coupling. We can thus convert to our conventions by replacing 
T° — > \[2 T°, and identifying 1/N = Ca — 2Cf and N = Ca- The one-loop amplitudes are given in the FDH scheme 
in Ref. (l3| . which we convert to the HV scheme using Eqs. (11001) and (11011) . 

a. q'q'qqV and qqqqV 

The tree-level amplitudes for q'q'qq V and qq qq V entering the Wilson coefficients in Eqs. (11731) and (11741) are given 
by 


4® (1J, 2,7; 3+, 4r ; 5+, 6 J) = - B <°> (1+, 2^; 3+ 4r; 5+, 67 ) 
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S12S56 
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(C13) 


Due to the length of the one-loop q'q'qq V amplitudes, we only show how to translate the decomposition of the 
amplitude in Ref. [Tl[ to our notation. The one-loop amplitudes are given in terms of the bare partial amplitudes 

A i;j (3q,2Q,lQ,4q ) of Ref. [ 3 ] as 


AW(V,2 r ;3„4 ff ;6+ 57 ) 
BW(l q ,,2 q -,;3 q ,4 g -;6+, 5J) 
AW(V,2 ? ;3 91 4,-;6+ 5J) 


-1327T 2 iV A 6 ; 1 (3 ? , 2g, 1 q, 4 9 -) - (y + 2 C f - 3„4 g -; 6 + 57 ), 

-1327T 2 iV A 6 ; 2 (3 9 , 2 q, 1q,4 9 -) - + 2C F - ^C A )B^{l q ,,2~ q '-,3 q Aq\ 6 +, 5^), 
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B^ = 0 . (C14) 


The overall factor —\32ir 2 is due to our different normalization conventions. We have not included helicity labels, as 
these relations are true for all helicity combinations. Note that the partial amplitudes do not include labels for 
the lepton momenta, which are implicitly taken as 6^, 57 . The terms in the first two lines proportional to A q °^ and 
Bq °' 1 come from the UV renormalization and switching from FDH to HV. 


b- gg qq v 


The tree-level amplitudes for gg qq V entering the matching coefficients in Eq. (I182D are given by 
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Due to the length of the one-loop gg qq V amplitudes, we again only show how to translate the decomposition of 
the amplitude in Ref. flij to our notation. The one-loop amplitudes are given in terms of the bare partial amplitudes 
A l - J (3 q , 1, 2, 4q). Al :j (3 q , 4q-; 1, 2), and A%(3 q , 4 q ; 1, 2) of Ref. [3 as 


AfH 1, 2 ; 3 q , 4 q - 6 + 57 ) = -i 64w 2 N A 6;1 ( 3„ 1, 2 , 4 q ) - (^ + C F ) A<°> ( 1 , 2 ; 3„ 4 q ; 6 + 57 ), 

(1, 2; 3„ 4 q - 6 + 57 ) = -i 64 tt 2 A 6 ; 3 (3 g , 4 q ; 1,2), 

4 1 ] (1,2; 3„ 4 q \ 6 + 57 ) = -i647r 2 Ag ; 4 (3„ 4 q ; 1,2), 

Bp( 1, 2; 3„ 4q- 6 + 57 ) = +i647r 2 1 A^. 4 (3„ 4 ? ; 1, 2), 

A^l, 2; 3„ 4 q - 6 + 57 ) = -i64 tt 2 A^(3 ? , 4 q ; 1,2), 

BW( 1, 2; 3„ 4 q - 6 + 57 ) = -A64tt 2 A [^(3,, 4 q ; 1 , 2 ) - A“ 4 (3 9 ,4,-; 1 , 2 ) - A^(3 q , 4 q ; 2 , 1 )] . (C16) 
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The overall factor —i647r 2 is due to our different normalization conventions. We have not included helicity labels, as 
these relations are true for all helicity combinations. Note that the partial amplitudes A i; j do not include labels for 

the lepton momenta, which are implicitly taken as 6jl~, 5^T. The term in the first line proportional to comes from 
the UV renormalization and switching from FDH to HV. 


Appendix D: Helicity Amplitudes for pp —► Jets 
1. pp —> 2 Jets 

In this appendix we give explicit expressions for all partial amplitudes that are required in Eqs. (11891) . (11901) . (11941) . 
and (I203D . for the various partonic channels of the pp —> 2 jets process. Since this process is planar, we can write all 
amplitudes for a given set of helicities with a common overall phase extracted, which is determined by the phases of 
the external particles. In this way, we do not need to worry about relative phases between the Wilson coefficients for 
different color structures when they mix under renormalization. The cross section does not depend on this overall 
phase. This simplifies the numerical implementation considerably for this process, as it avoids having to implement 
the complex spinor algebra. To extract the overall phase from the amplitudes, the following relations for the relative 
phases between the spinor products are useful: 

(12) ^ (34) = (14) = (23) = (13) = (24) 

[34] [12] [23] [14] [24] [13] ' 1 ’ 

These relations follow from Eq. (IA22I) with n = 4. 

We split the partial amplitudes into their IR-divergent and IR-finite parts, 

A = Aiiv + Aa n , B = .Bdiv + Bn n , (D2) 

where the IR-finite parts enter the matching coefficients. We expand the amplitudes and Wilson coefficients in powers 
of a s (p ) as 

X=[g s (p)} 2 f2x^( K ^y , (D3) 

n =0 

where X stands for any of Adiv.fln, -Sdiv.fin, and X^ and XA) are the tree-level and one-loop contributions, respec¬ 
tively. For simplicity, we drop the subscript “fin” for those amplitudes that have no divergent parts, e.g., for the 
tree-level amplitudes A^l = 0 and Ag^ = A^°K For the logarithms we use the notation 

La = ln(—^ — io) , iy /fc i=Lii-£fci=ln(-^-i0)-ln(-^-i0). (D4) 


a. qqq'if andqqqq 


Here we list all partial amplitudes up to one loop entering the Wilson coefficients in Eqs. (11891) and (11901) . The 
one-loop helicity amplitu des for qq q’q' and qqqq were first calculated in Ref. [23], and the two-loop helicity amplitudes 
were computed in Refs. 


We find agreement between the one-loop results of Refs. 19l| and [1921 ]. from 


which we take our resultsjjj Our one-loop matching coefficients agree with the calculation of Ref. [881 ]. 
The tree-level amplitudes are 


A(°)(l+, 21- ; 3+, 4~) = —B<°>(1+ 2r; 3+, 4r,) = - 


(24) [13] 


512 


£13 P i* 

512 


(+;+) 


10 Note that there is a minor disagreement here with the earlier 
calculation in Ref. [23| . presumably due to typos. Specifically, 
in Ref. [23| the factors (log 2 + 7r 2 ) and (log 2 + 7r 2 ) in 


Eqs. (5.10) and (5.12) respectively, m ust be sw apped to ac hieve 
agreement with the results of Refs. [8 Rll9lMl92(| . Reference [192(| 
also has a minor typo, having a flipped overall sign for the IR- 
divergent terms. 
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q ’ q ’ q' 1 q 

where the phases are given by 


A (0) (l+ 2r;3-,4+) = -B(°)(l+ 2r;3-,4+) = _MM = £il e ^( +; -) , 


S 12 


S 12 


e i4>( +;+ ) _ 


(24) 

(13)’ 


e i4> (+;-) 


(23) 

(14) 


(D5) 


(D 6 ) 


We have chosen to express all the one-loop amplitudes in terms of A ( 0 ^(l+, 2^ ; 3^,4-,) and A^°)(l+, 2^ ; 3 q/ , 4^). The 
divergent parts of the one-loop amplitudes are 

; 3^,4-,) = A ( 0 ) (l+, 2g ; 3j, 4-,)| —— Cf + ^ [C(f( 2 £i 2 — 4 £ 13 / 14 — 3) + C A {L 13/14 — £12/13)] j > 
-®div(-*-g ) ; 3j, 4-,) = 7l ( 0 )(Tq", 2^ ; 3+,, 4-,) |^Cf — - [Cf( 2£ 4 2 — 2L 13 / 14 — 3) + C A (L 13 / 14 — £12/13)] | , 

4 il(-*-g~ 1 2 , ; 3 ?/ , 4 ^) = ^ 0 ) ( 4 , 2 ^ ; 3 ^,, 4 ^) | — Cf + - [Cf( 2 £i 2 — 4 £i 3 /i 4 — 3 ) + Ca (£13/14 — £12/13)] j > 

£^(l+,2r; 3 -„4+) = AW( 1 + 2r ; 3 -,4+) ||(7f - ^[Cf( 2£ 12 - 2£ 13/14 - 3) + C A (L 13/U - L 12/13 )\ J . (D7) 

The finite parts entering the Wilson coefficients are 

4 n ( 1 +, 2 r; 3+ 4y) = A(°)(l+ 2 r; 3+ 4y) [/(a 12j s 13 , su, g) + ( 4 C F - C A ) g(s 12 , s 13 , *i 4 )] , 

£^( 1 + 2r ;3 + 4r ; ) = A<°\ 1+ 2r ;3 + 4 =) [ 4 C F £i 2 £i 3 /14 - /(ai 2 ,«i 3 ,«i 4 ,A*) + (C A - 2C F ) g(s 12 , s 13 , si 4 )] , 


4]] (1+ 2r ; 3",4+) = ^ (0) (1+, 2r ; 3 -, 4+) [/(si 2 , a 13 , s 14 , g) + 2 (C A - 2 C F ) g(s 12 , s 14 , s 13 )] , 

4n (lj, 2=-; 37 ,, 4+) = A ( 0 ) (l+ 2r ; 3 -, 4+) [4CV L 12 L 13/14 - f(s 12 , s 13 , s 14 , g) + 2{C F - C A ) g(s 12 ,s 1A , S13)] , 

10 


/(s 12 ; ^i 3 , Si 4 , g) — Cf — 2 £^ 2 + 2£i2(3 + 4£ 13 / 14 ) — 16 + — + C A ( 2 £i 2 (£12/13 — £13/14) + + 71-2 ) 


J (£ 12 _ I ) ’ 


S , (si2,Sl3,Sl 4 ) = 


£12 

Sl3 


(l - ^) (£?2/14 + 7r “) + £l2/14 


(D 8 ) 


b- ggqq 

The one-loop helicity amplitudes for gg qq were first calculated in Ref. ;23|, and the two-loop helicity amplitudes 
were computed in Refs. jl87l Il88| . We take our results from Ref. m 3 , converted to our conventions F 4 ! 

Here we list all partial amplitudes up to one loop entering the Wilson coefficients in Eq. (119411 . We start with 
the partial amplitudes where the gluons have opposite helicity, which are the only ones having a nonzero tree-level 
contribution. The tree-level amplitudes are given by 


H(°>(1+ 2-;3+ 4r) = -2 

A(°)(2-,l+;3+ 4r) = -2 
£?(°)(l + ,2-;3+4r)=0. 


(23) (24) 3 


(12) (24) (43) (31) 
(23) (24) 3 
(21)(14)(43)(32) 


= 2 V^ 13514 ! e i4>+-(+) 


S12 


= 2 S 13\/|S13S14| e i$ + _ (+) 


s 12 514 


(D9) 


11 We find a slight disagreement with the earlier results of Ref. [2^1 
for their subleading color amplitude in Eq. (5.24). This ampli¬ 
tude appears to have typos since it does not have the correct IR 


structure, as determined by the general formula [215f or by the 
SCET result in Eq. (|G4l) . Comparing with the matching calcula¬ 
tion of Ref. [ 88 |L we find a typo in the 7 r 2 term in W 4 in Eq. (54), 
which should have 37r 2 u 2 /(2ts) —>■ —37r 2 u/(4t). 
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In the second step we extracted a common overall phase from the amplitudes, which is given by 

(24) [13] [14] 


e i$ +-<+) = 


[24] ^/|S13 S 14 

The divergent parts of the corresponding one-loop amplitudes are 
4w( 1+ > 2"; 3+ 4r) = A<°>(1 + , 2"; 3+ 4r) 

(2 \ 1+; 3+ 4r) = 1+; 3+ 4r) 


(DIO) 


o (Ca + Cf) H (2 C F £12 + 2CU L13 — 3CV — /?o) 


—"o (Ca + C F ) H (2 C F £12 + 2CU L14 — 3CV — /?o ) 


(1+, 2"; 3+ 4r) = #(lt; 2“; 3+ 4r) ± 4T F (i 12/14 + ^ £12/13) • 

The corresponding finite parts entering the Wilson coefficient C + _(+) at one loop are 

^lin (^ + > 2 ;3+,4g) = A (0) (l + , 2 ; 3+, 4^ )|c(4 L\ z + £ 42 /13 + 1 4 g - ) + Cf £12 + 3£i2 — 8 + 

+ (Ca — Cf )— (£12/13 + 7r2 )l > 

Sl 4 J 

" 4 fin (2 -) 1 + ; 3 g , 4 g ) = 2 ,1 + ; 3 +, 4 ^ ^-2Lj 4 + £42/14 - 3 £12/14 + 1 + 

+ C F (-£?2 + 3£i2 - 8 + y) - U [ l - U 

+ (it - cv ) S + S ii2/i4 ) ■ Li2/i4 + w 3 n \} ’ 


(Dll) 


6 / 


+ £l 2/14 + -r-K 2 
s l3 . 


4 1 n ) (l+, 2"; 3+ 4r) = A<® (1+, 2"; 3+ 4r) 4T> 


— £l2£l3/14 -|-£l 4 £l2/13-7 —“ (£ 12/14 + 7T ~) 

S14 '4 S13 7 


3 Si2 


(D12) 


The partial amplitudes where both gluons have the same helicity vanish at tree level, 

4l (0) (l + , 2+; 3+ 4r) = B<°) (1+, 2+; 3+ 4j) = 0 , 

24<°>(1-, 2"; 3+ 4r) = B(°)(1", 2"; 3+ 4j) = 0 . (D13) 

The corresponding one-loop amplitudes entering the Wilson coefficients C ++ ( + ) and C _( + ) are IR finite. They are 

A( 1 )(l+,2+;3+ 4j) = 2^^e i$ ++<+> \(C A -C F )— + \{C A - 2T F n / )—1 , 

L S13 6 Sl2 J 

AW(2+, l+;3+ 4r) = -2^^e i$ ++(+> \{C A -C F )—+ \[C A - 2T F n f ) — 

L S14 4 S12 

B< 1 >(l + ,2+;3+ 4r)=Q, 


(D14) 


and 


A«(l%2-;3+ 4r) = 2^^e i$ -<+> [(CU - C F )— + ±(CU - 27>n / ) — 

L Sis o Si2-i 


A«(2“, l-;3+,4r) = -2^^e i$ -<+> (C A - C F )— + -(C A - 2 T F n f ) — 

L S14 4 S12 


1 1 


1 


B«(l-,2-;3+4r) =0 , 
with the overall phases 


,'*++(+) = 


[12] [13] (14) 

( 12 ) y/\ 


e i*__ (+) = (12) [13] (14) 


s 13 5 14| 


[12] /I 


s 13 5 14| 


(D15) 


(D16) 
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c- gggg 


The one-loop heli city am plitudes for gggg were first calculated in R ef. 1 2 31 . and the two-loop amplitudes were 
computed in Refs. jl93l Il94l |. The results given here are taken from Ref. [l94j . and converted to our conventions. We 
also find complete agreement with the expressions given in Ref. [H1EEI 

The amplitudes inherit the cyclic symmetry of the traces, which means that many of the amplitudes appearing in 
Eq. (12031) are related, for example 


A(l+,3-,4-,2+)=2L(2+,l+,3-,4-). 


(D17) 


For the convenience of the reader, we will explicitly give all amplitudes needed in Eq. (12031) . We start with the partial 
amplitudes with two positive-lielicity and two negative-helicity gluons, which are the only nonvanishing amplitudes 
at tree level. We have 


A ( °)(l + ,2 + ,3-,4-) = 4 
A (0) (l + ,3 _ ,4",2 + ) = 4 
^(°)(l+,4-,2+,3") =4 
with the common overall phase 


(34) 4 


(12)(23)(34)(41) s 14 


= 4£H e i<* > ++— 


(34) 4 


= 4— e i3>+ 


(13) (34) (42) (21) s 13 

^ 34 ) 4 = j S 12 i4>++— 

(14) (42)(23)(31) s 13 s 14 


i<3> 


++-- — _ 


[12] (34) 


(12) [34] ' 

The corresponding B all vanish, 

b (°)( 1+,2+,3-,4-) = #(o)( l+,3“,4-,2 + ) = J3 (0) (l + , 4", 2+, 3“) = 0. 

At one loop the B^> amplitudes can be expressed in terms of the n/-independent part of the A^ l \ 

R (1) (l + ,2+,3 _ ,4“) = ^2[A (1) (l + ,2+,3",4“) + A (1) (l+, 3", 4", 2+) + A (1) (l+,4",2 + ,3 _ )] 


nf—0 


(D18) 


(D19) 


(D20) 


(D21) 


The same relation also holds for the other helicity assignments. Using the cyclic symmetries of the amplitudes, it 
follows that the last three entries in the Wilson coefficients in Eq. (12031) at one loop are all equal to each other and are 
given by 2 Tf/Ca times the sum of the first three entries at n/ = 0. The divergent parts of the one-loop amplitudes 
are 


4iv(l + > 2+ > 3”, 4") = #(1+ 2+ 3", 4“) 
4’(l + ,3-,4-,2 + ) = #(l + ,3-,4-,2+) 
A«(l+,4-,2+,3-) = A(°)(l+,4-,2+,3-) 
B';>(l+,2+,3-,4-) = #(l + ,2+3-,4-) 
4>(l + ,3-,4-,2+) = #(l+,2+,3-,4-) 
R^(l+,4-,2+,3-)=A(°)(l+,2+,3-,4-) 


-2 H-( Ca L±2 + Ca Li 4 — / 3 o) 

. 6 6 

—2 ^a H— {Ca L12 + Ca Ai 3 — / 3 o) 

e e 

—2 Ca 4 — {C A L\z + Ca A14 ~ Ad) 
e e 

' 8 T f 

e 


T , Sl4 T 

-^12/13 H-^12/14 

Si 3 


8 Tp I S 14 Si2 

- -^13/14 H-^13/12 

e \ s 13 S13 


8 T f 


J 14/13 


S -^L 

Sl3 


14/12 


(D22) 


12 We have also compared with the matching calculation of 
Ref. [88| , which has a minor typo. In particular, in u) in 

Eq. (61) the rif terms must be dropped and /3 q set to IICa/ 3. 


Also as noted in Ref. |89|| , the last column of Table 5 in Ref. I88(| 
applies to helicities 7, 8, while the second-to-last column applies 
to helicities 9-16. 
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The finite parts entering the Wilson coefficient C+H_at one loop are 

C A (—2 L 12 L 14 - l + A) (iu - |) 

4 47T 2 ' 


4 1 „ ) (l + ,2+,3-,4-) = AW(l + ,2+,3-,4-) 


4n(l + >3-,4-,2+) = A(°)(l + ,3-,4- ! 2+) 


Ca (- 2 l 12 l 13 - | + A>(iis - 3 ) 


4^(1+,4“, 2 +, 3“) = A(°>(l+4- 2 + 3"){CU (- 2 L 14 Li 3 + - A) (| + 


- (C A - 2 T F n f ) 


513^14 


5 12 


1 + ( £!i _ £! 1 , i ls/I4 
S12 S12 / 


+ 2 - 


S13S14 \ 


’12 


( 


r 2 1 2 

^13/14 + ^ 


,rp S13S14 

- 3r F n/—- 2 — 


’12 


(-^13/14 + 7f2 ) | 


4 (1 + , 2 + , 3-, 4-) = flW (1+, 3-, 4-, 2+) = 441+, 4-, 2+, 3-) 
= -4T F y4 ( 0 ) (l + ,2 + ,3~,4“) 


5 ]_4 ci t t i o r r 

514 

S14 

<Si 3 

Sl4\ 

- Zivi3 -^12/14 i ^14 -^12/13 i 


+ — 



.513 

Sl2 

Sl2 

^Sl2 

Sl2' 


+ 


S14 

Sl 2 


( 2 iZ-Z) ( L 13/14 + 7I ' 2 ) 


(D23) 


Due to Eq. (ID 1 71) . the first two amplitudes in Eq. (ID18I) . as well as the first two in Eq. (ID23I) . can be obtained from 
each other by interchanging l" 1 " -o- 2 + which corresponds to S 13 O S 14 without an effect on the overall phase. 

The amplitudes with only one or no gluon with negative helicity vanish at tree level, 


A ( 0 ) (l + ,2+,3 + ,4 ± ) = A ( 0 ) (l + ,3+,4 ± ,2 + ) = A ( 0 ) (l+,4 ± ,2+,3 + ) = 0, 
B ( 0 ) (l + ,2+,3 + ,4 ± ) = B ( 0 ) (l + ,3 + ,4 ± ,2+) = B ( 0 ) (l+, 4 ± , 2+, 3+) =0. 


(D24) 


The corresponding one-loop amplitudes are infrared finite. Those entering C+-H_ are given by 

4<‘>(1 + ,2+3+4-) = 4 p J| 5M i(C, - 2 T r n,) (»„ + , M ) = 4e— i{C A - 2T r „,) (g + |f) , 

4P)(l + ,3+,4-,2+) = i«7, - 2T r n f ) (» IS + »„) = 4e»+++- \(C A - 2 T F n f ) (g + |f) , 

• 4< ‘ , (1 + . V,2 + ,3 + ) = - 2T r n f ) (» I3 + ,u) = 4e»— l (C A - 2T r n,) (JH + |f) , 

B ( 1 ) (l + ,2 + ,3 + ,4") = B ( 1 ) (l + ,3 + ,4",2 + ) = J5 (1) (1+, 4“, 2+, 3+) = -16 T F e^+++~ , (D25) 


and those for C++++ are 

24 (1) (l + ,2 + ,3 + ,4+) = A ( 1 ) (l + ,3 + ,4 + ,2 + ) = A ( 1 ) (l + ,4 + ,2 + ,3 + ) = 4e i4> ++++ \{C A - 2 T F n f ), 

o 

B ( 1 ) (l + ,2 + ,3 + ,4+) = B ( 1 ) (l + ,3 + ,4 + ,2+) = B ( 1 ) (l+, 4+, 2+, 3+) = 16 Tp e i<s, ++++ , 
where for convenience we have extracted the overall phases 

e i* +++ _ = M M M e i* ++++ = _ M M 

(12) (13) [14] ’ (12) (34) • 


(D26) 


(D27) 


2. pp —f 3 Jets 

In this appendix we give explicit expressions for all partial amplitudes that are required in Eqs. (12091) . (12101) . (12171) . 
and (|226L for the various partonic channels for the pp —> 3 jets process. The one-loop amplitudes for these processes 
were calculated in Refs. [HI, [24j, HH, respectively. These papers use Tp = 1 and g s T a /y/ 2 for the qqg coupling. Thus, 
we can convert to our conventions by replacing T° —> \/2 T a , and identifying 1/iV = C A — 2 Cp and N = C A . Below 
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we restrict ourselves to giving explicit expressions for the tree-level amplitudes, since the one-loop expressions are 
fairly lengthy. 

For each partonic channel, we expand the amplitude as 

OO / \ 

A' = [5 S (M)] 3 E A ' (n) (^r)”’ (D28) 

n —0 


where X stands for any of Adi v ,fin> -Bdiv.fin- 


a■ gqq q'q' and g qq qq 


The tree-level amplitudes entering the Wilson coefficients in Eqs. (12091) and (12101) are given by 

= 3 f ) = -V2 

/1»>(C;2+ 3r;4-„5j) = 5j;2+ 3J) = V~2 


(13) (14) (23) (45) 
(34) 3 


f?(°)(l+;2+ 3r;4+,5r,) = - v / 2 


£?W(l + ;2+ 3r;4-„5+) = V2 


(12)(15)(23)(45) 
(23)(35) 2 
' (12)(13)(23)(45) 
(23) (34) 2 


B (0) (l+;4+,5r,;2+3r) = -y2 


(13)(14)(23)(45) 
(35) 2 (45) 


BW(l + ;47„5+;2+,3r) = v^ 


(14)(15)(23)(45) 
(34) 2 (45) 


(12)(13)(23)(45) v (14)(15)(23)(45) 

Of these helicity amplitudes only 4 are independent. The one-loop amplitudes were computed in Ref. [2~i[ 


(D29) 


b- gggqq 

The three independent tree-level partial amplitudes which enter the Wilson coefficients in Eq. (12171) are given by 

(34) (35) 2 


A (0) (l+,2+,3 ;4+5 g -) 2v/ 2-( 12 )(i4)(23)(45) ’ 

7LW(2 + ,3-,1+;4+ 5 g ) = -2v / 2- 


(34) (35) 3 


'(13) (15) (23) (24) (45)’ 


A (0) ( 3-,l + ,2+;4+ 5 g ) = -2^ 


(35) 3 


“ v “(12)(13)(25)(45) ’ 
B (o) = C (0) = 0 _ 


(D30) 


At tree level, the partial amplitudes for the other color structures vanish, = C^ = 0. The one-loop amplitudes 
were computed in Ref. [25;]. 


c- ggggg 


The two independent partial amplitudes that enter the Wilson coefficients in Eq. (12261) are given by the Parke-Taylor 
formula [216l | 


A (0) (l+,2 + ,3 + ,4-,5“) =4V2 
A (0) (l + ,2+,4“,3 + ,5“) =4V2 


(45) 4 

(12) (23) (34) (45) (51) ’ 
(45) 4 

(12) (15) (24) (34) (35) ’ 


R (0) = 0. 


(D31) 


All other amplitudes can be obtained by cyclic permutations. The double-trace color structure does not appear at 
tree level, so = 0. The one-loop amplitudes were calculated in Ref. |22| . 
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Appendix E: RGE Ingredients 


In this appendix, we collect explicit results required for the running of the hard matching coefficients required to 
NNLL order. We expand the /? function and cusp anomalous dimension in powers of a s as 


OO 

/3(a s ) = -2a s y^j3 n ( 


Qs \ n+1 

4-7T/ 


OO 

Tcusp^s) = ^ ^ r n ( 


Os \ n+1 

47T/ 


(El) 


n =0 n—0 

Up to three-loop order in the MS scheme, the coefficients of the /3 function are [2rR (218i| 


A)=y C A -^T F n f , A = + 4U F ) T F n f , 


B _2S57 3 

P2 -~mT Ca 


(' 


205 


Cl - ^ C F C A - C 2 a ) 2 T f n f + (y C F + g C A ) 4T| . 


18 


(E2) 


and for the cusp anomalous dimension they are i2li l220l 

r 0 = 4, r, = (^-±d 

^ ^ + 44^ + 8 ^) C5+ ( y _ L _ ^ )ca2Tf „ i+ (B6 . ™ )Cr2Trni _ g ri 




To = 


2 2 
F n / • 


(E3) 


Note that here r cusp does not include an overall color factor; it differs from the usual qq case by a factor of C F . 
For the noncusp anomalous dimension of the Wilson coefficient, which is color diagonal to two loops, we write 


7(a s ) = K7c( a «) + n glc^ a ^\ 1 + °( a l) > 
as in Eq. (123811 . The quark and gluon noncusp anomalous dimensions, 

= (^)t'co + (y) lei, 7 &(<*,) = (y)7& 0 

have the following coefficients 

7c o = ~ 3 C F , 




(E4) 


(E5) 


7 q cl = ~C F 
7c o = “A>, 


41 


—— 26^3 C A + ( — — 27t + 24^3 ) C F + ( — + — ) /?o 


65 7T Z 


18 


59 


19 7U 


7ci — ( v ^y + 2 C3jC'i+^y + yJC'A/3o-/3i- 

The evolution kernels required for the resummation were defined in Eq. (12371) by the integrals 

O) r ( n \ r a ‘ \ 


(E6) 


Kr(no,n) = 


“»(«>) 
“«( m) 


Vr{po,p) = / da 
J a.(un) 


i b CUS p(o; s ) 

das ~57—r - 

P(a s ) 


(a*) 


£«) ’ 


cusp 


A' 7 (/r 0 ,M) = 


“s(Mo) 

*s(/i) 

da » 

»(«>) P( a 


/3(a s ) ’ 

7(« s ) 


Up to two loops, we can simplify the noncusp evolution kernel as 

A 7 (/x 0 ,/x) = (n q K*(no, ft) +n g K°{n 0 ,n)) 1. 
Explicit results to NNLL order are given by 

(l - i - tar) + (h - |)(1 -r + ln^r) + |-ln ! 


(E7) 


(E8) 
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a s (Mo) 


47r 


Vr{lMhV) = - 


£o_ 

2/3o 

A> 0 ,rt = - 2 | 


lnr - 

lnr ■ 

lnr ■ 


3? 

3 0 2 A> 

« s (mo) 

47T 

« s (/^o) 

47T 

Q s (Mo) 

47T 


)(¥*.).<»-» 


2 

£i _ A 

To Pa 


(1 — r + r lnr) — 

PiTi 


t ~i\ , a 2 (Mo) /r 2 PiTi pi p 2 \ 

y 167t 2 v r o A,r 0 /3 0 2 /V 


^_Ar^(i-r) 

r 0 

r 2 — 1 


21 


Adm/ 


[A 

( 1 - — ) 

A 

flnX 

+ 1-1)1 

[Pa 

K X) 

pi 

K X 

X )\ 


with r = a s (^)/a s (fj, o). The running coupling in the above equations is given by the three-loop expression 

—^ = f + lnX + 

a s (W a s (no) 47t/3 0 I67H 

with X = 1 + a s (/i 0 )/3olii(Ai/A i o)/(27r). 

Appendix F: Color Sum Matrices 

For each specific process considered in the text we de¬ 
composed the Wilson coefficients in a color basis as 


I- 

For gg qq in the basis Eq. (11921) , we have 

C A C F 


T - — 

± ssm ~ 


s-iai ■■■a n 
°+-(— 


E ^jk Jifll'-Qn — rj-iai---Oi n pj 

+ "( ) & 




(FI) 

where T a 1 "' (Xn is a row vector of color structures which 
form a complete basis of the allowed color structures for 
the particular process. Since convenient color bases are 
generically not orthogonal, the scalar product between 
Wilson coefficients is nontrivial. The Ct is given by 

= [< C 111 " 01 "]* T ai '"“” =C* T T, (F2) 



and for gggg in the basis Eq. (11971) . we have 


T — 

± 9399 ~ 


C A C F 


where 


p _ ^ ^ ^pai--a n ypa 1 ---a n (F3) 


(a b b c d c\ 
b a b c c d 
b b a d c c 
c c d e f f 
d c c f e f 
\c d c f f ej 


is the matrix of color sums. 

In this appendix we give explicit expressions for T for 
all the processes in this paper, both for general SU(N), 
as well as a numerical result for the specific case of IV = 3. 
For simplicity, in this section we restrict ourselves to the 
normalization convention Tp = 1/2, and Ca = N, and 
write the results for general SU(N) in terms of only Ca 
and Cf- 

For qq and gg in the basis in Eq. we have 


where 


a = C\ - -C A C F 


6 C| 


23 

12 


b = C\ — 5 C A C F + 6 C 2 f = -- 


1 


c = C F = - , d = 
e = C F C A = 4, / = 


(2 C F -C a ) 1 


6 


T qr] = C A = 3 , T gg = 2 C a Cf = 8. 

For gqq and ggg in the basis Eq. (11101) . we have 
fgq q = C A C F = 4, 

T 999 = Wf 


(F4) 
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(E10) 


2 Cp 2 Cf ~ Ca 1 
2Cf ~ Ca 2Cf 1 

1 1 CV 


(F7) 
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For g qq qq and g qq q'q' in the basis Eq. (I206D we have 
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For qqqq and qqq'q' in the basis Eq. (11861) . we have 


Tqqqq ~ Tqqq'q' 
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For ggg qq in the basis Eq. (12151) we have 
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where 

a = 4 C a C 2 f = y , b = C A ~ 2C f = ± , 

10 8 

c=(C 2 A + l)(C A -2C F ) = -, d = -2C F = -~, 

e = -1, / = 2 C a Cf = 8, g = 2 C\C F = 24, 
h = C A = 3, i = C\- 2 = 7, j = —2 . (F12) 


For ggggg in the basis Eq. (I221[) we have 
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a = C\ - 4 C\ + 10 = 55 , b = 2C\ - 4 = 14 , 

c = 2, d = 2 C\C F = 24, e = C A = 3, 

f = 2C 3 A C F = 72, g = C 2 A = 9. (F15) 


Appendix G: IR Divergences 


In this appendix, we explicitly check that the IR diver¬ 
gences of QCD are reproduced by SCET. This ensures 
that they drop out in the one-loop matching, and that 
the resulting Wilson coefficients are IR finite. They also 
provide a very useful cross check when converting from 
the different conventions used in the literature to ours. 

The one-loop matching equation relating the SCET 
operators and their Wilson coefficients to the QCD am¬ 
plitude is 


where 
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First we determine the residues of the propagators en¬ 
tering the LSZ reduction formula. Regulating both UV 
and IR divergences in dimensional regularization, all bare 
loop integrals in SCET are scaleless and vanish, i.e. the 
UV and IR divergences cancel. In particular, for the self¬ 
energy diagrams, we have 

E = E uv + Sir = 0. (G2) 


The UV divergences Euv plus possible additional UV 
finite terms T, x (as dictated by the renormalization 
scheme) determine the wave function renormalization Z £. 
The remainder Eir — E^, enters the residue 


ry-\ _ -] d(E UV + E x ) 

^ ^ d^ 

n -i i d(EiR — E x ) 

R e =1 d^- 


0 


(G3) 


At one loop in pure dimensional regularization, we then 
have = Z, and similarly for gluons R A = Z^ 1 . In 
the on-shell scheme E x = Em, so with pure dimensional 
regularization Z^ = = Z A = R A = 1. 
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Since all loop diagrams contributing to (CP)*- 1 ) vanish, we find 
the only nonzero contributions come from the countert¬ 
erm in Eq. (12291) and the one-loop residues. At one loop 

I 


(ot)( i )(7(° ) = (d^)^\(z^ q/2 z^ s/2 z c - i) + - i)jc (0) = (d t } (0) (z c - i)c (0) 
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Os_ 

47T 


—( n g C A + riqC f ) + — (^——n g (3o — —n q CF + 2A(/j, 2 


d (0) 


(G4) 


where we used the explicit expression for Zq derived in Sec. IVIIIBI One can easily check that this exactly reproduces 
the IR-divergent parts of the QCD amplitudes. For example, for ggqq, we have 


"^2 (2 Ca + 2CV) + — o — 3 Cp + 2,A ggq q{g")^ C+}_^(pi,p2',P3,P4,) — 


/ 4 w ( 1 + > 2 -; 3 + 4q 
A«(2-, 1 + ;3+ 4r) | . (G5) 

VB&(1+ 2-;3+ 4T), 


Hence, the IR divergences in (Of)(b(7( 0 ) and A ^ cancel each other and do not enter in dw, as must be the case. 
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